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Abstract 

This article aims at an integrated formulation of BSE's for 2- and 3-quark 
hadrons under the Markov- Yukawa Transversality Principle (MYTP) which pro- 
vides a deep interconnection between the 3D and 4D BSE forms, and hence offers 
a unified treatment of 3D spectroscopy with 4D quark-loop integrals for hadronic 
transitions. For the actual dynamics, an NJL-type realization of DB\S is achieved 
via the interplay of Bethe-Salpeter (BSE) and Schwinger-Dyson (SDE) equations, 
which are simultaneously derivable from a chiral Lagrangian with a gluonic (Vector- 
exchange) 4-fermion interaction of 'current' uds quarks, specifically addressing the 
non-perturbative regime. A prior critique of the literature on various aspects of the 
non-perturbative QCD problem, on the basis of some standard criteria, helps con- 
verge on a BSE-SDE framework with a 3D-4D interconnection based on MYTP. 
This framework is then employed for a systematic self-contained presentation of 2- 
and 3-quark dynamics on the lines of MYTP-governed DBxS, with enough calcu- 
lational details illustrating the techniques involved. Specific topics include: 3D-4D 
interconnection of qq and qqq wave functions by Green's Function methods; pion 
form factor; 3-hadron form factors with unequal mass loops; SU(2) mass splittings; 
Vacuum condensates (direct and induced); Complex H.O. techniques and SO(2, 1) 
algebras for Baryon Spectroscopy; plus others. 
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1 Introduction: QCD-Type Confinement Models 



One of the biggest challenges in physics to-day is a viable theory of strong interactions 
for which QCD is the leading candidate. Unfortunately, despite many of its extremely 
attractive features, this theory is not yet available in a sufficiently tractable form so as 
to appeal instantly to all its practitioners in as universal a manner as, e.g., in QED. The 
bone of contention in this regard is the non-perturbative sector of QCD which shows 
up as the phenomenon of 'confinement' at low and moderate energies: As yet there is 
no visible evidence of a sort of minimum consensus on a common dynamical framework 
to incorporate this physical effect in QCD applications in the strong interaction sector 
in a sufficiently convincing yet doable manner. As a result, there exist a multiplicity of 
approaches which, while incorporating the QCD ideas in varying degrees of sophistication, 
nevertheless often need to resort to additional parametric assumptions to calculate various 
low energy hadronic properties. Some of the principal approaches are: 
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Bag models [1]; QCD-sum rules [2]; later adaptations [3] of the Nambu-Jona-lasino 
model [4]; QCD bosonization approaches [5]; Instanton methods [6]; Vacuum self-dual 
gluon fields [7] ; Quark confinement models [8] ; Schwinger-Dyson and Bethe-Salpeter mod- 
els [9]; adaptation of BSE-SDE to a 3D-4D hybrid form [10-11] in a 'two-tier' fashion (to 
incorporate the spectroscopy sector). To this list one should also add QCD-motivated 
'quarkonia' models which is one of the oldest types in existence, and whose state of the 
art may be found in a fairly recent collection [12]. Last not least, one must pay homage to 
"Lattice QCD" which addresses confinement at a more fundamental level, and has grown 
into a self-contained field of study by itself. However its philosophy and methodology 
have so little in common with the less ambitious approaches listed above [1-12], that it 
does not fall within the scope of the present study. 



1.1 A Short Critique of Models [1-12] 



If the pretence of implementing confinement through an exact solution of the QCD equa- 
tions of motion is given up in favour of an 'effective confinement' programme, the central 
issue boils down to the extent to which the same can be formulated in a manner which is 
both physically convincing as well as mathematically tractable enough to warrant wide- 
ranging applications, all the way from low-energy spectroscopy [12] to deep inelastic pro- 
cesses amenable to perturbative QCD [2]. Such a philosophy is reminiscent of Bethe's 
"Second Principle Theory" for effective nucleon-nucleon interactions, now reborn at the 
level of quark-quark interactions, with confinement addressed in a semi- empirical manner 
which incorporates the main features of QCD structure. Unlike Lattice-QCD, such pro- 
grammes are not meant to address confinement directly, but rather to take its role more 
or less for granted in anticipation of future developments. It is from this angle that most 
of the approaches listed above may be viewed, from putting in the QCD feature by hand 
[1,8-12], to a conscious effort to 'derive' its content more explicitly [5-7]. Of the last group, 
the model that comes closest to tackling confinement is perhaps [7], but its methodology 
has understandably more limitations for wider applicational purposes. A complementary 
role is that of [5] which is characterized by a 'chiral perturbation' approach to the mecha- 
nism of formation of hadronic states in QCD, and in the process gives rise to an effective 
chiral Lagrangian for low energy hadron physics. However the perturbative expansion in 
the momenta, deemed small in the low momentum limit, robs such a Lagrangian of a 
vital property: its capacity to predict the bound (confined) states of hadrons in the low 
momentum regime, due to the lack of a 'closed form' approach. (A closed form approach 
is best exhibited by some sort of form factors characterized by a confinement scale, a 
feature that gets lost in any expansion in the momenta). 

A 'two-tier' 3D-4D BSE approach like [10-11] is meant for a 'more conscious' incor- 
poration of the spectroscopy sector, i.e., an explicit recognition of the fact, often ignored 
in the more usual formulations of BSE-cum-SDE methods [9] , that the observed hadronic 
spectra are 0(3)-like [13], while a literal BSE formulation in euclidean form, with a stan- 
dard 4D support to the kernel leads to 0(4)-like spectra [14]. In this respect, the two-tier 
strategy [11] invokes the Markov- Yukawa Transversality Principle (MYTP) [15] wherein 
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the quark-quark interaction is in a hyperplane which is transverse to the 4-momentum 
P M of the composite hadron, so that the modified BSE has a (covariant) 3D support to 
its kernel. This feature in turn leads to an exact 3D reduction of the BSE from its 4D 
form, and an equally exact reconstruction of the 4D wave function in terms of 3D ingre- 
dients [16], thus implying an exact interconnection between the 3D and 4D BSE forms 
[16]. (A parallel formulation of MYTP [15] by the Pervushin group [17] gave rise to the 
3D reduction from the 4D BSE form, but the inverse connection (from 3D to 4D) was 
missing in their paper [17]). The 3D BSE form makes contact with the observed 0(3)- 
like spectra [13], while the reconstructed 4D BSE form provides a natural language for 
evaluating transition amplitudes via quark loops [16,11]. 

A 3D BSE form has its own logical basis which receives support from several indepen- 
dent angles in view of its crucial role in the understanding of physical processes in general 
and the theme of the present article in particular (see below). We list three supporting 
themes which have been developed over the decades from entirely different premises, all 
converging to a basically 3D picture for the effective qq/qq interaction: 

i) It gives a physical meaning to the interaction of the quark constituents in their 
respective mass shells, consistently with the tenets of local field theory [18]; ii) It arises 
from the concept of instantaneous interaction [17,16] among the quark constituents in 
accordance with the Markov- Yukawa picture [15] of transversality to the composite 4- 
momentum P M ; iii) It is the only structure of the BSE kernel which makes this equation 
compatible with a pair of Dirac equations for two particles under their mutual interaction 
[19]. 



1.2 Bethe's 'Second Principle' Criteria for Model Selection 



Next we consider certain guiding principles (criteria) which form the basis of this study, be- 
fore identifying the specific model/models for a more detailed and self-contained exposure. 
A clue is provided by the observation that most of these approaches [1-12], irrespective 
of their individual theoretical premises, have a common characteristic: Applicability to 
hadronic processes viewed as quark composites, limited only by their individual predic- 
tive powers, while a deeper understanding of the underlying models themselves is best 
left to future investigations. Perhaps the painfully slow progress of Lattice QCD results 
gives an inkling of this scenario: The formidable dimensions of the quark-gluon strong 
interaction physics leaves little alternative to the respective practitioners of QCD but to 
settle for a less ambitious approach to the problem. Thus the different models [1-12] are 
best regarded as alternative strategies, each with its own methodology and parametric 
limitations, aimed at selected sectors of hadron physics that are suited to their structural 
budgets before the 'final' theory unfolds itself, leaving their successes or otherwise to be 
judged in the interim by the depth and range of their respective predictions vis-a-vis the 
data. This is once again Bethe's 'second principle' philosophy in retrospect, which would 
presumably continue to operate perhaps as long as QCD remains a partially solved theory. 
Within this restricted philosophy, some obvious criteria for theme selection from among 
the available candidates [1-12] could be the following partial shopping list (not mutually 
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exclusive) : 

A) Maximum number of mutually compatible as well as time-tested ideas that can be 
extracted from out of [1-12] in the sense of an "HCF" ; B) Close proximity to QCD as the 
ideal theory which stands on the three pillars of Lorentz-, Gauge- and Chiral- invariance; 
C) Formal capacity to address several sectors of physics simultaneously, all the way from 
hadron spectroscopy to quark loop integrals for different types of transition amplitudes 
within a common dynamical framework; D) Sufficient flexibility of the core dynamical 
framework to permit smooth on-line incorporation of mutually compatible ideas, like 
Markov- Yukawa transversality MYTP [15], and Dynamical Breaking of Chiral Symmetry 
{DB\S) [4], and other similar principles if need be, without causing any structural (or 
parametric) damage to the basic framework itself; E) Natural capacity of the conceptual 
premises to include both 2- and 3-quark hadrons within a common dynamical framework, 
to give concrete shape to the (widely accepted) principle of meson-baryon duality; F) 
A built-in microcausality in the dynamical framework which takes in its stride sensitive 
items like the structure of the vacuum in strong interaction physics, without the need for 
fresh ansatze/parametrizations. 

Although these criteria are neither exhaustive nor mutually exclusive, their collective 
effect is nevertheless focussed enough to eliminate many prospective candidates in favour 
of a chosen few that would survive the tests (A)-(F) for a reasonably self-contained account 
of strong interaction physics within the tenets of the 'Second Principle' philosophy. Thus 
the too simplistic premises of Quarkonia models [11] which played a crucial role in the early 
stages of QCD- motivated investigations, would not stand the tests of (B),(C) and (D). 
Bag models [1] which had also played a similar role in the early phases, would not qualify 
under (B), (C) and (F). QCD sum rules [2] represent perhaps one of the most successful 
applications of perturbative QCD by relating the high energy quark-gluon sector to the 
low energy hadron sector through the principle of an FESR-like duality discovered in the 
Sixties [20]. Even to-day it is extensively used for many hadronic investigations. Yet it 
fails on count (C) mainly because of its failure to satisfy (F): The lack of microcausality 
in this model can be traced to the 'matching condition' between the quark-level and 
hadron-level amplitudes whose solution is far from unique. Indeed, while the borelisation 
technique suffices for the prediction of ground state hadron masses, that very mechanism 
also causes it to lose information on the spectra of the excited states, thus reducing its 
predictability on this vital (low energy) front. 

As for models [5-8], the Quark Confinement model [8] lacks enough microcausality - 
condition (F) - which shows up through a relatively poor satisfaction of (C), since the 
spectroscopy sector is badly neglected. QCD bosonization methods [5] have made very 
impressive strides in respect of transition amplitudes through the powerful technique of 
construction of effective Lagrangians in terms of the hadronic fields, which make them 
ideally suited for 'tree-diagrams'. However such structures hide from view the composite 
character of the hadrons in the 'soft' QCD regime, which is best exhibited in 'closed' form 
via quark-hadron form factors - the vehicle for sensitivity to various non-perturbative 
features of the theory. For the baryon dynamics too, the inadequacy of the formalism [5] 
to depict correctly the 3-quark form factors shows up through its excessive dependence on 
the quark-diquark description with a rigid diquark structure which results in an inevitable 
loss of information on the true 3-quark content of the baryon. Finally, while [6,7] satisfy 
the conditions (B),(C),(F) on separate counts, there is not enough published evidence of 
support from the other 3 quarters (A),(D),(E). Hence while their comparison with others is 
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useful for a comparative discussion of different models, their claims to a primary 'theme' 
status fall rather short of a good starting point. This leaves the BSE-SDE framework 
[9-11] for a more detailed scrutiny to follow. 



1.3 BSE-SDE Framework:D£xS and 'Soft'-QCD 



The last group [9-11] is characterized by an interplay of BSE and SDE, both derivable from 
a suitably chosem 4-fermion Lagrangian as input. It has a very wide canvas and is fully 
attuned to Bethe's Second Principle Theory (BSPT for short). Its general framework 
equips it with arms to meet most of the conditions (A)-(F), ranging from flexibility to 
wide-ranging predictivity, thus lending it some credibility within the broad premises of 
'BSPT'. In particular, its natural roots in field theory endow it with standard features like 
dynamical breaking of chiral symmetry (DBxS) via the non-trivial solution of the SDE, 
thus giving it the powers to subsume the contents of the NJL-model [4]. Indeed, soon after 
the discovery of the NJL model [4], a field-theoretic understanding of its underlying idea 
was achieved in the form of a generalized DB\S in the QED domain via the non-trivial 
solution of the SDE [21]. And after the advent of QCD [22], the same feature showed up 
through the solution of the BSE for qq interaction via one-gluon-exchange [3]. This is a 
typical non- perturbative effect, although it does not cover all its aspects. 

Subsequently the concept of DB\S was generalized to show that this feature is shared 
by any extended vector-type 4-fermion coupling [23-24] which preserves the chiral sym- 
metry of the Lagrangian but the same gets broken dynamically through the non-trivial 
solution of the SDE, derivable from such a Lagrangian. Indeed the sheer generality of a 
Lagrangian- based BSE-cum-SDE framework, by virtue of its firm roots in field theory, 
gives it a strong mandate, in terms of both predictivity and flexibility, to accommodate ad- 
ditional principles like Markov- Yukawa Transversality [15] (MYTP) while staying within 
a basically Lagrangian framework, so that criteria (A)-(F) are still satisfied. 

In particular, a basic proximity to QCD is ensured through a vector- type interaction 
(condition (B)) [10,23], which while maintaining the correct o.g.e. structure in the per- 
turbative region, may be fine-tuned to give any desired structure in the infrared domain 
as well. The latter part is admittedly empirical, but captures a good deal of physics in the 
non-perturbative domain while retaining a broad QCD orientation, and hence does not 
rule out a deeper understanding of the infrared part of the gluon propagator within the 
same framework. More importantly, the non-trivial solution of the SDE corresponding to 
this generalized gluon propagator [11] gives rise to a dynamical mass function m(p) [11] 
as a result of DBxS, even while the input Lagrangian has chiral invariance due to the 
vector-type 4-fermion interaction [23,24] between almost massless u — d quarks. These 
considerations further strengthen the case of a Lagrangian-based BSE-SDE framework for 
a theme choice. 
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1.4 3D-4D BSE: From Spectra To Loop Integrals 

Now the canvas of a (second-principle) BSE-SDE framework is broad enough to accom- 
modate a whole class of approaches, and facilitates further fine-tuning in response to the 
needs of major experimental findings such as the observed 0(3)-like spectra [13], which 
essentially amounts to treating the time-like momenta separately from the space-like ones, 
as has long been known since the classic work of Feynman et al [25]. In this regard, the 
MYTP constraint [15] seems to fit this bill, by imparting a 3D support to the pair- 
wise BSE kernel [15-17], an ansatz which can be motivated from several different angles 
[16,18,19]. As to the 'soft' non-perturbative part of the gluonic propagator, it still remains 
empirical since orthodox QCD theory does not yet provide a closed form representation 
out of the infinite chain of equations that connect the successively higher order Green's 
functions in the standard fashion [26], thus necessitating parametric representations [27]. 
Parametrization is also compatible with MYTP [15] (see [11]), wherein the key constants 
are attuned to the hadron spectra of both 2- [28] and 3- body [29] types, within a common 
framework. The 3D support ansatz of MYTP [15] in turn gives a characteristic 'two-tier' 
[16] structure to the entire BSE formalism, wherein the first stage (3D BSE) addresses the 
meson [28] and baryon spectra [29], while the reconstructed 4D wave (vertex) functions 
[16] fit in naturally with the Feynman language of 4D quark loop diagrams for various 
types of transition amplitudes [11,30-33] in a unified fashion. 

A BSE-SDE formulation [9] represents a 4D field-theoretic generalization of 'potential 
models' [12], and is thus equipped to deal with a wider network of processes (e.g., high 
energy processes) not accessible to potential models [12]. In this way, the BSE-SDE 
approach occupies an intermediate position, sharing the off-shell feature with potential 
models [12], as well as the high energy flavour of QCD-SR [2], but its dynamical spirit is 
much nearer to [12] than to [2]). Indeed the role of the 'potential' [12] is played by the 
generalized 4-fermion kernel [23] (which is a paraphrase for the non-perturbative gluon 
propagator [11]). The 4D feature of BSE-SDE gives this framework a ready access to high 
energy amplitudes, as in QCD-SR [2] as well as in other models [5-8], while its 'off-shell' 
feature gives it a natural access to hadronic spectra [13], in company with (potential- 
oriented) quarkonia models [12] . (In contrast, certain models [2,5-8] do not have a 
basic infrastructure to address spectroscopy). Now with this twin feature of off-shellness 
and Lorentz-covariance, the BSE-SDE framework formally overcomes the shortcomings of 
'potential' models [23] in obtaining numerically 'correct' values for the various condensates 
which are employed as inputs in QCD-SR calculations [2]. This was indeed confirmed by 
a later derivation of similar results [30] in terms of a Lorentz-covariant formulation [11] of 
the BSE-SDE framework, which showed that vacuum condensates are calculable within 
a spectroscopy-rooted [28,29] framework. 

1.5 Off-Shellness in BSE: Parametric Links with QCD-SR 

The calculations [23,30] raise the interesting question of the possibility of a basic con- 
nection among the input parameters of different models, although conceived within very 
different premises. Thus in QCD-SR [2], the 'free' parameters of the theory are the 
condensates themselves as input, while in the BSE-cum-SDE methods [9-11], the corre- 
sponding parameters are contained in the input structure of the infrared part of the gluon 
propagator [11]. Now since the condensate parameters of QCD-SR [2] are explicitly calcu- 
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lable in BSE-SDE models by quark loop techniques [23, 30], using the gluonic parameters 
[11,30], this result at least settles the issue of a one-way connection: from BSE [11] to 
QCD-SR [2]. 

To pursue this question a bit further, let us compare the features of potential mod- 
els (of which BSE is a 4D generalization) with those of QCD-SR: Potential models are 
characterized by 'off-shell' features, whose parameters (corresponding to given 'potential' 
forms) are primarily attuned to low energy spectroscopy, so that their predictions tend 
to work upwards on the energy scale, starting from the low energy end. QCD sum rules 
on the other hand are attuned to the perturbative QCD regime, so that their predictions 
tend to work downwards on the energy scale, starting from the high energy end. The 
'softness' aspects of QCD-SR are typically simulated via the Wilson OPE expansion in 
inverse powers of 4- momentum Q 2 where the 'twist' terms of successively higher dimen- 
sions are symbolized by the corresponding 'vacuum condensates' which are thus the free 
parameters of the theory. Therefore prima facie it appears that the two methods are 
largely complementary to each other. The former, by virtue of its low energy/off-shell 
emphasis, is particularly successful on the spectroscopic front, but its techniques do not 
find easy access to transition amplitudes due to inadequate treatment of the high energy 
front (lack of covariance). The latter (QCD-SR) on the other hand, is ideally suited to the 
high energy regime, but does not find ready access to areas involving soft QCD physics, 
especially the spectroscopic regime. This is at least partly attributable to the method- 
ology of QCD-SR [2] which makes use of the 'quark-hadron duality' for 'matching' the 
respective amplitudes [2,20]: Because of the relatively 'macroscopic' nature of the 'match- 
ing' which is effected with a 'Borelization' technique [2], the predictions are reliable only 
for the hadronic ground states, but do not readily extend to the spectra of excited states. 

Due to the complementary nature of the two descriptions, it is not generally easy to 
relate the parameters of one to those of the other. However, the off-shellness feature of 
potential models gives it access to information on the interaction of the quark pair with 
the environment; in particular they possess signatures on the structure of the degenerate 
vacuum in the form of 'vacuum condensates'. That the vacuum condensates of QCD-SR 
[2] can be be expressed in terms of potential [23] and BSE [30] models is a reflection of their 
crucial 'off-shellness' property. As to the converse question, there is no published evidence 
of a corresponding exercise in the opposite direction viz., a derivation of the parameters 
of the BSE kernel/gluon propagator in terms of the various vacuum condensates that 
characterize QCD-SR [2]. A possible reason may lie in the role of microcausality (condition 
(F)) which is well satisfied by potential models, but perhaps not by QCD-SR [2]. Thus 
it would appear that 'microcausality' which underlies the 'off-shellness' feature of the 
'potential models' enhances their predictive powers vis-a-vis those which do not possess 
this crucial property. 

Now the off-shell structures of all 'potential-oriented' models [9-12] have a fairly direct 
connection with the 'spectral' predictions, unlike other types of confinement models [2,5- 
8] , which do not permit such predictions in an equally natural way. And for 3-quark states 
[29], the dichotomy seems to be even sharper, inasmuch as there is a strong tendency in 
the literature to simplify the 3-quark systems as quark-diquark systems [2,5-8], thus partly 
"freezing" some genuine 3-body d.o.f.'s and causing a loss of information on the spectra 
of L-excited states. 

The off-shell characteristics of the BSE-SDE framework [9-12] are perhaps the most 
important single feature responsible for extending their predictive powers all the way from 
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3D spectra to 4D transition amplitudes (of diverse types) via 4D quark loop integrals, 
under one broad canvas. The key to this capacity lies in the vehicle of the BS wave 
(vertex) function which has at its command the entire 'off-shell' information noted above. 
Here it is important to stress that this wave function is a genuine solution of the BS 
dynamics [11], so that it leaves no scope for any free parametrization beyond what is 
already contained in the (input) gluon propagator. (Potential models [12] also have this 
capacity in principle, but their 3D structure does not allow full play to the 'loop' aspect). 

1.6 Markov- Yukawa Transversality on the Null Plane 

Covariant Instaneity Ansatz (CIA) on the BSE [16] is not the only form of invoking 
MYTP to achieve an exact interconnection between the 3D and 4D structures of BSE. 
As will be found later (see Sec. 4 for details), the CIA which makes use of the local cm. 
frame of the qq composite, has a disadvantage: The 4D loop integrals are ill-defined due 
to the presence of time-like momentum components in the exponential/gaussian factors 
(associated with the vertex functions) caused by a 'Lorentz-mismatch' among the rest- 
frames of the participating hadrons. This is especially so for triangle loops and above, 
such as the pion form factor, while 2-quark loops [32] just escape this pathology. This 
problem is probably absent if the null-plane ansatz (NPA) is invoked, as found in an 
earlier study of 4D triangle loop integrals [33], except for possible problems of covariance 
[34]. The CIA approach [16] which makes use of the TP [15], was an attempt to rectify 
the Lorentz covariance defect, but the presence of time-like components in the gaussian 
factors inside triangle loop integrals [31] impeded further progress on CIA lines. 

Is it possible to enjoy the best of both the worlds, i.e., ensure a formal covariance 
without having to encounter the time-like components in the gaussian wave functions 
inside the 4D loop integrals ? Indeed the problem boils down to a covariant formulation 
of the null-plane approach. Now the null-plane approach (NPA) itself has a long history 
[35], and it is not in the scope of this article to dwell on this vast subject as such. Instead 
our concern is limited to the covariance aspects of NPA, a subject which is of relatively 
recent origin [18,36-38]. However in all these approaches [38], the primary concern has 
been with the NP-dynamics in 3D form only, as in the other familiar 3D BSE approaches 
[39] over the decades. On the other hand, the aspect of NPA which is of primary concern 
for this article, is on the possibility of invoking MYTP for achieving a 3D-4D BSE 
interconnection on the covariant Null Plane, on similar lines to Covariant Instantaneity 
(CIA) for the pairwise interaction [16]. Now it seems that a certain practical form of the 
null-plane formalism [33] had all along enjoyed both 3D-4D interconnection and a sort of 
'pedagogical covariance '(albeit implicitly) [40]. This basic feature can be given a formal 
shape by merely extending the Transversality Principle [15] from the covariant rest frame 
of the (hadron) composite [16], to a covariantly defined null-plane (NP) [41]. Because of 
its obvious relevance, the subject of 3D-4D interlinkage on the covariant Null-Plane [41] 
will be covered in Sec. (4.2), with a parallel CIA treatment in Sec. (4.1). 
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1.7 Scope of the Article: Outline of Contents 



We now focus on a BSE-cum-SDE form of dynamics derivable from a chirally invariant 
Lagrangian with an effective gluon-exchange-like interaction (pairwise), as the central 
theme of this study for a reasonably self-contained presentation, under the further con- 
straint of Markov- Yukawa Transversality Principle (MYTP). The emphasis is on a peda- 
gogical perspective on the problem of effective color confinement, converging on a vector 
exchange mediated Lagrangian whose chiral symmetry gets broken dynamically, after giv- 
ing a bird's eye view of the main approaches to effective confinement [1-11]. Indeed the 
(DBxS) theme, although originating from the NJL-model [4] for contact pairwise inter- 
action, admits a simple generalization to a (space-time extended) vector exchange qq / qq 
interaction which exhibit chiral symmetry at the input Lagrangian level, but get broken 
dynamically via the solution of the Schwinger-Dyson Equation (SDE) [23-27]. A more 
explicit QCD motivation must be achieved by hand, e.g., identification of the pairwise 
interaction with the entire gluon propagator (perturbative and non-perturbative [28a]), 
which in turn has several desirable consequences, such as the color effect which ensures 
that the strength of the qq force is half that of qq, within a common parametrization. 

The second item of emphasis concerns the remarkable facility of an exact interconnec- 
tion between the 3D and 4D BSE forms [16], that is provided by MYTP, a facility that 
other 3D approaches to BSE [39], or (basically 3D) Null-Plane approaches [35-38], do not 
seem to possess. This property allows the exposiition of the BSE-cum-SDE techniques in 
a very simple way, so as to provide the reader with a quick working knowledge of their 
applications to a wide class of problems which may be broadly classified in a two-tier form: 
A) Mass spectra; B) Quark-Loop diagrams. Such a division is natural since investigations 
of types (A) and (B) are mainly governed by the 3D and 4D aspects of the BSE respec- 
tively. Therefore after an introductory phase on the general BSE-SDE formulation, an 
early specialization to its MYTP-governed 3D-4D form (from Sec. 4 onwards) will form 
the basis for this (application oriented) article. 

A third item of emphasis is on the second stage of the 3D-4D BSE framework, viz., 
techniques of 4D quark-loop amplitudes, with a comparative study of CIA [16] vs CNPA 
[41], to bring out their relative (strong/weak) features. 

This article has been built on the infrastructure of one with a similar theme [40] writ- 
ten about a decade ago; it incorporates major advances through the present decade on 
the 3D-4D BSE front [42a-b], viz., the Covariant Instaneity Ansatz (CIA) and its more 
recent Null-Plane counterpart (CNPA) [41], both under the umbrella of MYTP [15]. The 
background of ref.[40] will be freely used, but the details on (3D) spectra on which CIA 
[28] and CNPA [41] have similar predictions, will now be omitted, except for drawing 
attention to their structural similarities. Instead more attention will be paid to the struc- 
ture of 4D quark loop integrals of selected types to bring out the applicational potential 
of this MYTP-governed formalism [41-2]. These types include i) certain hadronic form 
factors built out of triangle loops; ii) typical self-energy problems dealing with SU(2)- 
mass splittings among hadrons; iii) vacuum condensates which are inputs in QCD-SR [2], 
but calculable in the 3D-4D BSE-SDE formalism [30]. 

While giving the details of this article, we repeat at the outset that, except for the 
contents of Sections 1-3, it is not intended as a conventional 'review' of the BSE-SDE 
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framework such as [9]. Nor are conventional 3D BSE approaches [39], or the conventional 
NPA formalisms [35-38] the subjects of our detailed description. Aspects of contact NJL- 
type 4- and 6-fermion couplings (often employed in the'nuclear' field), are also not of 
interest here. 

As to the actual details, the Table of Contents, preceding the Introduction (Section 
1) gives a fair cross section of the included items: Sect. 2 gives a panoramic view of the 
NJL-Model [4] and its aftermath. Sect. 3 gives a general derivation of BSE and SDE 
in an interlinked fashion, with a gluon-like (Vector-exchange) propagator whose mass 
function m(p) stems via DB\S from a spatially extended 4-fermion interaction in the 
input Lagrangian. With this general background of SDE-BSE as well as of DBxS, the 
rest (Sect. 4-11) deals with different facets of the 3D-4D BSE-SDE framework under the 
Markov- Yukawa Transversality Principle [15] at two distinct levels of operation, viz., CIA 
[16, 17] which has been around for some time, and CNPA [41] which is formally a new 
proposal, although in effective (practical) use for quite some time [33,40]. 

Of the subsequent Sections, Sect. 4 collects the background for interlinked 3D-4D 
BSE techniques for qq hadrons. For the fermionic BSE, we have preferred to stick to 
its 'Gordon-reduced' version [lOa-b] adapted to the off-shell constituents [10a]. This is 
a conscious departure [10b] from the standard BSE-form [26] to make the BSE more 
tractable for wider applications, as in other BSE approaches [lOc-d], and does not violate 
the 'Bethe Second Principle' spirit, since the input 4-fermion coupling is an effective 
description of the pairwise interaction). 

Sects. (5-8) deal with some selected applications of triangle loops (form factors), two- 
loops (self-energy), and one-loop (vacuum condensates) techniques respectively. These 
include, among other things, a technique to include QED gauge insertions in arbitrary 
momentum-dependent vertex functions for the e.m. self energy and form factors. Wher- 
ever possible, a parallel treatment is provided for CIA and CNPA for a comparative view 
of the two distinct MYTP-governed BSE formalisms, but some technical problems with 
CIA [16] often lead to a preference for CNPA. Some calculational details on the form 
factor plus normalization are given in Appendix A. 

Sect. 6 gives a general method for triangular quark-loop integrals applicable to a large 
class of transition amplitudes for 3-hadron coupling [31], to bring out a major simplifying 
feature of the resulting structure arising out of a 'cancellation' mechanism between the 4D 
quark propagators and the 3D D-functions in the hadron-quark vertices of the two-tier 
BS formalism [16]. This prevents free propagation of quarks by eliminating the Landau- 
Cutkowsky (overlapping) singularities [16,31]. 

Sects. 7,8 give results for self-energy diagrams [32] and of vacuum condensates [11,30], 
requiring two and one .Sp-functions respectively. The self-energy calculations in Sect. 7 are 
illustrated with SU(2) mass splittings of pseudoscalar mesons [32b]. A general method 
to deal with QED gauge corrections to the e.m. mass differences is outlined in Appendix 
B. For the vacuum condensates [11,30], Sect. 8 offers a new gauge invariant technique for 
loop integrations, on the lines of Schwinger [43]. We reiterate that such predictions are 
intimately linked with spectroscopy via the infrared structure of the gluon propagator 
[11,30]. 

The third part Sects. 9-11 deals with the BSE formalism for a 3-quark baryon, with 
emphasis on the qqq structure taking into view that in most approaches, including other 
BSE models [9b], the dynamical treatment has often relied heavily on the quark-diquark 
approximation [5b, 8b, 9b, 11], which amounts to a "freezing" of the 3-body degrees of 
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freedom. It has also been recognized in the literature that with a 3-body BSE treatment to 
the baryon, there are some technical problems associated with the status of the spectator 
[44]. In the Two-tier BSE model this problem has been regularly addressed at various 
stages of its development [10,40,29]. The 3-quark dynamics is described in 3 Sections 
(9-11) as follows. 

Sect. 9: A panoramic view of the baryon dynamics as a general 3-body problem with 
full permutation symmetries [45] in all the relevant d.o.f.'s incorporated; a detailed cor- 
respondence with the quark-diquark model; Complex HO techniques for the qqq problem 
[46]; problems of 3D reduction and 4D reconstruction for qqq BSE [47]; and fermionic 
BSE with gluonic interactions in pairs [29]. 

Sect. 10: Green's function techniques for 3D reduction of the BSE, and reconstruction 
of the 4D qqq wave function [47] ; see Table of Contents. 

Sect. 11: A summary of the relativistic fermionic qqq BSE with the same gluonic 
propagator as employed for the qq problem; the 3D reduction [29] of the qqq BSE is on 
closely parallel lines to the two-body case [28]. The derivation of an explicit mass formula 
is greatly facilitated by taking a complex HO basis [46]. However loop techniques for 
baryonic amplitudes are not included for explicit presentation. 

2 NJL Model: Recent Developments (Nambu) 

The precursor of the NJL-model [4] was the 'Nambu-Goldstone' picture of the pion as a 
zero mass particle arising from the chiral non-invariance of the vacuum [48]. This view of 
the pion received quantitative shape at the hands of Gell-Mann and Levy [49a] who started 
with a SU(2) x SU(2) symmetry of the Lagrangian (termed SU(2) cr-model) involving an 
1 = 1 pseudoscalar n and an I = scalar a field. Due to spontaneous symmetry breaking 
of the vacuum, the a-field shifts to a minimum < a >= —f w ^ 0, while the pion field 
remains unshifted (< n >= 0) and stays at zero mass. 

The Gell-Mann Levy a-model set the stage for modern chiral theories, stimulated by 
an important paper due to Skyrme [49b], to describe pseudoscalar mesons and baryons 
through a solitonic picture wherein baryons are generated as bound states of weakly 
interacting mesons. These models were developed in the QCD context wherein, in the 
large N limit, QCD becomes equivalent to a non-linear meson theory. The underlying 
logic is that although the QCD Lagrangian has chiral symmetry for massless quarks, this 
symmetry is spontaneously broken, giving rise to massless pions, etc. These methods give 
rise to effective Lagrangian descriptions at the tree level [5], but will not concern us any 
further in this article. 

The NJL-model [4] on the other hand, which is the very raison d'etre of this article, is 
characterized by chiral symmetry breaking in a dynamical fashion, and allows a formally 
composite structure of the pion, in company with other hadrons. This is a distinct ad- 
vance over the elementary field picture of the pion [48-49] , and facilitates a more natural 
understanding of many of its observational properties (form factor, L-excited states, etc) . 
A short summary of the NJL model follows. 
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2.1 Outline of NJL Model 

The NJL Lagrangian [4] may be written in two different ways: 

L NJL = L + L t = (L + L s ) + (Li - L s ) = L' + L{; (2.1) 

where Lq = —ip'j.dif} and L{ (see below) are chirally invariant, but L s = —mipip which 
stands for the observed fermion, is not, and represents the symmetry breaking effect. The 
interaction term Lj is given by 

U = g [(^) 2 - {^f] = -g [(^l^) 2 - (^7m75^) 2 ]/2 (2.2) 

The rearrangement in (2.1) is meant to diagonalize L' Q , and treat L\ as a perturbation; 
this implies a redefinition of the vacuum by introducing a complete set of 'quasi-particle' 
states which are eigen states of L' . The L s is now determined from the requirement that 
L- shall not yield additional self-energy effects. This gives the standard Schwinger- Dyson 
Equation (SDE) for m in terms of the loop self energy 

m = S| (i7 . p+m=0) = -8img f (2n)~ 4 d 4 p f^ P '^ . (2.3) 

J m z + p 2 — le 

where F(p, A) is a cut-off factor. The trivial solution m = corresponds to the usual 
chiral invariant vacuum characteristic of perturbation theory. The non-trivial solution 
m = m NJL is found from 

(2tt) 4 = -8ig J d A p{m 2 + p 2 - ie) _1 F(p, A) (2.4) 

in terms of g and A. It is alos called the 'gap' equation, and is based on a shifted vacuum 
Q m which is chiral non- invariant. With a fixed Lorentz-invariant cut-off A in Euclidean 
space (and F — 1), eq.(2.4) reduces to 

27i 2 /g = A 2 - m 2 ln(A 2 /m 2 + 1); < 27r 2 g Q - 1 A- 2 < 1 (2.5) 

The two vacuua Q$ (chiral invariant) and Q m (non-invariant) are fully orthogonal to 
each other, and correspond to two different worlds. Q m , with the lower energy, is the 
true ground state. The chirality operator defined as x — 1 4>l4lbi>d! i commutes with the 
original hamiltonian Hq with vacuum Qq, but not with H m with vacuum Q m . However 
X has no matrix elements connecting the two worlds Qq and Q m , a sort of superselection 
rule. Now the following paradox arises: The ^-conservation in the Qq basis implies the 
existence of a conserved current 

j M 5 = #7a»75^; did/* = (2.6) 
On the other hand, for a massive Dirac particle in the fl m basis 

W7„75V0 = 2mV>75^ + (2.7) 

To reconcile these two statements, the x-current operator between physical states suffers 
radiative corrections w.r.t. the simple term 27^75 so that, on grounds of Lorentz invariance 

< P'\J^\P >= F(k 2 )u(p')[t llll5 + ^^]u(p); (k = p- p') (2.8) 
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Thus the real fermion (quark) is not a point particle since its \- current has an anomalous 
7 5 term. This in turn implies a pole at k 2 = for the 75 term, corresponding to a zero 
mass pseudoscalar, whose natural identification is the pion. 

The pion which arises here as the lowest qq bound state, has clearly the nature of a 
collective excotation, thus also implying the existence of higher excitations in the same 
package, something which the elementary field model [48-49] could not provide. Indeed 
the BS amplitude \I> for the bound state composite is 

^(x,y)=<0\T(^(x)^(y))\B> (2.9) 

which is related to the vertex function T in momentum space as 

*(pi,p 2 ) = S F (q + P/2)T(q, P)S F (q - P/2) (2.10) 

where the individual quark momenta pi^ = P/2 ± g in terms of the total (P) and relative 
(q) 4-momenta. For the pseudoscalar state in question, the BSE for T(p 1 ,p 2 ), viz., 

(27r) 4 r(g, P) = 2^o75 / d 4 q'Tr[ l5 S F (q' + P/2)T(q', P)S F (q' - P/2)] (2.11) 

which for P M = 0, has a self-consistent solution T = C75, C being a constant, provided go 
satisfies eq.(2.4), which is just the 'gap' equation (SDE) for the mass m. 

This crucial result of the NJL model, which shows that in the chiral limit P M = 0, 
the BSE and the SDE are identical, is a direct consequence of the 75-invariance of the 
input Lagrangian. It also tells us that in the P M = limit, the qq vertex function and 
the quark mass function m have the same (constant) structure. The constancy of each 
is of course a consequence of the contact interaction, but the basic equality of these two 
quantities is also valid for an extended 4-fermion chirally invariant Lagrangian such as a 
vector mediated one [24]. 

The true significance of NJL was realized in the QCD context [22], through the study 
of non-perturbative solutions of SDE as a DB\S mechanism for more general 4-fermion 
couplings [3,8,23,24], including reformulations of the bag model [8a], and renormalization 
group equations [24d]. And it was eventually subsumed in the generalized BSE-SDE 
formalism [9-11], which is of course the subject of this review. 

2.2 BCS Mechanism, Mass Relations, SYSY, Etc 

During the last decade, Nambu [50] has abstracted the findings of a new symmetry from 
BCS-type theories of dynamical symmetry breaking (due to short-range attraction), re- 
sulting in a new vacuum state. The residual symmetry in question is a remarkably simple 
relation among the fermion mass mj and the composite boson (n, a) masses as low en- 
ergy modes in the new vacuum, viz., (M n : rrif : M a ) = (0:1:2). In more complex 
the fermion mass rrif and the (composite) boson masses (M 1? M 2 ) obey the generalized 
relation M\ + M| = 4mj. The low energy properties of the system can be represented 
by an effective Hamiltonian like in the a-model [49a] where the coupling constants are so 
related as to yield such mass relations automatically. 

Coming to the SUSY aspects, the essential thrust of Nambu's discovery [50b] is a hid- 
den SUSY in the BCS mechanism, manifesting via two physical scenarios: i) a cascading 
chain of symmetry breaking (tumbling); ii) a bootstrap mechanism in which the symmetry 
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sustains itself among a set of effective fields without the need to refer to a substructure. 
The main ideas are the following. 

A BCS mechanism has two energy scales: i) The high energy scale corresponds to the 
force responsible for the formation of Cooper pairs; its analogue in particle physics is the 
pion decay constant /„-. ii) The low energy scale is the pairing energy, and one of its 
manifestations is the quasi-fermion mass rrif which corresponds to the constituent quark 
mass m q = m^jL- More explicitly, there are both fermionic and bosonic excitations 
in the low energy scale: the quasi-fermion (rrif), the Goldstone boson (pion) and the 
Higgs (sigma) boson. In the simplest BCS (NJL) mechanism, their masses are in the 
ratio rrif : m n : m a = 1:0:2. This low energy picture can also be articulated by 
an effective Ginzburg-Landau-Gell Mann-Levy Hamiltonian involving these fermion and 
boson fields with Yukawa couplings and a Higgs potential. Their characteristic parameters 
are the 'high-energy' sigma-condensate c ~ /„-, and the 'low-energy' dimensionless Yukawa 
coupling constant G = rrif jc. To satisfy the mass ratio constraints, the Higgs self- coupling 
must be equal to G 2 . The non-relativistic analogue of the condensate c is \/N /2 where 
N is the density of states of the constituents at the fermi surface. The origin of the mass 
scales is a more dynamical question depending on the SUSY Hamiltonian structure, for 
a derivation of which the interested reader is referred to [50 a-d]. 

The physical scenario envisaged by Nambu for this broken SUSY structure is two-fold. 
The first is a cascading hierarchy of symmetry breaking ('tumbling') which in particle 
physics [51a] means something like the following. Suppose a symmetry breaking at a high 
energy level gives rise to a a-boson at the low energy scale. The latter, being a scalar, 
will induce attraction between the quasi-fermions, which in turn may generate a second 
generation symmetry-breaking, and so on. According to Nambu, a similar example of 
tumbling also exists in nuclear physics. Thus the a-boson, which is a fall-out of chiral 
symmetry breaking and quark-mass generation in the bulk of nuclear binding, also causes 
nuclear pairing which can be estimated quite accurately [50c]. 

The second scenario [50e] is the theoretical possibility of a (Chew-like) bootstrap, 
not at the hadron, but at the quark-lepton level, on the assumption that the t'Hooft 
self-consistency condition [51b] is satisfied between these two levels. This leads to the 
following bounds on the t-quark and Higgs masses: m t > 120GeV; m H > 200GeU. This 
and other details may be found in [50e]. 

3 Gauge Theoretic Formulation of SDE-BSE 

As seen in Sec. 2, the simple NJL-model [4] succinctly articulates the DB\S mechanism 
which gives rise to dynamical quark-mass generation on the one hand, and a Nambu- 
Goldstone [48] realization of the massless pion on the other. Another result is the formal 
identity of the mass-gap equation (SDE) with the homogeneous BSE for the vertex func- 
tion for a massless pseudoscalar qq composite. We are now in a position to pursue the 
same logic to give a formal theoretical basis to a gluon-exchange (vector)-like 4-fermion 
interaction (to simulate QCD effects) in the input Lagrangian by deriving from it an inter- 
linked BSE-SDE framawork [9-12] which is the backbone of this article. In this respect we 
shall skip an alternative non-perturbative treatment of the BCS-NJL pairing mechanism 
by the Bogoliubov-Valatin method [23] which is not easy to adapt to a Lorentz-invariant 
formulation. 
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3.1 Minimal Effective Action: SDE & BSE 



We outline a treatment due to Munczek [52] on the derivation of the equations of motion 
for composite fields. Consider an action functional 

S = J dx$(- r d - M)^ + ^p\(x) + h.c] ~\f J dxdy^ s G s (x - y)J s (x)J s (y) (3.1) 

where X(x) is an external source, G s is the propagator of the exchanged boson, and J s (x) 
=ijj(x)r s ijj(x) is the current function. This form is approximately derivable from the 
standard generating function for non-abelian QCD with r s = ry^A/2, when G s becomes 
the gluon propagator. The NJL-type contact interaction corresponds to G s = 5 4 (x — y), 
but the treatment is more generally valid for non-local interactions too. The standard 
approach is to introduce bilocal boson fields [53] which for several types of spin excitations 
has the form [52] 

r](x, y) = i: s T s ip(x)T s 4>(y)G s (x - y) (3.2) 

where 77 has a 4 x 4 matrix form. With a second auxiliary field B(x,y) [52], one gets the 
following generating functional 

Z = N~ 1 J DiljDipDriDBexp[iS(ilj, $,B,ri)+i J dx(ip\ + (3.3) 

S = J dxip{—^.d — M)il) — Tr J dxdyr)(x,y)[B(y,x) — ip(y)^i(x)] (3.4) 
+ \Tr J J dxdyJ2G s (x - y)B(x - y)r s B(y,x)r s 

When the functional integration is carried out over r)(x, y), it gives a 5-function 5[B(y, x) — 
if}(y)ip(x)]. Subsequent integration over B gives eq.(3.1). After this check, the order of 
integration may be reversed so as to integrate out over tp and ■0, and yield the effective 
action 

S = Tr[-iln(-rd-rj)-vB + BB/2]; B(x,y) = Z s G s (x - y)T s B(x,y)T s (3.5) 

Here r], B, (7-<9 + M), are matrices in spinor, internal symmetry, and configuration space 
indices, so that 

r\B = J dzr)(x, z)B(z,y) =< x\i]B\y >; Tr[i]B] —Tr J dx < x\r]B\x > (3.6) 

Varying S w.r.t. B and i] gives 

r}(x,y) = B(x,y); B = i(- 7 .d - M - rj)' 1 = i(-~f .d - M - B)' 1 (3.7) 
Replacing B in (3.7) by the vacuum expectation value < B > = iS F , gives the SDE 

S F = (- 7 .9 — M — tSp)- 1 = Z s G s (x - y)T s S F (x - y)T s (3.8) 
whose detailed form is 

(_ 7 .0 _ M)S F (x -y)-ij dzZ s G s (x - z)T s S F (x - z)T s S F (z - y) = 5\x - y) (3.9) 
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Next, for the quantum corrections to B, write 



B(x,y)=iS F (x-y)+<f>(x,y) (3.10) 
and obtain the homogeneous equation 

inf 

i S F (0S F ) n = tS F ^S F + S F ^; y) = E s G s {x - y)T s <P(x, y)Y s (3.11) 

1 

If the non-linear term in <fi in (3.11) is neglected, the result is the homogeneous BSE 

(f>(x,y)=ij f dzdtS F (x- z)E s G s (z-t)r s (f)(z-t)r s S F (t-y) (3.12) 

which must be solved along with the SDE (3.8) for the propagator. Note that the kernel of 
the BSE is G s , i.e., the same form factor as appears in the input Lagrangian itself. This is 
the basic logic of the interplay of the SDE with the BSE. Next we describe this interplay 
in momentum space for the case T s = ry M A a /2, to bring out the Nambu- Goldstone nature 
of a pseudoscalar state (0 proportional to 75), one in which the Ward identity plays a 
crucial role. 



3.2 Self-Energy vs Vertex Fn in Chiral Limit 

The formal equivalence of the mass-gap equation (SDE) and the BSE for a pseudoscalar 
meson in the chiral limit [24] will now be demonstrated for an arbitrary confining form 
D(k) (not just the perturbative form k~ 2 ). Denoting the mass operator by E(p) and the 
vertex function by T H , the SDE after replacing the color factor Ai.A 2 /4 by its Casimir 
value 4/3, reads as 

E(p) = ^(2vr)- 4 | d 4 kD^(k) lfI S' F (p- D^k) = (5^ - k»h v /tf)D(k) (3.13) 

S' F is the full propagator related to the mass operator E(p) by 

E(p) + i-f.p = Sp 1 (p) = A(p 2 ) [i-f.p + m(p 2 )] (3. 14) 

thus defining the mass function m(p 2 ) in the chiral limit m c = 0. In the same way 
the vertex function IV (g, P) for a qq hadron (H) of 4-momentum made up of quark 
4- momenta p± t 2 = P/2 ±q satisfies the BSE 

V H (q,P) = -U(2ixY A j d A q'D^(q - q'h,S F (q' + P/2)Y H {q\P)S F {q' - P/2) lv (3.15) 

The complete equivalence of (3.13) and (3.15) for the pion case in the chiral limit P M — > is 
easily established. Indeed, with the self-consistent ansatz Y H =7 5 r(g), eq.(3.15) simplifies 
to 

T(q) = ^(2tt)- 4 J d 4 k^S' F (k - q)T(q - k)S' F (q - k) lv (3.16) 

where the replacement q' = q — k has been made. Substitution for S' F from (3.14) in 
(3.16) gives 
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where we have relabelled q — > p. On the other hand substituting for S' F (3.14) in (3.13) 
gives for the mass term of E(p) the result 

A(p*)mtf) = - W< fd'k ^^yL (3-18) 
y ' y ' 3 V 7 J A 2 (q')(m 2 (q' 2 ) + q' 2 ) v 7 

where q' = p — k. A comparison of (3.17) and (3.18) shows their equivalence with the 
identification T(q) = A(q)m(q 2 ) , i.e. the identity of the vertex and mass functions in the 
chiral limit, provided A — 1. this last is a consequence of the Landau gauge for D^ v in 
eq.(3.13), since in this gauge, the function A(p) does not undergo renormalization [54], 
so that it may be set equal to unity. Note that this result is more general than in the 
contact type NJL model, since both quantities are now functions of momentum due to 
the extended nature of the 4-fermion coupling caused by the gluonic propagator D(k). 

3.3 S(p) vs V{q,P) via Ward Identities 

The connection between E(p) and T(q, P) away from the chiral limit (P M = 0) is achieved 
via a systematic use of the Ward identities for vector and axial vector types. The following 
derivation due to [24a] may be instructive for applications. Consider some approximation 
scheme (based on a BSE with a specified kernel) to determine E(p) via eq.(3.13), so as 
to obey the Ward-Takahashi identities. E.g., the quark-gluon vertex function T x satisfies 
the inhomogeneous equation 

T A = 7a - ^(2tt)- 4 J d 4 q' lu S' F (q' + P/2)Y x S' F {q' - Pffl^D^q - q') (3.19) 

Multiplying (3.19) by P\ and using the WT-identity 

P x T x (q, P) = S'^iq + P/2) - S'f\q - P/2) (3.20) 

gives the result 

S> F {P/2 + q)~ S' F {q-P/2) = (3 ' 21) 

rP-jJ 0j4 D ^ ~ JMSfU - P/2) - S F (q> + P/2)] 7 , 

which is entirely consistent with (3.13) when one uses the definition (3.14) for E(p). In a 
similar way, for the axial vector T^, the corresponding BSE obeying chiral symmetry is 

r A5 (g, P) = ?7 A75 - ^(2tt)- 4 J d 4 q'D^(q - q') lv S F (q> + P/2)r A5 (g', P)S F (q' - P/2)^ 

(3.22) 

It is again consistent with eq.(3.15) and the definition (3.14) for E(p) if one uses the axial 
WT identity 

- iP x T X5 (q, P) = S' F \q + P/2) 75 + l5 S' F \q - P/2) (3.23) 

The LHS of (3.23) must now be identified with the pseudoscalar vertex function T 5 (q.P), so 
that the corresponding RHS gives its full structure that is consistent with gauge invariance, 
viz., 

T 5 (q,P) l5 = (3.24) 
i-y.(q + P/2)A(q + P/2) - ij.(q - P/2)A(q - P/2) + B(q + P/2) + B(q - P/2); 
B(p) = A{p)m{p 2 ) 
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This equation checks with (3.18), in the Landau gauge (A = 1), in the chiral limit P M = 0, 
but now provides the corrections for P^ ^ as well. In the Landau gauge (3.24) simplifies 
to 

r 575 = i-y.P + m(q + P/2) + m(q - P/2) (3.25) 

In recent years, the determination of vertex functions via WT identities has become a 
fairly standard practice, although it is not always the most convenient method in practice 
for incorporating gauge-invariance within a given (semi-phenomenoligical) framework. 
For the present report, we shall have occasion to incorporate QED gauge invariance 
in arbitrary momentum- dependent form factors, and the method will be explained in 
Sec. (5), and in more detail in Appendix B, in connection with the P- meson e.m. self- 
energy calculations to be given in Sect. 7 

4 3D-4D SDE-BSE Formalism Under MYTP 

As per the programme outlined in Sect.l, we shall from now on specialize to a more 
practical form of SDE-BSE framework born out of 3D support (defined covariantly) to 
a vector-exchange mediated 4-fermion coupling at the input Lagrangian level with 'cur- 
rent' (almost massless) quarks. The vector exchange simulates the effect of a gluonic 
propagator, encompassing both the perturbative and non-perturbative regimes, and thus 
preserves the chiral character of the input coupling. The derived SDE and BSE, a la Chap 
3, automatically incorporates DP>xS and hence generates the dynamical mass function 
m{p) whose low momentum limit m(0) gives the bulk contribution to the it constituent 
mass m cons , while the current mass m curr for uds quarks (that enter the input Lagrangian) 
gives a small effect. This last is in keeping with Politzer's Additivity principle [55], viz., 
m cons = fn C urr + ^(0), providing a rationale for the quark masses usually employed in 
potential models [12]. 

Now to implement the covariant 3D constraint of MYTP [15] on the BSE kernel 
(which stems from one on the input Lagrangian), we shall consider two methods in parallel 
for a direct comparison: i) Covariant Instantaneity Ansatz (CIA) [16-17]; ii) Covariant 
Null-Plane Ansatz (CNPA) [41]. The latter [41] gives a formal 'covariance structure' to 
an earlier pragmatic formulation with essentially the same content [40], while the former 
[16] is already covariant as it is. We shall now outline a connected account of the 3D 
BSE reduction for both CIA and CNPA types (with scalar followed by fermion quarks), 
to bring out the structural identity of the resulting BSE's for a qq system. This will 
be followed by a reconstruction of the 4D BS vertex functions for both types [16, 41] 
which will serve as the basic framework for 4D quark loop calculations in the subsequent 
chapters. 

4.1 3D-4D BSE Under CIA: Spinless Quarks 

To keep the contents fairly self-contained, we start with a few definitions for unequal 
mass kinematics in the notation of [16,10b]. Let the quark constituents of masses m 12 
and 4-momenta pi^ interact to produce a composite hadron of mass M and 4-momentum 
Pp. The internal 4-momentum g M is related to these by 

Pi 1 2 = mi i2 P±g; P 2 = -M 2 ; 2m 1)2 = 1 ± (m\ - m\)jM 2 (4.1) 
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These Wightman-Garding definitions [56] of the fractional momenta mi i2 ensure that 
q.P = on the mass shells mj+pj = of the constituents, though not off-shell. Now define 
% — % ~ as the relative momentum transverse to the hadron 4-momentum P M 

which automatically gives q.P = 0, for all values of q^. If the BSE kernel K for the 2 
quarks is a function of only these transverse relative momenta, viz. K = K(q, q'), this is 
called the "Cov. Inst.Ansatz (CIA)" [16] which accords with MYTP [15]. For two scalar 
quarks with inverse propagators Ai i2 , this ansatz gives rise to the following BSE for the 
wave fn $(g,P) [16, 10b]: 

i(27r) 4 A 1 A 2 $(g, P)= j d 4 q'K(q, q')$(q', P); A 1>2 = m\ 2 + p\ 2 (4.2) 

The quantities TOi i2 are the 'constituent' masses which are strictly momentum dependent 
since they contain the mass function m(p) [55], but may be regarded as almost constant for 
low energy phenomena m(p) = m(0). Further, under CIA, m(p) = m(p), a momentum- 
dependence which is governed by the DP>xS condition [4] (see below). 

To make a 3D reduction of eq.(4.2), define the 3D wave function <j)(q) in terms of the 
longitudinal momentum Ma as 

(j)(q) = J Mda<S>(q,P); Ma = Mq.P/P 2 (4.3) 
using which, eq.(4.2) may be recast as 

i(27r) 4 A 1 A 2 $(g,P) = j d z q'K(q,q')<P(q')- d 4 q' = d 3 q' Mda' (4.4) 

Next, divide out by AiA 2 in (4.4) and use once again (4.3) to reduce the 4D BSE form 

(4.4) to the 3D form 

(2vr) 3 P(g)0(g) = / d*q'K(q, q')^ ^ = / ^ (4-5) 

Here D(q) is the 3D denominator function associated with the like wave function <p{q). 
The integration over da is carried out by noting pole positions of A 12 in the a-plane, 
where 

Ai, 2 =coi,2 2 - M 2 (m h2 ± a) 2 ; cu h2 2 = m 2 12 + q 2 (4.6) 
The pole positions are given for A 12 = respectively by 

M(a + mi) = ±cu 1 =f ie; M(a - rh 2 ) = ±tu 2 =F ie (4.7) 

where the (±) indices refer to the lower/upper halves of the a- plane. The final result for 
D(q) is expressible symmetrically [16]: 

D(q) = M u DM\ W = ~ + - ^ 

^o(g)=g 2 ~ AK 4 ^ M2) ; A = M 4 - 2M 2 (m\ + m 2 ) + (ml — ml) 2 (4.9) 

The crucial thing for the MYTP is now to observe the equality of the RHS of eqs (4.4) and 

(4.5) , thus leading to an exact interconnection between the 3D and 4D BS wave functions: 

r(4 )s A 1 A^(,,P) = ^« (4.10) 
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Eq.(4.10) determines the hadron- quark vertex function F(q) as a product D(p of the 3D 
denominator and wave functions, satisfying a relativistic 3D Schroedinger-like equation 
(4.5). 

Some comments on the entire BSE structure are now in order. The 'two-tier' character 
of the formalism is seen from the simultaneous appearance of the 3D form (4.5) and the 
4D form (4.4), leading to their interconnection (4.10). The 3D form (4.5) gives the 
basis for making contact with the 3D spectra [13], while the reconstructed 4D wave 
(vertex) function (4.10) in terms of 3D ingredients D and enables the evaluation of 4D 
quark- loop integrals in the standard Feynman fashion [40]. Note that the vertex function 
T = Dip/ {2m) has quite a general structure, and independent of the details of the input 
kernel K. Further, the D- function, eq.(4.8), is universal and well-defined off the mass shell 
of either quark. The 3D wave function is admittedly model-dependent, but together 
with D(q), it controls the 3D spectra via (4.5), so as to offer a direct experimental check 
on its structure. Both functions depend on the single 3D Lorentz-covariant quantity q 2 
whose most important property is its positive definiteness for time-Ike hadron momenta 
(M 2 > 0). 

4.2 CNPA for 3D-4D BSE: Spinless Quarks 

As a preliminary to defining a 3D support to the BS kernel on the null-plane (NP), on 
the lines of CIA [16], a covariant NP orientation [41] may be represented by the 4- vector 
n^, as well as its dual n M , obeying the normalizations n 2 = n 2 = and n.h — 1. In the 
standard NP scheme (in euclidean notation), these quantities are n = (001; -i)/V2 and 
n = (001; i) / \/2, while the two other perpendicular directions are collectively denoted by 
the subscript _L on the concerned momenta. We shall try to maintain the n-dependence 
of various momenta to ensure explicit covariance; and to keep track of the old NP no- 
tation p± = p ± p 3 , our covariant notation is normalized to the latter as p + = n.p\/2; 
p_ = —n.p^/2, while the perpendicular components continue to be denoted by p± in both 
notations. 

In the same nitation as for CIA [16], the 4th component of the relative momentum q = 
rh'iPi —rhip 2 , that should be eliminated for obtaining a 3D equation, is now proportional 
to q n = n.q, as the NP analogue [40] of P.qP/P 2 in CIA [16], where P = p 1 + p 2 is the 
total 4-momentum of the hadron. However the quantity q — q n n is still only q±, since its 
square is q 2 — 2n.qn.q, as befits q]_ (readily checked against the 'special' NP frame). We 
still need a third component p 3 , for which a first guess is zP, where z = n.q/n.P. And for 
calculational convenience we shall need to (temporarily) invoke the 'collinear frame' which 
amounts to P±.q± = 0, a restriction which will be removed later by a simple prescription 
of 'Lorentz completion'. Unfortunately the definition g M = (q±^,zP^) does not quite fit 
the bill for a covariant 3- vector, since a short calculation shows again that q 2 =q\. The 
correct definition is seen as g 3M = zP n n^, where P n =P.n, giving q 2 = q\ + z 2 M 2 , as 
required. We now collect the following definitions/results: 

= q — q n n; q = q± + xP n n; x = q.n/P.n; P 2 = — M 2 ; (4.11) 
Qn,Pn = n.(q, P);q.n — q.n; q.n = 0; P±.q± — 0; 

P.q = P n q.n + P.nq n ; P.q = P n q.n; q 2 = q\ + M 2 x 2 

Now in analogy to CIA, the reduced 3D BSE (wave-fn 0) may be derived from the 4D 
BSE (4.2) for spinless quarks (wave-fn <3>) when its kernel K is decreed to be independent 
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of the component q n , i.e., K = K(q,q'), with q = (q±,P n n), in accordance with the TP 
[15] condition imposed on the null-plane (NP), so that d A q = d 2 qj_dq 3 dq n . Now define a 
3D wave-fn 0(g) = / dq n §(q), as the CNPA counterpart of the CIA definition (4.3), and 
use this result on the RHS of (4.2) to give 

t(2n)^(q) = Ar^- 1 J d 3 q'K(q, q')cf>(q') (4.12) 

which is formally the same as eq.(4.4) for CIA above. Now integrate both sides of eq.(4.12) 
w.r.t. dq n to give a 3D BSE in the variable q: 

(2n) 3 D n (q)<f>(q) = J d 2 q ± 'dq 3 'K(q, q')<f>(q') (4.13) 

which again corresponds to the CIA eq.(4.5), except that the function D n (q) is now defined 
by 

Jdq n A l - 1 A 2 - 1 = 2mD- 1 (q) (4.14) 

and may be obtained by standard NP techniques [40] (Chaps 5-7) as follows. In the q n 
plane, the poles of A 12 lie on opposite sides of the real axis, so that only one pole will 
contribute at a time. Taking the A 2 -pole, which gives 

2 5 „ = = "■S+fo-'^)' (4.15) 

m 2 P.n — q.n 

the residue of Ai works out, after a routine simplification, to just 2P.q = 2P.nq n + 2P n q.n, 
after using the collinearity condition P±.q± = from (4.11). And when the value (4.15) 
of q n is substituted in (4.14), one obtains (with P n P.n = — M 2 /2): 

DM=ZP-n{?- KM \^ ml) Y g 2 = gi + MV; x = q.n/Pn (4.16) 
Now a comparison of (4.12) with (4.13) relates the 4D and 3D wave-fns: 

2m$(q) = A t (g)Ar 1 A 2 - 1 0(g) (4.17) 

as the CIA counterpart of (4.10) which is valid near the bound state pole. The BS vertex 
function now becomes T = D n x (fi/(2ni). This result, though dependent on the NP 
orientation, is nevertheless formally covariant, and closely corresponds to the pedagogical 
result of the old NPA formulation [40], with D n <^> D + . 

A 3D equation similar to the covariant eq.(4.13) above, also obtains in alternative 
NP formulations such as in Kadychevsky-Karmanov [38] (see their eq.(3.48)). Both are 
'covariantly' dependent on the orientation n M of the NP, i.e., have certain n-dependent 
3-scalars, in addition to genuine 4-scalars. However the independent 4-vector n M which 
has a dual interplay with n M in the above CNPA formulation, does not seem to have a 
counterpart in [38]. Secondly this manifestly covariant 4D formulation needs no 3- vector 
like n, or explicit Lorentz transformations, as in such alternative NP formulations [38]. 
As to the 'angular condition', a question first raised by Leutwyler-Stern [35d], no special 
effort has been made to satisfy this requirement, since the very appearance of the 'effective' 
3-vector g M in the 3D BSE in a rotationally invariant manner is an automatic guarantee 
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(in the sense of a 'proof of the pudding') of the satisfaction of this condition [35d] without 
further assumptions. 

A second aspect of the above 3D-4D BSE under CNPA (which allows for off-shell mo- 
menta) is that it has no further need for 'spurions' [38] (to make up for energy- momentum 
balance due to on-shellness of the momenta in such formulations [38]), so that normal 4D 
Feynman techniques suffice, as in the old-fashioned NPA formulation [33,40]. However, to 
rid the physical amplitudes of n M -dependent terms in the external (hadron) momenta, after 
integration over the internal loop momenta, one still needs to employ a simple technique 
of 'Lorentz-completion' (to be illustrated in Sec. 5 for the pion form factor calculation) as 
an alternative to other NP prescriptions [37,38] to remove n-dependent terms. 

A more succinct comparison with other null-plane approaches concerns the inverse 
process of reconstruction of the 4D hadron-quark vertex, eq.(4.17)), which has no coun- 
terpart in them [37,38], as these are basically 3D oriented. Thus in [38], the nearest 
analogue is to express the 3D NP wave function in terms of the 4D BS wave function (see 
eq.(3.58) of [38]), but not vice versa. This problem of 'loss of Hilbert space information' 
inherent in such a process of reconstruction, has been discussed recently in the context of 
the qqq problem [47]; (see Section 10 for details). 

4.3 Fermion Quarks with QCD-Motivated BSE 

We are now in a position to give a corresponding description of the 3D-4D BSE for fermion 
quarks, for both CIA and CNPA cases taken together, just as for spinless quarks above. 
The 4D BSE for fermion quarks under a gluonic (vector-type) interaction kernel with 3D 
support has the form [10 a-b]: 

z(2tt)^(P, q) = S F1 ( Pl )S F2 (p 2 ) J d*q'K(q, q')V(P, q'); K = F 12 i^i^V(q, q') 

(4.18) 

where F 12 is the color factor Ai.A 2 /4 and the V- function expresses the scalar structure 
of the gluon propagator in the perturbative (o.g.e.) plus non-perturbative regimes. The 
'hat' notation on the momenta covers both CIA and CNPA cases simultaneously, where 
the longitudinal component q^ is defined for the CNPA case as = zPnii^, with P n =P.h. 
The full structure of V (used in actual calculations [28,40]) is collected as under, using 
the simplified notations k for q — q', and V(k) for the V fn: 

V(k) = Ana s /k 2 + -^J dv[r 2 {\ + 4A m 1 m 2 M > 2 r 2 )~ 1/2 - C /uj 2 }e lLr ; (4.19) 

u 2 qq = 4M > m 1 m 2Wo 2 a s (M > 2 ); a s (Q 2 ) = ^-gL_ ln(M > /A)- 1 ; (4.20) 

m 1>2 = [1 ± (m\ -m 2 2 )/M 2 ]/2; M> = Max(M, mi + m 2 ); C = 0.27; A = 0.0283 

(4.21) 

And the values of the basic constants (all in MeV) are [28,40] 

co = 158; m ud = 265; m s = 415; m c = 1530; m b = 4900. (4.22) 

The BSE form (4.18) is however not the most convenient one for wider applications 
in practice, since the Dirac matrices entail several coupled integral equations. Indeed, as 
noted long ago [10 a-b], a considerable simplification is effected by expressing them in 



24 



'Gordon-reduced' form, (permissible on the quark mass shells, or better on the surface 
P.q = 0), a step which may be regarded as a fresh starting point of our dynamics, in 
the sense of an 'analytic continuation' of the 7- matrices to 'off-shell' regions (i.e., away 
from the surface P.q = 0). Admittedly this constitutes a conscious departure from the 
original BSE structure (4.18), but such technical modifications are not unknown in the 
BS literature [lOc-d] in the interest of greater manoeuvreability, without giving up the 
essentials, in view of the "effective" nature of the BS kernel (see Chap 1 sec. 6). 
The 'Gordon-reduced' BSE form of (4.18) is given by [lOa-b] 

A 1 A 2 $(P, q) = - l (2nr 4 F 12 J d 4 q'V^V^V(q, q')^P, q'); (4.23) 

where the connection between the \t r - and $- functions is 

V(P,q) = (mi -i7 (1) .p 1 )(m 2 + i 7 (2) .p 2 )<l>(P,g); p lj2 = m lj2 P ± q (4.24) 

= ±2mi, 27 i 1 ' 2) ; V® = p ¥ + + Ur®{p* ~ Pi) (4.25) 

Now to implement the Transversality Condition [15] for the entire kernel of eq.(4.23), 
all time-like components a, a' in the product V^.V^ must first be replaced by their 
on-shell values. Substituting from (4.25) and simplifying gives 

(Pi+Pi)-(P2+P 2 ) = 4mim 2 P 2 ~(q + q') 2 - 2(m 1 -m 2 )P.(q + q') + ll spin -Terms"; (4.26) 

"SpmTerms" = -i{2m x P + q + q')^}k v + *(2m 2 P - q - q')^}k v + a$a$ (4.27) 

This is identical to eq.(7.1.9) of Ref. [40], via the correspondence q 44> q x , q y , qs(= q+/P+), 
so that both CIA [16] and CNPA [41] have formally the same structures as the 'old- 
fashioned' NPA [40], and hence give identical predictions on the 2-body spectra [28]. 

The 3D reduction of eq.(4.23) now goes through exactly as in the spin-0 case, eqs.(4.2- 
8), so that without further ado, the full structure of the 3D BSE can be literally taken 
over from Ref.[40]-Chap 7 (derived under old-fashioned NPA). In particular, for harmonic 
confinement, obtained by dropping the Aq term in the 'potential' U(r) of (4.23) (a very 
good approximation for light (ud) quarks), the 3D BSE works out as 

P(g)0(g)=^-P(g)0(g); (4.28) 

D n (q)m = (4.29) 

for the CIA and CNPA cases respectively, where D n is given by (4.16) and D(q) by (4.8). 
The other quantities retain the same meaning for both. Thus 

D(q) = 4m!m 2 M 2 (V 2 + C /cj 2 ) + 4g 2 V 2 + 8g.V + 18 - 8J.S + (4C /cJ 2 )<f (4.30) 

For the spectroscopic predictions on qq hadrons, vis-a-vis data, the reader is referred 
to [28]. For algebraic completeness however, the (gaussian) parameter /3 of the 3D wave 
function 0(g) = exp(—q 2 /2/3 2 ), which is the solution of (4.28-29) for a ground state hadron 
[40,28] is: 

. 8m 2 mnM 2 u;ntt(M 2 ) 9 „ n A t . ^ ^ , , ;/l9 /.„,n 

P 4 = h * \ > 1,, ~ ; < o > 2 = 1 + 24A m 1 m 2 M > f f3 2 4.31 

[1 — 86 mim 2 Q; s (M^)J < a > 
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Note that (3 is a 4D invariant quantity, independent of n/i, etc. (For an L-excited hadron 
wave function, see [40]). The full 4D BS wave function ty(P,q) in a 4x4 matrix form [40] 
is then reconstructed from (4.23-24) as in the scalar case, eq.(4.10), viz., [40,16,44] 

V(P,q) = S F ( Pl )T(q) lD S F (-p 2 )- T(q) = N H [1; P n /M]D{q)<t>{q)/^ (4-32) 

where 7^ is a Dirac matrix which equals 75 for a P-meson, 27^ for a V-meson, 27^75 
for an A-meson, etc. The factors in square brackets stand for CIA and CNPA values 
respectively. N H represents the hadron normalization given by (see Appendix A): 

Nx 2 = 2(2tt) 3 J d^[M w ;M}(f) 2 (q)[(l+5m 2 /M 2 )(f-X/4M 2 )+2m 1 m 2 (M 2 -5m 2 )} (4.33) 

where M u is given in (4.8), and 5m = mi — m 2 , and again the factors in square brackets 
represent CIA/ CNPA values. 

4.4 Dynamical Mass As DB\S Solution of SDE 

We end this Section with the definition of the 'dynamical' mass function of the quark 
in the chiral limit (M n = 0) of the pion-quark vertex function T(q), in the 3D-4D BSE 
framework [11]. The logic of this follows from the BSE-SDE formalism outlined in Sec. 3, 
eqs.(3. 12-18), for the connection between eq.(3.18) for m(p) and eq.(3.17) for T(q) in the 
limit of zero mass of the pseudoscalar. So, setting M = in (4.28-29), the scalar part of 
the (unnormalized) vertex function may be identified with the mass function m(p), in the 
limit P^ = 0, where p^ is the 4- momentum of either quark; (note the appearance of the 
'hatted' momentum). The result is [11] 

m(p) = [u;(p); V2p.n] q / 0(]3) (4.34) 

m\ 

under CIA and CNPA respectively. The normalization is such that in the low momentum 
limit, the constituent ud mass m q is restored under CIA [11], while the corresponding 
'mass' under CNPA is p + [35c]. 

A more important aspect of the 'dynamical' mass function is its appearance as the 
non-trivial of the SDE under DP>xS [9,24]. We now give a derivation of the 3D-4D 
counterpart [11] of this basic result [24]. To that end we start with the non-perturbative 
part of the gluon propagator D /JU/ (k) = D{k)[5 liv — k tl k v /k 2 ] for the (harmonic) interaction 
of ud quarks (forming the 'pion'structure), where the scalar factor D{k) has the form [11] 

D(k) = - A {2Kfujl2m q a s {Am 2 q )[Vl + C /u 2 ]5\k) (4.35) 

This form is immediately derivable from the structure of the 'potential' function Vk) of 
eq. (4. 19-20), with the A -term set equal to zero, and taking M> = 2m q for the 'pion' case. 
Note that D(k) has a directional dependence n M = P^jP 2 on the pion 4-momentum P M , 
so that k 2 > over all 4D space; it also possesses a well-defined limit for P M — > 0. This 
structure may now be substituted in the 'gap-equation' (3.17/18) foa a self-consistent 
solution in the low momentum limit. This exercise has been carried out in [11], wherein 
the SDE (3.18) in the Landau gauge A{p 2 ) = 1 reduces to the form 

m(p 2 ) = ^J d*kdkom q a s x [^ 2 V, 2 + C ]P(k) {p ,^ 2 \ p>2) ( 4 - 36 ) 
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where p' = p — k is 4D, and (k, ko) are (3D, ID) respectively. The integration is essentially 
over the time-like k , with the 'pole' position at p ' = m{p' Q ) = m NJL , leading finally to 

[11] 

3m q a s 2^21 6tt 

m NJL = 3w - C m NJL ; a s = — r 4.37 

rriNjL 29ln(10m q ) 

after substituting the value A = 200MeV for the QCD constant. The further identification 
of ra 9 with rriNjL in this equation, yields an independent self-consistent estimate m^jL ~ 
SOOMeV", which may be compared to the input value 265MeV, eq.(4.22), employed for 
spectra [28]. This analysis so far ignores the Politzer relation [55] m ud = m c + m NJL , 
for the constituent mass m q away from the chiral limit. The derivation of the pion and 
cr-meson masses away from the chiral limit, may be found in [11]. 



5 CNPA Applications: Gauge-Inv Pion F.F. 

The first example of our applications of the 3D-4D BSE structure developed in Chap 4 is 
to 4D triangle loop integrals. This example has been chosen to illustrate the difficulties of 
CIA (as noted in Chap 1.4) in tackling their ill-defined nature as a result of acquiring time- 
like momentum components in the exponential/gaussian factors associated with the vertex 
functions (4.10) due to a 'Lorentz-mismatch' among the rest-frames of the concerned 
hadronic composites, for triangle loops and above, such as the pion form factor, while 2- 
quark loops [32] just escape this pathology. This problem was not explicitly encountered 
in the old-fashioned NPA treatment [33] of the pion form factor, except for lack of explicit 
covariance. The CIA approach [16] to MYTP [15] enjoys covariance, but its application 
to triangle loop integrals causes other problems such as complexity of the corresponding 
amplitudes [57], apparently without good reason. On the other hand, the apparent success 
of the old-fashioned NPA [33] in circumventing this (complexity) problem [57], gives the 
hope that with its 'covariant' formulation [41], CNPA, the powers of this method should 
stand a better chance of testing via the form factor problem. 

To recall a short background, the pion form factor has through the ages been a good 
laboratory for subjecting theoretical models and ideas on strong interactions to observa- 
tional test. Among the crucial parameters are the squared radius < r 2 xp = 0.43±.014/m 2 
[58a], and the scaled form factor at high k 2 , viz., k 2 F(k 2 ) 0.5 ± O.lGeV 2 [58b] that 
represent important check points for theoretical candidates such as QCD-sum rules [2b], 
Finite Energy sum rules [59], perturbative QCD [60], covariant null-plane approaches 
[37,38], Euclidean SDE [61], etc. The issue interface of perturbative and non-perturbative 
QCD regimes has been studied in terms of the relative importance longitudinal vs trans- 
verse components [63], but this is the subject of a full-fledged dynamical theory (such as 
[9-12]), and not of some intuitive ansatze [62]. 

To that end, we outline a calculation of the P-meson form factor for unequal mass 
kinematics with full gauge invariance, including correction terms arising from QED gauge 
invariance, and illustrate the techniques of 'Lorentz-completion' to obtain an explicitly 
Lorentz invariant quantity. As a check on the consistency of the formalism, the expected 
k~ 2 behaviour of the pion form factor at high k 2 is realized. Some calculational details 
on the triangle-loop integral for the P-meson form factor are given in Appendix 5. A. 
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5.1 P-Meson Form Factor F(k 2 ) for Unequal Masses 

Using the two diagrams (figs, la and lb) of ref.[33c], and in the same notation, the Feyn- 
man amplitude for the h — > h! + 7 transition contributed by fig. la ( quark 2 as spectator 
is given by [33c] 

2P,F(P) = A{2nYN n {P)N n {P')em 1 J 2] . (5-1) 

4T <i) = y r [ 75 ( mi _ 27.^)17^(7711 - ^7-^1)75(^2 + 17-^2)]; Aj = m- + p 2 ; (5.2) 

Pi,2 = mi, 2 P ± g; pi 2 = m 1)2 P' ± (/ p 2 =p' 2 ; P - P' = Pl - p[ = k; 2P = P + P'. 

(5.3) 

After evaluating the traces and simplifying via (5.2-3), T M becomes 

T« = (p2„ - P M )[W — M 2 — A 2 ] - A; 2 p 2M /2 + (A x - (5.4) 

The last term in (5.4) is non-gauge invariant, but it does not survive the integration in 
(5.1), since the coefficient of k^, viz., Ai — A[ is antisymmetric in pi and p[, while the 
rest of the integrand in (5.1) is symmetric in these two variables. Next, to bring out the 
proportionality of the integral (4.1) to P M , it is necessary to resolve p 2 into the mutually 
perpendicular components p 2 j_, (p2-k/k 2 )k and (p 2 .P/P 2 )P, of which the first two will 
again not survive the integration, the first due to the angular integration, and the second 
due to the antisymmetry of k — p\ —p[ in Pl and p[, just as in the last term of (5.4). The 
third term is explicitly proportional to P M , and is of course gauge invariant since P.k = 0. 
(This fact had been anticipated while writing the LHS of (5.4)). Now with the help of 
the results 

p 2 .P = -m 2 M 2 -A 1 /4-A' 1 /4; 2m 2 = 1 - [m\ - m 2 )/M 2 ; P 2 = — M 2 — /c 2 /4, (5.5) 

it is a simple matter to integrate (5.1), on the lines of Sec. 4, noting that terms proportional 
to Ai A 2 and A' : A 2 will give zero, while the non-vanishing terms will get contributions only 
from the residues of the A 2 -pole, eq.(4.15). Before collecting the various pieces, note that 
the 3D gaussian wave functions 0, 0', as well as the 3D denominator functions D n , D' n , do 
not depend on the time-like components p 2n , so that no further pole contributions accrue 
from these sources. (It is this problem of time-like components of the internal 4-momenta 
inside the gaussian ^-functions under the CIA approach [16], that had plagued a earlier 
CIA study of triangle diagrams [57]). To proceed further, it is now convenient to define 
the quantity q.n = p 2 .n — m 2 P.n to simplify the 0- and D n - functions. To that end define 
the symbols: 

(<L <?') — Q =t ?7i2&/2; Z2 = q.n/P.n; k = k.n/P.n; (9k, rjk) — 1 ± k/4 (5.6) 
and note the following results of pole integration w.r.t. p 2 „ [40]: 

J dp 2n i-[l/Ai; 1/Ai; 1/(A 1 A' 1 )] = [1/D n , 1/D' n , 2p 2 .n/(D n D' n )} (5.7) 

The details of further calculation of the form factor are given in Appendix A. An essential 
result is the normalizer N n (P) of the hadron, obtained by setting k^ = 0, and demanding 
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that F(0) = 1. The reduced normalizer Nh = N n (P)P.n/M, which is Lorentz- invariant, 
is given via eq.(A.9) by: 

N H 2 = 2M(2tt) 3 J d 3 ^ 2 ^ 2 [(1 + 5m 2 /M 2 )(q 2 - A/4M 2 ) + 2m 1 m 2 (M 2 - 5m 2 )} (5.8) 

where the internal momentum q = (q±,Mz2) is formally a 3-vector, in conformity with 
the 'angular condition' [35d]. The corresponding expression for the form factor is (see 
Appendix A): 

F{k 2 ) = 2MN 2 H (2ixfexp[-(Mm 2 k/P) 2 /A6 k ](nP 2 f /2 ^m l G{k) + [1 2] (5.9) 
where G(k) is defined by eqs.(A.12-13) of Appendix A. 
5.2 'Lorentz Completion' for F{k 2 ) 

The expression (5.9) for F(k 2 ) still depends on the null-plane orientation n M via the 
dimensionless quantity k = k.n/P.n which while having simple Lorentz transformation 
properties, is nevertheless not Lorentz invariant by itself. To make it explicitly Lorentz 
invariant, we shall employ a simple method of 'Lorentz completion' which is merely an 
extension of the 'collinearity trick' empolyed at the quark level, viz., P±.q± = 0; see 
eq.(4.11). Note that this collinearity ansatz has already become reduntant at the level of 
the Normalizer Nh, eq.(5.8), which owes its Lorentz invariance to the integrating out of the 
null-plane dependent quantity z 2 in (5.8). This is of course because N H depends only on 
one 4-momentum (that of a single hadrori), so that the collinearity assumption is exactly 
valid. However the form factor F(k 2 ) depends on two independent 4-momenta P, P', for 
which the collinearity assumption is non-trivial, since the existence of the perpendicular 
components cannot be wished away! Actually the quark- level assumption P±.q± = 
has, so to say, got transferred, via the g-integration in eq.(5.9), to the hadron level, as 
evidenced from the /c-dependence of F(k 2 ); therefore an obvious logical inference is to 
suppose this /c-dependence to be the result of the collinearity ansatz P±.P' ± = at the 
hadron level. Now, under the collinearity condition, one has 

P.P f = P ± .P' ± + P.nP'.n + P'.nP.n = P.nP^ + P'.nP n ; P.n = P n . (5.10) 

Therefore 'Lorentz completion' (the opposite of the collinearity ansatz) merely amounts 
to reversing the direction of the above equation by supplying the (zero term) Pj_.P' ± to 
a 3-scalar product to render it a 4-scalar! Indeed the process is quite unique for 3-point 
functions such as the form factor under study, although for more involved cases (e.g., 
4-point functions), further assumptions may be needed. 

In the present case, the prescription of Lorentz completion is relatively simple, being 
already contained in eq.(5.10). Thus since P, P'= P ± k/2, a simple application of (5.10) 
gives 

Ak 2 

k.nk n = +k 2 ; P.nP n = -M 2 -k 2 /4; k 2 = ^ - ^ = A6 k - 4 = 4 - 4 Vk (5.11) 

This simple prescription for k automatically ensures the 4D (Lorentz) invariance of F(k 2 ) 
at the hadron level. (It may be instructive to compare this to the Cov. LF prescription 
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of 'recognizing' the n- dependent terms (unphysical) of F(k 2 ) and then dropping them). 
For more involved amplitudes (e.g., 4-point functions) too, this prescription works fairly 
unambiguously, if their diagrams can be analyzed in terms of more elementary 3-point 
vertices (which is often possible). We hasten to add however that strictly speaking, a 
'Lorentz completion' goes beyond the original premises of restricting the (pairwise qq) in- 
teraction to the covariant null-plane (in accordance with MYTP [15]), but such 'analytic 
continuations' are not unwarranted, since in Cov. LF theories too [38], implementing the 
angular condition [35d] involves the introduction of 'derivative' terms, implying a tacit 
enlargement of the Hilbert space beyond the null-plane (see Chap 2 of [38]). 

5.3 QED Gauge Corrections to F(k 2 ) 

While the 'kinematic' gauge invariance of F(k 2 ) has already been ensured in Sec. 5.1 above, 
there are additional contributions to the triangle loops - figs. la and lb of [33c] - obtained 
by inserting the photon lines at each of the two vertex blobs instead of on the quark lines 
themselves. These terms arise from the demands of QED gauge invariance, as pointed out 
by Kisslinger and Li (KL) [63] in the context of two-point functions, and are simulated by 
inserting exponential phase integrals with the e.m. currents. However, this method (which 
works ideally for point interactions) is not amenable to extended (momentum-dependent) 
vertex functions, and an alternative strategy is needed, which is described below. 

The way to an effective QED gauge invariance lies in the simple-minded substitution 
Pi — eiA{xi) for each 4-momentum pi (in a mixed p, x representation) occurring in the 
structure of the vertex function. This amounts to replacing each occurring in T(q) = 
D(q)(j)(q), by q^ — e q A^, where e q = rh 2 ei — and keeping only first order terms in 

after due expansion. Now the first order correction to q 2 is —e q q.A — e q A.q, which 
simplifies on substitution from eq.(5.11) to 

- 2e q q.A = -2e g A M [g M - q.nn^ + P.nq.nn^/ P.n] (5.12) 

The net result is a first order correction to T(q) of amount e q j(q).A where 

j(q), = -4M>g M 0(g)(l - (q 2 - \/AM 2 )/2(3 2 ) (5.13) 

The contribution to the P-meson form factor from this hadron-quark-photon vertex (4- 
point) now gives the QED gauge correction to the triangle loops, figs.(la,lb) of [33c], to 
the main term F(k 2 ), eq.(5.1), of an amount which, after a simple trace evaluation (and 
anticipating the vanishing of all A-terms remaining in the trace, as a result of contour 
integration over q n ) simplifies to (0 = 0(g), etc) 

F^k 2 ) = A{2^N 2 H e q m l M 2 > J d 4 q(M 2 - 6m 2 W[^j^ + + [1 =► 2]; 

(5.14) 

In writing down this term, the proportionality of the current to 2P M has been incorporated 
on both sides, on identical lines to that of (5.1), using results from (5.2-5.7) as well as 
from Appendix A. Note that e q is antisymmetric in '1' and '2', signifying a change of 
sign when the second term [1 =>- 2] is added to the first. The term q.P/P 2 simplifies 
to 2q.n(l — k/2)/P.n, after extracting the proportionality to P^. Next, after the pole 
integrations over q n , q' n in accordance with (5.7), it is useful to club together the results 
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of photon insertions on both blobs for either index ('1' or '2'); this step generates two 
independent combinations for the '1' terms (and similarly for '2' terms): 

A n = q.n(l - fc/2); B n = q.n(l - k/2)(q 2 - A/4M 2 )/2/3 2 (5.15) 

Collecting all these contributions the result of q n - integration is 

F^k 2 ) = ^fNle.m.Ml [ d 3 q(M 2 - 6m 2 )<M>'[ K + K =3 ~ K ] + [1 =► 2] (5.16) 

The rest of the calculation is routine and follows closely the steps of Appendix A for the 
(main) F{k 2 ) term, including the translation z 2 — > z 2 + rh 2 k 2 /29k, and is omitted for 
brevity. The final result for Fi{k 2 ) is 



.(Ml - 5m 2 )%fc 4 (M > m 2 ) 
8G(0)9 7 k /2 lP 



Fl {k 2 ) = -e^rh^m + k 2 ) r> z:::^ ^ i + [i 2] (5.i 7) 



where we have dropped some terms which vanish on including the [1 =>• 2] terms, noting 
the (1, 2) antisymmetry of e q . 

5.4 Large and Small k 2 Limits of F(k 2 ) 

We close this section with the large and small k 2 limits of the form factors F(k 2 ) and 
F^k 2 ). For large k 2 , eq.(5.11) gives k= 2, 9 k = 2, and r] k = 4M 2 /k 2 , so that 

AM 2 

F{k 2 ) = 2MN 2 {2nfm 1 ^{7i(3 2 /2) 3 / 2 G{mf)exp[-{Mm 2 /(3) 2 /2] + [1 2] (5.18) 

where, from eqs.(A.ll-12), 

G(inf) = (1 + <5m 2 /M 2 )(/3 2 - A/4M 2 + M 2 m 2 ) + (M 2 - <5m 2 )m 2 - 2m 2 M 2 (5.19) 
Similarly from eq.(5.14), the large k 2 limit of F\(k 2 ) is 

F.ik 2 ) = 2V2M 2 k- 2 e q m 1 m 2 (M 2 > - W)[ M >^~f a) ] (5.20) 

where we have taken account of the (1,2) antisymmetry of e q in simplifying the effect 
of the [1 =^ 2] term on the RHS. As a check, both F(k 2 ) and F 1 (k 2 ) are seen to satisfy 
the 'scaling' requirement of a k~ 2 variation for large k 2 . This result can be traced to the 
input dynamics of the (non-perturbative) gluonic interaction, eq.(4.19), on the structure 
of the vertex function, eq.(4.32). Perturbative QCD of course gives a k~ 2 behaviour [61]. 
The covariant NP(LF) approach [38] also gives a similar behaviour, but extracted in a 
somewhat different way from the present 'Lorentz completion' treatment. Note that for 
the pion case the QED gauge correction term Fi(k 2 ) gives zero contribution in the large 
k 2 limit. 

For small k 2 , on the other hand, we have from eq.(5.11) 

k 2 = k 2 /M 2 ; [9 k , Tfc] = 1 ± k 2 /4M 2 (5.21) 
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In this limit, the form factor, after substituting for Nh from (5.8), and summing over the 
'!' and '2' terms, works out as 



F(k 2 ) = (1 - 3k 2 /8M 2 )[l - mim 2 ( 



k 2 k 2 5m 2 



3k 2 p 2 (l + 5m 2 /M 2 ) 
8M 2 G(0) 



] (5.22) 



4/3 2 M 2 G(0) 



where G(0) is formally given by eq.(A.lO), except for the replacement of q 2 by 3/3 2 /2. As a 
check, F(k 2 ) is symmetrical in '1' '2', as well as satisfies the consistency condition F(0) = 
1. Similarly the small k 2 value of Fi(k 2 ), after taking account of the (1, 2) antisymmetry of 
e q , is of minimum order k 4 , so that it contributes neither to the normalization (i*i(0) = 0), 
nor to the P-meson radius. 

For completeness we record some numerical results for large and small k 2 limits. For 
the pion case, in the large k 2 limit, eqs.(5.12-13) yield after a little simplification the 
simple result 



where m q = 265MeV stands for mi = m 2 ; and M> stands for the bigger of mi + m 2 
and M. Substituting for j3 2 =§M§3GeV 2 [32] and G(0)=0.166Ge\/ 2 , yields the result 
C=0.35GeV 2 , vs the expt value of 0.50 ±0.10 [58b]. For comparison, we also list the per- 
turbative QCD value [60] of Snagf^ 2 = 0.296GCV 2 , with f 7T =133MeV , and the argument 
Q 2 of a s taken as M 2 . 

For low k 2 , eqs.(5.14-15) yield values of the pion and kaon radii, in accordance with 
the relation < R 2 > = — Vfc 2 F(/c 2 ) in the k 2 =0 limit. Substitution of numerical values 
from (4.21-22) yields 

R K = 0.629fm(vs : 0.53 - expt[58a]); R w = 0.661fm(vs : 0.656 - expt[58a]) (5.24) 

We end this Section with the remark that a simple-minded, conventional NP approach 
[33, 40] to BS dynamics had already produced most of the results of this form factor 
calculation, but had been critized [34] on grounds of 'non-covariance'. The CNPA with 
an explicit formulation of the Transversality Principle (TP) [10-8] on a covariant null 
plane (NP), hopefully, keeps both the advantages, since the 4D loop integrals are now 
not only perfectly well-defined, but a major part of the dependence has got eliminated 
in the process of q integration, while the remaining NP orientation dependence has been 
transferred to the external (hadron) 4-momenta. In this regard the present approach is 
already in the company of a wider NP(LF) community [37-38] which has also to contend 
with some n M dependence. The solution offered here to overcome this problem is a simple 
prescription of 'Lorentz completion' wherein a 'collinear frame' ansatz P±.q± = is lifted 
on the external hadron momenta P, P' etc, after doing the internal q integration, so as to 
yield an explicitly Lorentz- invariant result. The prescription, though different from other 
LF approaches [38], is nevertheless self-consistent, at least for 3-point hadron vertices, 
(and amenable to extension to higher-point vertices provided the latter can be expressed 
as a combination of simpler 3-point vertices). (It may be added parenthetically that the 
old-fashioned NP treatment [33c] had yielded a slightly better curve for the pion form 
factor, but this was due to the use of the "half-off- shell" form of the NP wave function 
[40], which however did not come out naturally from the present 'covariant' treatment). 



F(k 2 ) = C/k 2 ; C = 2v / 2-^r(/5 2 + m 2 q )e- M >^ 2 




(5.23) 
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6 Three-Hadron Couplings Via Triangle-Loops 



For a large class of hadronic processes like H — > H' + H" and — > if' + 7, the quark 
triangle loop [31] represents the lowest order "tree" diagram for their evaluation. Criss- 
cross gluonic exchanges inside the triangle (see Fig.l of [31] ) are not important for this 
kind of description in which the hadron-quark vertices, as well as the quark propagators 
are both non-perturbative, and thus take up a lion's share of non-perturbative effects. This 
is somewhat similar to the "dynamical perturbation theory" of Pagels-Stokar [64] in which 
such criss-cross diagrams are neglected. 

In this Section we shall give an outline of the calculational techniques for such diagrams 
for the most general case of unequal mass kinematics mi 7^ m 2 7^ m 3 , but with spinless 
quarks only, since the 'spin' d.o.f. does not introduce any new singularities over the 
spin-0 case. In this we shall closely follow the method of ref. [31], which is a 3-hadron 
generalization of Sec. 5 for the e.m. form factor of a pseudoscalar meson. However, as 
already noted therein, the CIA form [16] of 3D-4D BSE is fraught with problems of ill- 
defined integrals (and hence complexity of amplitudes) due to the presence of time-like 
momentum components [25] in the (gaussian) wave functions of the participating hadrons. 
So we shall work only with CNPA [41] structures, as derived in Sec. 4. 2. 



6.1 Kinematical Preliminaries 

According to Fig.l of ref [31], and in the same notation, the 3 hadrons with all incoming 
4-momenta Pj, with masses M$, interact via the quark triangle loop wherein P k dissociates 
into the quark pair with 4-momenta {—Pi,Pj) and masses {m^rrij) respectively, so that 
P k = - p . + Pj: and P k + Pi + Pj = 0. Thus [31]: 

- P k = P t + Pj = P ij; Pi = Ml P k = -Pij = p 3 - Vl (6.1) 

where (i,j,k) are cyclic permutations of 1,2,3). Similarly the relative 4-momenta q k 
between quarks i,j corresponding to the break-up P k = Pj — p%, and Q k between hadrons 
i,j for the break-up P k = —Pi — Pj are: 

qk = fcjPj + frjiPi, Qk = rhijPj - rhjiPi; -P k = Pi + Pj (6.2) 

The fractional momenta fiij at the quark level, and rhij at the hadron level, are given by 
the Wightman-Gaerding [56] definitions 

m\-m 2 . Mf — M 2 , 

2£tt = l+ M 2 J ; 2m tJ = l+ * J (6.3) 

The relative signs are determined by the phase convention of Fig.l of [31]. 

Now to define the 'hatted' relative 3-momenta q k and Q k , we must follow the CNPA 
procedure [41] instead of CIA [31]. Further, since the content of CNPA is for all practical 
purposes identical with that of the old-fashioned NPA [40], considerable simplification is 
achieved by adopting the latter notation [40], which is what is already done in [31], albeit 
with CIA content. Indeed, with the collinear ansatz (Sec. 4) the NPA values of ^ are 
simply qi±,q i3 , where q i3 = Miq i+ /P i+ , etc [40]; and Q i3 = MiQ i+ /P i+ . However, since 
the (?j's are not all independent, it is necessary to take a basis momentum (say p 2 ) in 
terms of which to express others. Now in a fixed pi basis, we have 

Qk=Pi + p-ijPk] Qj = Pj ~ kkPj] q% = Pi ~ frjkPj + hj p k (6.4) 
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For later purposes we shall consider a p 2 basis, for which 

qi=P2 + ^2?,Pl\ 93=^2-/^21^3 (6.5) 

We also record some useful results for the kinematics of external particles, if they are 
on-shell (Qi-Pi = 0), under the collinearity condition [31]: 

P i± = ±^Q, ± ; M? = P t+ P t+ - Q\ = -Q.^- = ml ^ Ml) (6.6) 

which lead to the further symmetry relations 

Q 1 M 1 = Q 2 M 2 = Q 3 M 3 = y/\{M*,Ml,M$)/4 (6.7) 

Further we can define a 3 x 3 matrix structure n^, with n^- = P i+ /Pj + , which satisfy the 
relations 

riijUji = 1; ni 2 n 23 n 31 = 1; Ej7i a = (6.8) 

showing a definite phase relation among these quantities, which are more explicitly ex- 
pressed by the matrix structure 

K\ = l;-m 13 ±Q 2 /M 2 ;-m 12 ^Q 3 /M 3 (6.9) 
-rh 23 T Qi/Mi; 1; -m 21 ± Q 3 /M 3 (6.10) 
-m 32 ±Q 1 /M 1 ; -m 31 =p Q 2 /M 2 ; 1 (6.11) 

with a two-fold sign ambiguity expressed by the statement that only the upper, or only 
the lower signs, must be taken. It is easily verified that eqs (6.8) (hence phase relations) 
are satisfied by the matrix (6.9). 

6.2 Structure of HHH Form Factor 

The full structure of the 3-hadron amplitude due to Fig.l of [31] is 

AW-^r/***^ (6.12) 

exhibiting cyclic symmetry, where the normalized vertex function Tj in CNPA [41] is given 
from Sec. 4-5 as 

m) = N^y^D.mm a = 2m^ 2 - mh ^ ml \ (6.13) 

where we have defined the 'reduced' denominator function Di as Di + Mi/Pi + and written 
the (invariant) normalizer NiH as iVj. The color factor and the effect of reversing the loop 
direction are given by 2/v^3, while (2ir) 8 is the overall BS normalizer [40]. Aj = mf + pf 
=ojf 1 _ —p i+ pi-. Spin and flavour d.o.f. will give rise to a standard 'trace' factor [31] [TR] 
which is skipped here for simplicity. 

To evaluate (6.10), we first write the cyclically invariant measure: 

<fp i = d p J^d{xl)Mfdy i] x i =p i+ /P i+ ; y i =p i -/P i . (6.14) 
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The cyclic invariance of this quantity ensures that it is enough to take any index, say 2, and 
first do the pole integration w.r.t. the y 2 variable which has a pole at y 2 = $, 2 =u 2± / (M 2 x 2 ). 
The process can be repeated, by turn, over all the indices and the results added. Note 
that the ^-functions do not include the time-like yi variables under CNPA [41], so that the 
residues from the poles arise from only the propagators. The crucial thing to note is that 
the denominator functions Di and D 3 sitting at the opposite ends of the p 2 Tine in Fig.l 
of [31] will cancel out the residues from the complementary (inverse) propagators A 3 and 
Ai respectively. Indeed by substituting the pole value y 2 = £2, in Ai ;3 , the corresponding 
residues in an obvious notation work out as: 

Ai ;2 = 6n 32 M 2 2 + x 2 n 2Z Ml - 2/} 21 M 3 2 ; A 3;2 = -£ 2 n l2 Mj - x 2 n 2X M\ - 2jl 23 M* (6.15) 

It is then found, with a short calculation using (6.5), that 

^M = 2M 3 x 2 n 23 ; = 2M 1 x 2 n 21 (6.16) 

Al;2 A 3;2 

which shows the precise cancellation mechanism between the Dj-functions and the residues 
of the propagators A« at the A2 pole. This mechanism thus eliminates [16] the (overlap- 
ping) Landau-Cutkowsky poles that would otherwise have caused free propagation of 
quarks in the loops. The same procedure is then repeated cyclically for the other two 
terms arising from the A 3j i poles. Collecting the factors, the result of all the 3 contribu- 
tions is compactly expressible as (c.f. [31]): 

A(3H) = 8^S 123 J J M 2 n 2 3n 21 n 2 dx 2 d&xl[TR} 2 D 2 (q 2 )n 123 M i N i cj) i (6.17) 

where the limits of integration for both variables are — inf < (£ 2 , x 2 ) < +inf, since these 
are governed, not by the on-shell dynamics of standard LF methods [37-38], but by off- 
shell 3D-4D BSE. The difference from [31] (under CIA [16]) arises from using CNPA [41] 
here. 



6.3 Discussion on Applicability 

Eq.(6.14) is the central result of this Section. Its general nature stems from the use of 
unequal mass kinematics at both the quark and hadron levels, which greatly enhances 
its applicability to a wide class of problems which involve 3-hadron couplings, either as 
complete process by themselves (such as in decay processes) or as part of bigger diagrams 
in which 3-hadron couplings serve as basic building blocks. What makes the formula 
particularly useful for general applications is its explicit Lorentz invariance which has 
been achieved through the simple method of 'Lorentz Completion' on the lines of Sec. 5 
for the e.m. form factor of P-mesons (pion). 

How much of this derivation is model independent, except for the use of the MYTP 
[15] ? The answer lies in the structure T# = D x <fi for the hadron-quark vertex function, 
which is a direct consequence of the 3D support ansatz which in turn receives support 
from several angles [17-19], although this specific form [16] does not seem to have been 
used elsewhere. Its factorable structure in which the denominator function D is quite 
universal and depends only on kinematics, has helped reduce the 4D loop-integral to a 3D 
form, and in so doing, has succeeded in eliminating the Landau-Cutkowsky (overlapping) 
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singularities in a very simple and transparent manner, thus preventing the free propagation 
of quarks in the intermediate (loop) stages. (The only model dependent entity in the 3D 
wave function 0, but it has been related to the (observable) spectra [13]). 

In the spirit of this generality, this article was not intended for specific allpications per 
se, but some possibilities are readily listed. The simplest class is that of strong decay of a 
resonance (Hi) into two lighter hadrons (Hi,H 2 ) under kinematically allowed conditions 
(whose signature is carried by the external variables Qi,Mi inside the integral (6.14)). 
The amplitude A(3H) can also be adapted, via See. 5, to include e.m. or semi-leptonic 
processes, expressed by Hi — > H 2 + 7, (where 7 is real or virtual), where the signature 
of virtuality is carried by M| = t. Non-leptonic weak decays (see Fig.2 of [31]) are also 
amenable to this treatment. As an example one may cite the experimental discrepancy [65] 
of the vector form factors in the semi-leptonic process D K*eu with theoretical models 
that prefer to represent the intermediate states through effective meson propagators [66] . 
On the other hand, a crucial role of the appropriate quark-triangle, with a considerable 
effect of the unequal masses of the participating quarks, seemed to be strongly indicated 
in resolving the discrepancy [67]. Other applications include the so-called Sullivan process 
[68], details of which the interested reader may find in [69]. 

7 Two-Quark Loops: SU(2)-Breaking Problems 

To illustrate other applications of the 3D-4D BSE formalism, we now turn to the (simpler) 
problem of two-quark loops which are useful for estimating SU(2)-breaking effects in 
phenomena like i) mass-splittings in P-mesons [18a], and ii) p — uo mixing in meson- 
exchange forces [18b]. Simpler 2-quark loops, such as those involved in the weak and e.m. 
decay constants of hadrons, are already available in previous studies of this formalism 
[40], and will not be the subject of this semi-review. Further, its scope does not include 
detailed analysis of these 2-loop phenomena [18], but only their essential physics, and a 
quick derivation of their core structures, leaving the reader to ref [18] for numerical results 
plus more references. 

7.1 Strong SU(2)-Breaking in P-meson Multiplets 

To recapitulate the essential physics of hadronic mass splittings within SU(2)multiplets 
(I = 0.5, 1.0), these were for long thought to be of e.m. origin, until the advent of QCD 
[22] when the possibility of strong breaking of SU(2) due to the intrinsic u — d mass dif- 
ference started being taken seriously. (This was despite the prior existence of the GMOR- 
mechanism [70] which had sought to relate the pseudoscalar masses to the current quark 
masses and the vacuum condensates ). In this respect the trend was set largely by Wein- 
berg's analysis [71], characterized by the 'Weinberg ratios' m u /md = 0.55, and m s /m u d = 
20.1, confirmed by a recent analysis [72]. A conservative estimate of the u-d mass differ- 
ence is believed to be d — u = 3 — AMeV [71-72]. On the other hand the absolute values 
of the current masses are not as well known, but the SU(2) mass splittings [13] among 
the known pseudoascalar multiplets (ir, K, D, B) is a useful mathematical laboratory to 
determine the d — u mass difference from the corresponding 'constituent mass' difference, 
via Politzer additivity [55]. The issue is basically a dynamical one (in view of the sensi- 
tive nature of this labiratory), necessitating a high degree of parametric control on the 
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strong vertex functions involved in the concerned Feynmam diagrams (Figs.l(a,b,c) of 
ref.[18a]). The problem clearly goes beyond mere additivity in the quark masses, as the 
observed pattern of mass splittings [13] seems to suggest a basically decreasing trend from 
the lightest (pion) to the heaviest (beauty) flavour, tapering off almost to zero for the 
B — B + mass difference. 

For the meson self-energy, there are 3 basic contributions, a la Figs.l(a,b,c) of [18a]: 
i) In Fig. la, a 2-point 5m U( i vertex inserted at each propagator by turn represents the 
principal source of strong SU(2) breaking; ii) Fig. lb simulates the e.m. breaking effect 
by joining the two quark lines internally by a photon propagator; iii) Finally Fig.lc 
simulates the effect of the difference of the quark condensates < uu > and < dd > on the 
strong SU(2) breaking of hadron masses. Figs.l(a,c) are one loop diagrams, while Fig. lb 
represents a two-loop process, which is moreover sensitive to QED gauge constraints [63], 
as in the e.m. form factor case (see Sec. 5). 

Using the dynamical framework collected in Sec. 4. (2-3), it is fairly straightforward to 
write down the integrals accruing from these diagrams. Both CIA [16] and CNPA [41] 
are valid mechanisms for evaluating these diagrams, but in view of a prior exposure [18a] 
of CIA for this problem, it may be instructive to adopt the CNPA alternative here. Thus 
Fig. la gives in terms of the results of the previous sections, 

n a (M 2 ) = t(27r)- 1 N 2 H J d^qD^q^^Tr^Sp^P + q) l5 S F (-m 2 P + q)} (7.1) 

where we have used the representation of the normalized vertex function given in eq.(6.11), 
and D is the reduced denominator fn in CNPA. After evaluating the traces this expression 
simplifies to 

II a (M 2 ) = -2iWN% J A^ 2 (g)0 2 (g) Al + A2 A ^ 2 " W (7.2) 

where 8m = mi-m2, and all kinematical quantities are as defined in Sec. 4 The integration 
over the time-like component g_ is carried out very simply, using the result (5.7). Note 
that the vertex function D x does not involve this variable, and also the D-fn exactly 
cancels out the residues arising from the propagators, as shown generally in Sec. 6. The 
resultant 3D integration over d 3 q is expressible simply as 

II a (M 2 ) = 2N 2 H J d^mm^ + H| + M 2 - 5m 2 } (7.3) 

where Xi= p i+ /P + = rhi ± x for i= 1,2 respectively. The third component q 3 of CNPA 
[40-41] is simply Mx, so that q = q±,q3- The normalizer N H is given by eq.(5.8). The 
parallel CIA result is [18a] 

IT a (M 2 ) = 2N 2 H J d 3 q[D 2 (q)(^- + -L) + D((q)(M 2 - 5m 2 )] (7.4) 

A comparison between the CNPA and CIA forms of Il a is now in order. In CIA, 
eq.(7.4), the D 2 term is well defined and is amenable to simple quardature. On the other 
hand, the CNPA form, eq.(7.3), encounters singularities at x± >2 = 0, on integration w.r.t. 
x, taking account of the relations Xi i2 =m 12 — x, and q 2 = q\ + M 2 x 2 . The final results 
are quite similar for both cases. 
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The formulae (7.3-4) for n a and (5.8) for Nh, show explicit dependence on the masses 
777,1,2, and facilitate the evaluation of mass splittings within the SU(2) isospin multiplets as 
follows: For the K, D, B mesons, take m 2 as the mass of the -urf-quark with m 1 > m 2 , and 
while differentiating w.r.t. m 2 , consider the increment 8 C . A little reflection then shows (by 
virtue of the Politzer [55] additivity relation) that this quantity may be directly identified 
with the difference m dc — m uc between the current d- and u -masses provided the hadron 
mass with the w-quark gets subtracted from that with the d-quark (e.g.,i^o — K- , etc). 
Of course the normal rules of differentiation apply, viz., Sf(m 2 ) = f'(m 2 )8 c , where the 
argument of f'(m 2 ) must use the average 'constituent mass' of -ud-quark, viz., 265MeV, 
eq(4.22). 

For the pion case, some extra care is necessary since both the constituents are now 
■u/rf-quarks so that both m\ and m 2 must be subjected to differentiation in turn . On 
the other hand these two contributions come with just equal but opposite signs, so that 
they cancel out exactly, giving a net vanishing contribution, as seen more directly from 
the fact that the AI = 1 field uu — dd in the Lagrangian cannot contribute to 7r + — 7r 
anyway. For the details of numerical results on <5II(M 2 ), see [18a]. 



7.2 E.M. Contribution to Self-Energy 

The e.m. contribution to the hadron self-energy is given by a 2-loop diagram (Fig lb of 
[18a]) in a slightly simplified notation as follows: 

r r DrhD'6' 

U b (M 2 ) = N 2 H e 2 ei e 2 J J dW j^j^Tr[ lb S F { Pl ) i^S F (p[) l5 S F (-p' 2 ) i^S F (-p 2 )} 

(7.5) 

where k=q — q'=Pi — p'\=p' 2 — P2 is the exchanged quantum ; e\ and e 2 are the charges of 
the quarks involved ( in units of e ) and q, q' are the internal 4-momenta of the LHS and 
RHS hadrons respectively. This integral involves simultaneous (pole) integrations over the 
time-like components of q and q' which do not figure in the respective vertex functions 
and therefore can be carried out exactly. However the rest of the 3D integrations (two 
sets) do not quite factor out, so they need some strategy before they can be carried out 
without much tears. To that end a simple device that suggests itself naturally is based on 
the following observation : By the very topology of the diagram it is fairly clear that the 
time-like components of both the 4-vectors q and q' are quantitatively similar, so that 
their effects largely "cancel out" in the factor k~ 2 in eq.(7.5). As a result the quantity k 
=q — q' effectively reduces to the space-like quantitty (q — q') 2 which can be manipulated 
to desired numerical accuracy in the 3D integrations over q and q' . We list both CIA [18a] 
and CNPA (new) results in the form of 3D integrals in q, q' jointly as follows. 

Tl b (M 2 ) = 4N 2 H e 2 eM^r 3 j j d^hat q '[..]^^- (7.6) 
the quantity [...] is first listed for CIA as follows [18a]: 

[..]cia = (M 2 -5m 2 ) 2 -2m 1 m 2 (M 2 -m 2 12 )-5m 2 (q-q') 2 + (7.7) 
' +7^-)(7^-l + 7^-l)D(q)D(a') 



1m\ 2m 2 1vri\ 2m 2 
+ (M 2 - ml - 2m 2 2 + m 1 m 2 )D(g)/cJ 2 + (M 2 - m 2 - 2m\ + m^^D^/ui, 
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The dual quantity [...]cnpa may be simply read from the above merely by the replacements 
^1,2— >2(M#i, Mx 2 ) respectively, where 2+2=m li2 ± x - 

To convert the mass shifts from quadratic to linear, it is of course necessary to divide 
both U a b in the avobe equations by 2M. In the e.m. case, no further 'differentiation' 
w.r.t. m 2 is necessary, since (7.6) is already of second order in e. As regards the factor 
eie 2 , its differential 5(eie 2 ) is easily found as +(1/2), —(1/3), +(2/3) and —(1/3) for the 
differences ir + — 7r , K° — K + ,D + — D°, and B° — B~ in this order .It turns out [18a] that 
this alternating sign pattern is of great help in reinforcing and reducing the net differences 
within the K,D,B multiplets (after taking account of the strong breaking effects, Figs 
la and lc), towards a good pattern of accord [18a] with the data [13]. (For results, see 



We now consider QED gauge corrections [63] to the e.m. value, eq.(7.6), arising from 
Fig. lb of [32a], on the lines of corresponding corrections to the e.m. form factor derived 
in Sec. 5. This correction is sketched in Appendix B for P-mesons, using diagrams listed in 
ref.[63], and in their notation for the contributing figures. The resulting QED correction 
for the kaon e.m. mass difference turns out to be nearly a 60 percent increase over the 
CIA result 1.032Mel / [18a] arising from the main term (7.6). We omit the corresponding 
CNPA treatment for brevity. 

7.3 Effect of Quark Condensates 

Another source of mass splittings arises from the difference between the uj rf-quark conden- 
sates, in accordance with Fig.(lc) of [18a]. Indeed some recent calculations via QCD-sum 
rules have used this as the principal mechanism [73] for the mass splittings, with much less 
contribution from Fig. (la). Indeed the value of S < qq > in itselh has been the subject of 
separate investigations in chiral perturbation theory [74] as well as in QCD-sum rules [74]. 
On the other hand, the BSE-SDE formalism [11,23] provides a 'direct' ab initio estimate 
[11] of this condensate (as well as others [30]). 

To recapitulate the logic of the condensate calculation by the 'direct' method [11], 
in terms of the quark's non-perturbative mass function, m(p), note that the latter is the 
chiral (M n = 0) limit of the pion-quark vertex function T(q), given by eq.(4.35), and must 
be used in the expression of the full propagator, Sf(p), Sec. 7, which appears in the formal 
definition of the condensate as follows [11]: 



after doing the pole-integration over the time-like component of p M . Here N c = 3, and 
Nf = 1 (since each separate flavour (u/d) is counted). Now to evaluate the 3D inte- 
gral (7.8), substitute the CIA structure (4.35) for m(p), with 4>(p)=exp(—p 2 /2f3 2 ). This 
integral formula has an analytic structure in terms of the constituent mass m q of the 
«/(i-quark, so that it is now a matter of simple differentiation to give an explicit form of 
its increment w.r.t. Sm q which equals 5 C . The final formula is [18a]: 



Using the inputs from (4.21-22) gives f3 2 = 0.0603, and the final results under CIA are 



[18a]). 




(7.8) 




(7.9) 



< qq >= -(266Me\/) 3 ; 5 < qq > = +0.06645c 



(7.10) 
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These values are fully rooted in spectroscopy but are otherwise free from adjustable param- 
eters, except for the quantity 5 C . They have a fair overlap with QCD-SR determinations 
[76]. 

For completeness we now give the condensate results under CNPA, substituting the 
CNPA mass function (4.35) in (7.8). This gives 

rW2 f , 3 „, P-n^ + ^lW) 



< qq >= ttttj" / d pdpn 2 , 2 — r (7.ii) 

(2tt) 4 J m\+p\- 2p.np n 

The integration over p n is trivial and yields 

K qq >= / d3p[1 + = -( 242MeV ) 3 ( 7 - 12 ) 
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Substituting the gaussian form (as above) for (ft and integrating, yields an analytic struc- 
ture useful for calculating 5 < qq >: 

< qq >= -3V2(P 2 /2nf 2 [l + 3[3 2 /m 2 q ] = -{2A2MeVf (7.13) 

a value which seems to be even closer to the estimate — (240) 3 of QCD-SR [2] than the 
CIA result -(266) 3 of [18a]. 

As to the contribution of 5 < qq > to the strong SU(2) mass splittings, a la Fig.lc of 
[18a], we skip the detailed derivation in favour of [18a], since it turns out to be rather 
small within this BSE framework. This is in sharp contrast to the QCD-SR findings 
[73] wherein the condensate contribution seems to dominate. This is not too surprising 
since within a BSE-cum-SDE framework, most of the non-perturbative effects are already 
contained in the hadron-quark vertex function, with a correspondingly smaller role for the 
condensates. On the other hand in QCD-SR [2] these represent major non-perturbative 
effects when seen from the high energy perturbative QCD end. 

A few comments on the main results of this exercise are in order. The e.m. contribu- 
tions alternate in sign in the mass splittings between the charged and neutral components 
in the sequence n, K, D, B. The condensate contribution to strong SU(2)-breaking being 
small, the sensitivity to the d — u mass difference comes almost entirely from Fig. 1(a) 
of [18a]. Next, the feature of unequal mass kinematics has played a big role in the for- 
malism, being mainly responsible for a sustematic decrease in mass splittings as one goes 
up on the mass scale. This aspect has come about mainly from the properties of the 
/^-functions (mostly model independent) The numerical values show a good overall pat- 
tern of agreement with data [13], (within less than half MeV), for the parameter S c in the 
range (3.5 — 4.0)MeV for all the 3 cases K, D, B. This value of S c appears well within the 
phenomenological limits of acceptability [72]. However, as the results of Appendix B on 
QED gauge corrections indicate, inclusion of these tends to decrease the effective value 
of 5 C . Finally the calculational technique seems to conform to the spirit of 'Dynamical 
Perturbation Theory' of Pagels-Stokar [64] (neglect of 'criss-cross' diagrams) which must 
be carefully distinguished from a naive interpration of perturbative QCD. 

7.4 Off-Shell p-u Mixing 

Before concluding this Section, we shall briefly draw attention to a similar SU(2)-breaking 
phenomenon which has proved to be of considerable interest for the understanding of 
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certain anomalies in nuclear forces [77]: off-shell p — uj mixing. Although nuclear topics 
are not of direct concern for this article, the basic logic of charge-symmetry-breaking 
(CSB) to explain the Nolen-Schiffer anomaly [77] via p — uj mixing [78], stimulated by 
new experiments [79] on polarized n — p scattering, comes directly under the theme of this 
Section. Indeed, the sensitivity of p — uj mixing to the d — u mass difference 8 C , especially 
off-shell [78a, 18b], is as strong as that of P-meson masses [18a]. 

To recall the basic logic, the small difference between the proton vs neutron analyzing 
powers at an angle 6q corresponding to the vanishing of the average analyzing power [79] , 
is proportional to the CSB potential Vqsb whose contribution from p — uj mixing may be 
schematically expressed as [78a] 

V£sb =< NN\H int \NNu >G < uj\H CS b\p° > G < p°NN\H mt \NN > +(p° & uj) 

(7-14) 

Here Go is the appropriate V-meson propagator, and < uj\Hcsb\p° > gets its dominant 
theoretical contribution from the d — u mass difference <5 C , with H C sb = P-^8 C 2 , and a 
partial contribution from the e.m. chain p =>- 7 =>- uj via vector dominance and/or 2- 
quark loops. Alternatively, the matrix element can be estimated from the experimental 
e + e~ =>- 7r + 7r~ amplitude at the u;-pole, which gives the on-shell value 6(M 2 ) of the p — uj 
mixing amplitude [80]. On the other hand, it is its off- shell value 6(q 2 ) which is relevant 
to the CSB potential, eq.(7.14), for the V^-meson exchange in a space-like region where 
its effect on Vcsb has been claimed to be greatly suppressed [78a]. This question in turn 
requires a theoretical model for the necessary extrapolation which can be defined in terms 
of a dimensionless parameter A as [18e]: 

6(q 2 ) = 6(M 2 )[l-(l + q 2 /M 2 )\] (7.15) 

A calculation of this parameter A is the central issue of any investigation of the CSB 
effect, wherein its value has been variously estimated to be within the (0 — 1) range [78]. 
In particular, the function 6(q 2 ) is also amenable to the 3D-4D formalism [16], using the 
self- energy techniques [18b] outlined in this Section. Its on-shell value 9(M 2 ) [18b] agrees 
with the data [80], while the off-shell parameter A comes close to unity, signifying a change 
of sign for 9{q 2 ) in the transition region between the space-like and time-like momenta, 
in agreement with a 'nucleonic' self-energy calculation [78d]. 

8 QCD Parameters from Hadron Spectroscopy 

In this Section we outline a simple method of calculation [30] of QCD Condensates in 
terms of the (spectroscopy-oriented) parameters of the 3D-4D BSE framework. These 
parameters of QCD simulate non-perturbative effects as coefficients in Wilson's operator 
product expansions {OPE) [81,55]. The method of QCD sum rules represented the first 
practical attempt [2a] to relate these quark-gluon quantities to hadronic amplitudes by 
employing a duality principle [20] between the quark-gluon and meson-baryon pictures. 
Basically the idea is to find a Q 2 region 1 GeV 2 ) where one may incorporate non- 
perturbative physics, generated via OPE [81], into the perturbative QCD treatment of 
physical processes involving hadrons. The QCD-SR ansatz [2] for the evaluation of a cer- 
tain correlation function n(p) , is to replace the free quark (or gluon) propagator by one 
more suitable for the nontrivial vacuum, and on the other hand to express, via dispersion 
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relations, the same correlation function in terms of the variables of the meson-baryon 
picture. The two results are then equated to yield sum rules connecting the variables of 
the two physical descriptions. 

8.1 Field-Theoretic Definition of Condensates 

While QCD-SR per se [2,82] is not the subject of this review, its basic building blocks- 
the condensates, are the main concern of this Section. These may be defined in terms of 
quark- and gluon- fields [82,30]: 

< qO iQ >= £ < 0| : #(0)0^(0) : |0 >= -J -0^TrS$(p)O t , (8.1) 

where is an operator representing the nature of condensate, the index A represents the 
effect of a background field, and Sp~(p) is the quark propagator with the perturbative part 
suitably subtracted. At this stage, we must distinguish between the gluonic background 
field and other external ones (electromagnetic, axial, etc.): The latter can be taken per- 
turbatively, but the former, with its characteristic problem of color gauge invariance, must 
be addressed more fully, a subject on which there exists a vast literature [83]. However 
it is possible to incorporate in practice a major fraction of this effect through the simple 
device of changing the variable of integration in eq.(8.1) from p^ to = p^ — ^g s X a G^, 
where G a is the gluon field. This would in general not be possible if one were to evalu- 
ate complicated integrals involving more propagators and vertex functions, but since the 
integral in (8.1) "sees" only one such quantity, the trick should work, especially since we 
are mainly interested in a constant background G^ v — field, i.e. G®(x) = — |rr M G^„. This 
is basically a non-abelian adaptation of the famous Schwinger method [43] to the present 
situation but the details of the available methods [83] are not necessary for justifying this 
step. With this understanding, we shall not use any additional subscript or superscript 
in (8.1) to specify the gluonic background, but rather take the integration variable p^ to 
represent = p M - \g s X a G^. 

The principal quark condensate < qq >o corresponds to Oi — 1 and A = 0. The 
corresponding gluon condensate is defined as 

< g 2 s G 2 >= Tr(V M V, - 5^ 2 )g 2 D^(0), (8.2) 

where V M is the gauge covariant derivative and D^ v (x) is the non-perturbative part [11] 
of the gluon propagator. These quantities which are free parameters in QCD-SR, provide 
access to the non-perturbative domain of QCD, but except for the two principal conden- 
sates < qq > and < g 2 G 2 > , which are amenable to cross checks against many data, 
the determination of the higher order ones often leave ambiguities. A partial list is [82] 

_1 _1 
< qil^q >A, <q(T ilu q> F , <q-Xa-Gq> , < q-XG^q > F . (8.3) 

In the method of QCD-SR [2,82], there is no intrinsic mechanism to evaluate them from 
first principles but only an extrinsic 'matching' between the two sides of the duality rela- 
tion with the help of suitable parameters. And for condensates of still higher dimensions, 
additional assumptions, such as factorization, are needed. The BSE-SDE framework [9- 
12] on the other hand, has a more microscopic structure which gives it simultaneous access 
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to both high and low energy phenomena under one umbrella. Thus the condensates (8.3) 
as well as others, are calculable within such a framework with as much ease as the (low 
energy) spectroscopy is accessible to it (See Sec. 1.4 for discussions thereof). The same 
facility also holds for its 3D-4D adaptation which provides a two-tier structure, with the 
3D sector specifically attuned to spectroscopy, while the 4D structure is good for loop 
integrals, thus naturally giving rise to a spectroscopic linkage between the high and low 
energy descriptions of hadrons via QCD. To that end eqs (3.12-18) of the BSE-SDE inter- 
play [82], adapted to its 3D-4D form [11], are collected in Sec. 4. 4: i) an explicit expression 
(4.33) for the mass function m(p) derived from the condition that it is the pion-quark 
vertex function in the chiral limit of M n =0; ii) the non-perturbative gluon propagator 
D(k), eq.(4.34); iii) its more general form V(k), eqs. (4. 19-20); iv) and the formula (4.30) 
for the inverse range (5 of the 3d wave function 0. These are the main ingredients needed 
for the condensate calculations in this Section. 



8.2 The Gluon Condensate in 3D-4D Formalism 

We start by rewriting the gluon propagator in a more general form than (4.34) by making 
use of the more complete ^-function, eqs. (4. 19-20), as under: 



where a, b are the color indices in the adjoint representation. The logic of the connection 
between the V^-fn, and the \^-fn Eqs. (4. 19-20) comes about from its relation with the 
fermionic kernel K: 

K{q,q')^ llx V llv {q-q') lv . (8.5) 

, where the scalar part V{k) in the infrared region may be identified with the confining 
part of the ^-function as [11,28a]: 

V{k) = |(27r)^ 2 2m g a c (4m;)[ 2 + ^]5 3 (k), (8.6) 

4 J 1 - A to^V| w o 

using the full q — q potential V(k) which fits the spectroscopy for all flavors (light and 
heavy) [28], but specialized to the equal mass (m q ) case. The constants Cq, ujq, Aq are 
given by eqs. (4.21-22), while the QCD coupling constant a s is given by [28]: 

47T 1 

a » = ii _2/v ■ As = 200 MeV - ( 8 - 7 ) 

The coordinate representation T>(R) of the gluon propagator (8.6) is 

V(R) = -ul ■ 2m q a s (Q 2 )[^ - . ]. (8.8) 

4 ^l + A m 2 q R 2 

Note the, in the R — > oo limit, T>(R) is linear in R as well as flavor independent (the 
m q — factor cancels out), except for the a s (Q 2 ) effect. Thus the structure (8.8), despite 
its empiricity, respects the standard QCD constraint, but only in the strict confining 
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region. On the other hand, the smallness of Aq (= 0.0283) ensures that for light flavors 
its structure is dominated by the harmonic form, which amounts to setting A = 0. This is 
an excellent approximation for the pion- vertex function in the chiral limit (M n = 0), and 
hence for the quark mass function given by (4.33) [11], and normalized to m(0) = m q and 
identified as the constituent mass for the w<i-quarks only (ignoring their small 'current' 
values). The 3D wave function 0(g) is a gaussian with inverse range f3 given by (4.30), 
which for equal masses simplifies to 

#«) = «*(-W); ^ 2m ^W . (8 . 9 ) 
2 1 - 2a s (4m^)G 

For the inputs (4.21-22), (3 2 works out as (O.OQO)GeV 2 . 

We shall use the mass function (4.33)/(8.9) in the next subsection for the quark con- 
densates. Here we indicate briefly a derivation of the gluon condensate, eq.(8.2), by 
inserting the gluon propagator (8.8) in its definition. The Co— term may be dropped as 
it will not survive the subsequent differentiations in eq.(8.2). For the main term, the 
following integral representation is employed: 

R 2 2m, i R 2 

dR ° i I a ^ 2D2 > ( 8 - 10 ) 



l + A Q m 2 q R 2 2 ™ J l + Avm 2 q R 2 ' 

where R 2 = R 2 — R^ (Lorentz-invariant). The 4D expression T>(R) may now be inferred 
from its definition in terms of the 3D quantity T>(R): 

a s (Q 2 ) 2my VA R 2 
V{R) ~ 7T 1 + A m 2 q R 2 ■ (8 - U) 

This is as far as one goes by adopting the 3D form (8.8) for T>(R). However, it is sufficiently 
suggestive of the extrapolation needed to make it fully covariant, viz. R 2 — > R 2 in the 
numerator of eq.(8.10), which we adopt in what follows. [On the other hand, if this 
replacement it not adopted, then the resulting gluon condensate will be reduced by a factor 
3/4]. The full propagator in the Landau gauge is (8.4), where k^ is read as k^ = —idj}. 
To evaluate the gluon condensate we first note the result: 

^ < GlA >= i-^d^(R) + 2d*d*V%{R)\\ R=(i , (8.12) 

and obtain by straightforward differentiation 

< g 2 G 2 >= ^A A7ra s (Am 2 q )(6m q cu ) 2 /n 2 . (8.13) 

The remaining question concerns what value of the quark mass m q , i.e. what flavor, 
should be employed for evaluating the gluon condensate. The structure (8.8) does exhibit 
the desired features of linear confinement and flavor independence, but the extrapolation 
of these features in the opposite limit (R — > 0), as demanded by eq.(8.14), brings in 
an "effective flavor dependence" of the final formula (8.15). The heavier the flavor, the 
more important is the corresponding mass {m q ), vis-a-vis the A Q — term in the qq potential 
(8.8). Since, on the other hand, the full potential (8.8) fits all the flavor sectors rather 
well [28], a simple "weighting" procedure was chosen in [30], involving only the 3 flavor 
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sectors with a nontrivial flavor mass, viz. ss, cc, and bb with equal weights (in the sense 
of a geometric mean), taking account of the m q - dependence m^a s (4m^) of eq.(8.15). This 
gives the result 

< m 2 q a s (Am 2 q ) >= 13.91{m 2 a s (4m 2 )}, (8.14) 

in units of its value in the (-uci)-region, and its substitution in (8.15) yields the final 
estimate < g 2 G 2 > = 0.502 GeV\ versus the value of 0.47 GeV 4 adopted in the QCD 
sum rule literature [82]. 



8.3 < qOiq > Condensates 

We now substitute the mass function m(p 2 ), eqs.(4.33)/(8.9) into the general formula 
(8.1), to derive the various condensates for different choices of Oi. As already noted 
in Sec. 8.1 (in light of color gauge invariance), the quantity p M in eqs.(4.33)/(8.9), and 
everywhere else in the following, must be read as 11^ [43], with appropriate non-abelian 
corrections. The formula (8.1) now reads as 

< qO iq >o= Tr-^- [d*Tl m ^-^ U Qi (8 . 15 ) 

(2tt) 4 7 m 2 (n) + ( 7 .n) 2 1 ; 

in the absence of external fields. Note that the subtracted part with m(lV)—0 in this 
equation gives no effect on tracing in the absence of external fields. We first express the 
denominator in an alternative form: 

m 2 (n) + ( 7 .n) 2 = Cj 2 - s 9 - n 2 = a - s 3 ; n 2 = n 2 -n 2 (8.16) 

^ 2 =m 2 (fl) + fl 2 ; E 9 = i^lA a G>^, (8.17) 

where 11/ is the longitudinal component of ILj, d i Tl=d 3 TldTli, and the integration must 
first be carried out over II/. Because of the presence of the E s -term in (8.18), however, a 
further "rationalization" of eq.(8.15) is necessary according to the identity 

1 1 A + E « (8.18) 



A-E 9 i. 2 -n?-E 9 



At this stage it is probably adequate to replace E 2 in the denominator of (8.18) by its 
spin-color- averaged value < E 2 >: 

E 2 ^< E 2 >= 1 < g 2 GG >= /i 4 (= 8.48mJ) (8.19) 

after the necessary substitutions have been made from (8.14) and (4.21-22). Thus < E 2 > 
contains a strong signature of the gluon condensate whose large value introduces some bad 
analyticity properties in the denominator of the integrand in (8.15) or (8.18), for purposes 
of II;— integration, since the Cj 2 — term is numerically much smaller than /i 2 . It has been 
emphasized [30] that this feature has nothing to do with the 3D-4D BSE treatment, since 
we have not yet passed the barrier of the orthodox 4D quark propagator in the background 
of the gluon field. It is rather a very general manifestation of the strong spin-color effect 
of the quark-quark interaction via the color magnetic field. The problem is not so serious 
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in QED [43] where the smallness of the coupling constant leaves the counterpart of the 
/i 2 term well below the positivity limit ( i.e., Qj 2 — /i 2 > 0), but the large value of /i 2 in 
the present (QCD) case tends to invalidate the standard analyticity structure of (8.15) 
for purposes of further integration with respect to d 3 U. This problem could not be solved 
in [30], but it seems to deserve more serious attention from a wider community. (Taken 
literally, it would imply the introduction of a. phase in the condensates !) In the meantime, 
a conservative view [30] was taken that the maximum allowed value of < S 2 > (consistent 
with the positivity of the denominator after Ih-integration) should not exceed a) 4 for all 
values of II 2 , i.e. 

< S 2 >= a 2 < m\. (8.20) 
Thus eq.(8.18) should be understood as 

1 _ A + S *. A -,.2 A2 



A-S 3 A 2 - a- 2 ' 



to 2 - Uf. (8.21) 



For the numerator of eq.(8.21) which still carries the spin-dependent quantity S s , 
eq.(8.17), there is no restriction of magnitude for one S s — factor only, since it contributes 
to condensates only after contracting with another X— factor in eq.(8.15). [However, other 
factors which come from the rationalization of the denominator with higher powers of S ff 
must be subject to the same restriction]. With this precaution, eq.(8.15) serves to define 
two condensates simultaneously, viz, these with, Oi — 1 and Oi = g s (A a /2)G a a^, where 
the latter is expressible in the notation of Ref. [82] as 

< 0|g2S 3 g|0 >= m 2 < qq > . (8.22) 

To evaluate the integral over dlli, we have 

1 f m A; a 

where 



/^A^ = [/(<7);J(<7)1 ' (8 ' 23) 



£ Vuj z — a \ / uj z + a 
After collecting the necessary trace factors the final result for the two condensates is 
expressible as a simple quadrature (q = u or d): 

3 r°° - - 2 < S 2 > 

<gg> [l;m 2 ] = — / IT 2 dU m(U) [1(a); g —J(<r)\. (8.25) 

n z Jo a 

On insertion of the structure (4.33)/(8.9) for the mass function, and putting the "maxi- 
mum allowed value" of a, viz, m 2 eq.(8.20), the results under CIA are [32a] 

< qq > = (266 MeVf; ml = 0.130 GeV 2 ; (8.26) 

these results may be compared with the QCD-SR (input) values [82] of (240MeV) 3 and 
0.8GeV 2 respectively. The corresponding CNPA [41] result, as worked out in Sect. 7.3 
with m{p) obtained from Sect. 4. 4 is (242MeV) 3 for the first item. 

We next calculate three induced condensates Xi K i an d J, due to a constant external 
e.m. field F^, which are defined as [82]: 

< <7<V<7 >f= eeqxFvu, < qq > ; (8.27) 
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g s < q(\ a /2)G%q > F = ee qK F^ < qq > ; (8.28) 

9s < q(\ a /2)e lluafS G a afS q > F = ee^F^ < qq > . (8.29) 

In these equations the relative phases of the induced condensates are defined with respect 
to the main condensate < qq > , in accordance with (8.15), and this feature must be kept 
systematic track of. Like the two condensates (8.25), the quantities x an d k are in a sense 
dual to each other, and are best described together. The e.m. field is introduced through 
the substitution 

[rh + i 7 .n] -> [rh + i 7 „(n^ - eA„)} (8.30) 

in the propagator, eq.(8.1), and keeping only the first order term in A^. Thus we have to 
calculate 

J TrS F (U)(ie r A)S F (U)[a, u ; \\ a G%,\. (8.31) 
This is facilitated, for a constant e.m. field, by the representation 

1 d 

A„ = ~-x v F^ Xfl = i—. (8.32) 

The substitution in eq.(8.25) and subsequent trace evaluation is routine but lengthy. 
However certain precautions are necessary in the matter of extraction of two groups of 
terms, proportional to a^ u and G^ v respectively, before the trace evaluation, which will 
survive contraction with the external e.m. field F^ v . Thus, 

n M n, => l - 9s l -\ a G%- 7/l7 „ ia^. (8.33) 

In terms like ia^xTlxTlv, additional survivors come from the symmetrized product {IIa, 11^} 
for which we make the standard isotropy ansatz. In this respect, their association with 
(space-like) magnetic effects makes it more meaningful to do an effectively 3D averaging, 
viz. Il M n„ =^ |ll 2 ((5 M ^ — 17^17,,) where r) M is a unit vector whose direction need not be 
specified too precisely. After this step, the tracing process is straightforward, and we 
omit the details. But a useful formula is 

Tr[^\ a G;„g s Z g <j a pF a p] = ^< g 2 s G 2 > F, v . (8.34) 

The results for the three quantities Xi K i C are [30]: 

x = -3.56 GeV~ 2 ; k = -0.11; C = +0.06 GeV~ 2 ; (8.35) 
where the QCD-SR value for x is (6 ± 2)GeV- 2 [82]. 

8.4 Axial Condensates 

So far there has been no explicit need to subtract the perturbative contribution (m = 0) to 
the condensates calculated above, since their traces are zero. We now consider the axial 
condensate [0{— « 7m7 5) in a constant external axial field A^, where an explicit subtraction 
is necessary to ensure convergence of the integral. This condensate is connected with the 
axial isoscalar coupling which enters the Bjorken sum rule [84] for DIS of polarized 
electrons on a polarized proton [85] . It is defined through the relation 

< WlplfSl >a= AsA^, (8.36) 
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and its value was calculated in [85b] as f 2 ~ f 2 , on the assumption that the axial field 
interacts with the 8 th component (isoscalar) of the unitary octet current. In the present 
treatment it does not need any such extra assumption but can be simply calculated 
from eq.(8.1) with {Oi = 27^75), and introducing the axial field by the gauge substitution 
— > ILj — 75 A M in the propagator, and keeping only the first order term in the expansion. 
The result is expressed by 



—1 



AsA " = j2^ Tr I d ^ Sp W *^75 S F (U) * 7m75 ] - [«m = 0"] (8.37) 



where the term under quotes is the value of the main term for rh — 0. Evaluating the 
trace and using the isotropy condition < ILjII,, >=5 flu Il 2 /4 we obtain 

A - ^ Tr fd^ r 2 - n2/2 + ^ + n2 / 2 ~^ l (* 

As ~ j2^ Tr J d n[ (A-z g y + (rp-s^ ] ( 8 - 38 ) 

In this case however it is perhaps not as meaningful to keep track of the £ s -terms for 
numerical purposes as for the e.m. case; we shall drop them at this stage. Then with a 
simple rearrangement rh 2 — Il 2 /2= 3II 2 /2 — A/2, the A/2 term can be combined with the 
last term through a Feynman variable u (0 < u < 1) and the pole integration carried out. 
The final result is 



4vr 2 



/ U 2 dU / du fn 2 s + (8.39) 

Jo Jo L (^ 2 + n 2 ) 3/2 (m 2 M + n 2 ) 3 / 2j 



which yields 0.021 GeV 2 , to be compared with f 2 0.018, or perhaps better with f 2 
which is the relevant isoscalar quantity [83c] having a larger value [13] than f 2 . 

For a discussion of these results vis-a-vis QCD-SR, see [30]. Since the spectroscopic 
linkage of the QCD condensates has been the main theme of this Section, we should like to 
end it with the remark that the (MYTP- governed [15]) CIA [16] by itself does not carry 
information on the dynamics of spectroscopy which must be governed by other considera- 
tions (non- perturbative QCD simulated by DBxS [4,24]), but it certainly offers a broad 
enough framework to accommodate such dynamics, without having to look elsewhere. Of 
course, the importance of spectroscopy as an integral part of any 'dynamical equation 
based' approach merely reiterates a philosophy initiated long ago by Feynman et al 



9 qqq Dynamics: General Aspects 

The dynamics of baryons as qqq systems represents the third stage of the three-body 
problem in its journey from the atomic through nuclear to the hadronic level of compos- 
iteness. The first (atomic) stage had been relatively free of theoretical ambiguities due 
to its strong QED foundations in the domain of non-relativistic quantum mechanics. In 
contrast, the second (nuclear) stage, although providing the initial stimulus for few-body 
dynamics, has from the outset remained bogged down in a continual empiricity in the 
theoretical foundations of strong interaction dynamics. Indeed by the time the meson 
exchange picture started being taken seriously for a parallel treatment of meson-nucleon 
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systems on the lines of electron-photon systems, the carpet got quietly removed from 
under its feet, through the slow but sure realization of its tenuous character born out of 
the quark compositeness of the underlying (meson) fields. Indeed the quark-gluon picture 
which had taken firm shape by the end of the Seventies, told in no uncertain terms the 
futility of understanding the inter-hadronic forces directly in terms of their it own species, 
as if they were elementary fields! On the other hand, the emergence of nuclear three-body 
techniques in the Sixties had an instant impact on the quark-level 3-body problem, thus 
providing a big boost to its development in a language strongly reminiscient of the nuclear 
3-body problem, on the lines of Bethe's Second Principle Theory (see Seal.), except for 
the realization of its relativistic character which demands the input dynamics to to be 
Bethe-Salpeter-like (albeit with wide variations), rather than Schroedinger-like. In this 
Section we shall give a panoramic view of three general aspects governing the dynamics 
of qqq baryons: i) classification of baryonic states [86]; ii) problem of connectedness in 
3-body dynamics [87]; iii) BS-dynamics for fermionic qqq systems under DBxS, [29b], in 
parallel with qq, Sec. 4. 3. The details of topics (ii) and (iii) are taken up in Sections 10 
and 11 respectively. 

Yet another type of approach to the qqq problem, as available in the literature, con- 
cerns parametric representations attuned to effective Lagrangians for hadronic transitions 
to "constituent" quarks, with ad hoc assumptions on the hadron-ggg form factor [88a], sim- 
ilar (parametric) ansatze for the hadron- quark-diquark form factor [88b]; or more often 
simply direct gaussian parametrizations for the qqq wave functions as the starting point 
of the investigation [88c]. Such approaches are often quite effective for the investigations 
of some well-defined sectors of hadron physics with quark degrees of freedom, but are in 
general much less predictive than dynamical-equation-based methods like NJL- Faddeev 
[89] or BSE-SDE framework [9-12], when extended beyond their immediate domains of 
applicability. 



The initial qqq formulation was provided by a non-relativistic form of dynamics, and the 
first systematic classification [86a,b] of qqq states proved remarkably successful for the 
understanding of many details of hadronic spectra. On the other hand, the high degree 
of degeneracy of the h.o. model [86a] caused problems on the details of observed states, 
such as the absence of (the relatively low-lying) 20 states, in favour of more restricted 
types which, in a broad SU(6) x 0(3) classification, are all 'natural parity' states [86c] 



while the (complementary) 'unnatural parity' states like 20 + seemed to be missing from 
the data [13]! The natural parity baryons in turn are amenable to a simple quark- 
diquark picture [86d], with diquarks of the types 'scalar-isoscalar' D s and ' (axial) vector- 
isovector'/}^, as well as (complementary) diquarks of the types (pseudo)vector-isoscalar 
and scalar-isovector [90a], all of which go to make up the list (9.1) above. On the 
other hand the 'unnatural' parity baryons require diquark ingredients of opposite parity 
to above, viz., pseudoscalar-isoscalar, vector-isovector; vector-isoscalar, and pseudoscalar- 
isovector, respectively, to make up a complementary list of SU(Q) x 0(3) baryons [29b]: 



9.1 SU(Q) <g> 0(3) Classification 



[56, (even) + ],[70, (odd)-}; 



[70, (even) + ],[56, (odd)') 



(9.1) 



[20, (even) + ), [70, (even)-]; 



[70, (odd) + ],[56, (even)-) 



(9.2) 
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which have not yet been observed [13, 29b]. 

Despite the compactness and elegance of the quark-diquark description, a certain 
amount of dynamical 3-body information gets lost due to the 'freezing' of a quark d.o.f. 
in the (rigid) diquark structure. While a good part of the S3 (permutation) symmetry can 
be recovered by appropriate SU(6) classification [86d], the dynamical information in the 
full 3-body structure is not entirely retrievable, showing up, e.g., in the k 2 - dependence 
of the e.m. form factor of the qqq baryon. To see more clearly the interconnection 
between the two descriptions, let us write down the baryon wave function, with proper 
S'3-symmetry, in both the qqq and q — d q notations. To that end, its full wave function 
with S3- symmetry for 3 identical quarks, may be analyzed into its orbital ip a , spin x an d 
isospin (p components, where a, (3, 7 superscripts stand for the S 3 -symmetry types [25, 
45], which, in the Verde notation [45a], are (s; to', to"; a) for symmetric, mixed-symmetric, 
and antisymmetric respectively. Since only color singlet states are being considered, we 
suppress the (antisymmetric) color wave function C a for brevity, so that the 'active' part 
of the wave function is symmetric [86b]. The full structures of \I/ for 56 and 70 states are 
[90a] 

N\ 6 = (xV + x"0")^7v / 2; N« TO = (0V' + «")x7V2; (9.3) 

N d 70 = [( x >"+x"W + (xV-x"0")^"]/2; 

A*56 = XW; A d 7O = ( x V + x >")07v / 2; 

where the superscripts d and q stand for spin-doublet and spin-quartet respectively, and 
the product of the orbital (ip) and spin (x) functions for higher L-states must be read in 
the standard sense of adding angular momenta in terms of C.G coefficients. For strange 
baryon (A,E) states, the symmetry is reduced to S2, due to the higher mass of the s-quark, 
and the corresponding states have the following representations [90] 

A 56 = <t>W/V2; A 7O = , (xV'±X> , )/2 (9.4) 
£ 56 = fff/v^l £ 7O = 0"( X y±r>")/2 

To relate these structures to the quark(g)-diquark(/}) description, the qD contents 
of these wave functions in a lorentz-invariant form may be read off from the following 
correspondence [90a] 

X '<f>' & D s ; X "0"^i 757 ^r a ; (9.5) 

Here for simplicity, a basis spinor symbol \I/ on the RHS has been suppressed for all 
(baryon) states. However, the additional (Rarita- Schwinger) spin and isospin symbols 
needed for several such states to make up the full baryon structure have been supplied 
via the unit vectors and e M where necessary. [Of course orbital functions ip are needed 
to make up the spatial overlap for the g-D-pair]. This correspondence may be faithfully 
substituted in the set (9.3) to give the precise qD content in S77(6)-form to give the dif- 
ferent cases with correct normalizations. The dynamical effects are now entirely contained 
in the orbital wave function ifj. 
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9.2 Connectedness in a 3-Body Amplitude 

The problem of connectedness [88] in a 3-particle amplitude has been in the forefront of 
few-body dynamics since Faddeev's classic paper [87a] showed the proper perspective, by 
emphasizing the role of the 2-body T-matrix as a powerful tool for achieving the goal. 
The initial stimulus in this regard came from the separable assumption due to Mitra [87b] 
which provided a very simple realization of such connectedness via Faddeev's T-matrix 
structure, a result that was given a firmer basis by Lovelace [87c]. An alternative strategy 
for connectedness in a more general n-body amplitude was provcided by Weinberg [87d] 
through graphical equations which brought out the relative roles of the T-and ^-matrices 
in a more transparent manner. (In particular Weinberg showed that the T-matrix was 
not the only way to achieve connectedness). It was emphasized by both Weinberg and 
Lovelace [87c-d] that an important signal for connectedness in the 3-body (or n- body) 
amplitude is the absence of any 5-function in its structure, either explicitly or through 
its defining equation. This signal is valid irrespective of whether or not the V- or the T- 
matrix is employed for the said dynamical equation. 

The above equations were found for a non-relativistic n- body problem within a basi- 
cally 3D framework [88d] whose prototype dynamics is the Schroedinger equation. For the 
corresponding relativistic problem whose typical dynamics may be taken as the Bethe- 
Salpeter Equation (BSE) with pairwise kernels within a 4D framework, it should be 
possible in principle to follow a logic similar to Weinberg's, using the language of Green's 
functions with corresponding diagrammatic representations [88d] , leading to equations free 
from 5-functions. However there are other physical issues associated with a 4D support 
to the BSE kernel of a confining type, such as contradictions of the spectral predictions 
[14] with data [13]. Indeed, this very issue has been discussed in detail in Section 1, 
culminating in the 'two-tier' 3D-4D BSE approach as the central theme of this article, 
under the name of Covariant Instantaneity [16] for 3D support to the BSE kernel [17] 
which receives formal justification from the MYTP principle [15]. The principal result 
of this ansatz is the exact interconnection between the 3D and 4D forms of the BSE, at 
least for the 4D two-body problem [16]. 

One may now ask: Does a similar interconnection exist in the corresponding BS am- 
plitudes for a three-body system under the same conditions of 3D support for the pairwise 
BS kernel? The question is of great practical value since the 3D reduction of the 4D 
BSE already provides a fully connected integral equation [29b], leading to an approximate 
analytic solution (in gaussian form) [29b] for the corresponding 3D wave function, as a 
byproduct of its success on the baryon spectra [13]. Therefore a reconstruction of the 4D 
qqq wave function in terms of the corresponding 3D quantities should open up a vista 
of applications to various types of transition amplitudes involving qqq baryons, just as in 
the two- body case outlined in Section 4. This exercise is carried out in Section 10, using 
Green's function techniques for both the 2- and 3- body systems (the former for checking 
against the known results of Section 4). There is however a big difference between the two 
systems, born out of the 'truncation' of the Hilbert space due to the 3D support ansatz for 
the pairwise BSE kernel. Such truncation, while still allowing an unambiguous reduction 
of the BSE from the 4D to the 3D level, nevertheless leaves an 'undetermined element' in 
the reverse direction, viz., from 3D to 4D. This limitation for the reverse direction is quite 
general for any n-body system with n > 2; the only exception is the case of n = 2 where 
both transitions are reversible without extra assumptions (a sort of degenerate situation). 
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As will be shown in Sec. 10, the extra assumption (in its simplest form) needed to complete 
the reverse transition is facilitated by some ID 5-function which however has nothing to 
do with connectedness [87]. 

9.3 Fermionic qqq BSE with DB X S 

We now outline the essential logic of a BSE treatment for a fermionic qqq system, for 
pairwise kernels with covariant 3D support, under conditions of DBxS, on closely parallel 
lines to the qq case (Section 4). In Section 3, the derivation [52] of the equation of motion 
from an input Lagrangian for extended 4-fermion coupling shows that the BSE structure 
(4.2) emerges in the linear approximation to the 0- field. This immediately suggests that 
the BSE for a qqq system in the same (linear) approximation must be one with a linear 
sum over all the three pairs of interaction, which for spinless quarks reads as [10b]: 

t(2n) 4 <S>(p lP2 p 3 ) = J2 Afi A F2 [ d q 12 K(q 12 , q' 12 Mp' lP ' 2 p 3 ) (9.6) 

123 J 

where Api——iA i ~ 1 , etc. Under CIA [16], the relative momenta 'hatted', i.e., 

they are orthogonal to the total momenta Pjj of the ij pairs, as explained in Sec. 4. It 
is however more convenient for calculational purposes to take all these relative 'hatted' 
momenta qij to be perpendicular to a common 4- momentum P = pi + p 2 + p 3 , instead 
of to the individual pairs. Technically this amounts to the introduction of 3-body forces 
at the quark level. However the difference turns out to be small [10b], which suggests 
that the 3-body forces are (expectedly) small. Having thus checked this (small) 3-body 
effect, we shall from now on consider a common 'hat' symbol for all the three pairs, i.e., 
Qij;n—Qij — Qij-PPfj./ P 2 '■ It is this version for hatted symbols for the qqq problem that we 
shall consider in the next two sections. 

10 Interlinking 3D And 4D qqq Vertex Fns 

In this Section we outline a fairly detailed (self-contained) method of Green's functions 
for 2- and 3- particle scattering near the bound state pole, for the 3D-4D interconnection 
between the corresponding wave functions. For simplicity we consider identical spinless 
bosons, with pairwise BS kernels under CIA conditions [16], first for the 2-body case for 
calibration, (see Sects. 4. 1-2), and then for the corresponding 3-body case, on the basis of 
the Green's Fn counterpart of the general structure, eq.(9.6). 

10.1 Two-Quark Green's Function Under CIA 

In the notation and phase convention of Section 4, the 4D qq Green's fn G(pip 2 ;pi'p 2 ) 
near a bound state satisfies a 4D BSE (no inhomogeneous term): 

t(27r) 4 G(p 1 p 2 ;p 1 'p 2 ') = A^A^ 1 j dp 1 "dp 2 "K(p 1 p 2 ;p 1 "p 2 , ')G(p 1 "p 2 ";Pi'P2) (10.1.1) 
where 

A 1=Pl 2 + m g 2 , (10.1.2) 
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and m q is the mass of each quark. Now using the relative 4- momentum q= ij>\ — P2}/2 
and total 4-momentum P = p 1 + p 2 (similarly for the other sets), and removing a 5- 
function for overall 4-momentum conservation, from each of the G- and K- functions, 
eq.(lO.l.l) reduces to the simpler form 

i{2TifG{q.q') = A 1 - 1 A 2 - 1 J dq"Mda"K(q, q")G(q", q') (10.1.3) 

where g M = g M — <rP M , with a = (q.P)/P 2 , is effectively 3D in content (being orthogonal 
to P M ). Here we have incorporated the ansatz of a 3D support for the kernel K (inde- 
pendent of o and cr'), and broken up the 4D measure dq" arising from (10.1.1) into the 
product dq"Mda" of a 3D and a ID measure respectively. We have also suppressed the 
4-momentum P M label, with (P 2 = — M 2 ), in the notation for G(q.q'). 
Now define the fully 3D Green's function G(q, q') as [47] 

G(q,q') = J J M 2 dada'G(q,q') (10.1.4) 
and two (hybrid) 3D-4D Green's functions G(q,q'), G(q,q') as 

G(q, q') = J MdaG(q, q'); G(q, q') = J Mda'G(q, q'); (10.1.5) 
Next, use (10.1.5) in (10.1.3) to give 

i (2n) 4 G(q, q') = A 1 _1 A 2 _1 j dq"K(q, q")G(q", q') (10.1.6) 
Now integrate both sides of (10.1.3) w.r.t. Mda and use the result 

J MdaA l - 1 A 2 - 1 = 2mD- 1 (q); D(q) = 4uj(lo 2 - M 2 /4); 0? = m q 2 + q 2 (10.1.7) 
to give a 3D BSE w.r.t. the variable q, while keeping the other variable q' in a 4D form: 

(2vr) 3 G(g, q') = IT 1 J dq"K(q, q")G(q", q') (10.1.8) 

A comparison of (10.1.3) with (10.1.8) gives the desired connection between the full 4D 
G-function and the hybrid G(q, g')-function: 

2niG(q, q') = D(q)A 1 - 1 A 2 - 1 G(q, q') (10.1.9) 

Again, the symmetry of the left hand side of (10.1.9) w.r.t. q and q' allows rewriting the 
right hand side with the roles of q and q' interchanged. This gives the dual form 

2mG{q, q') = DfflA^A^Giq, q') (10.1.10) 

which on integrating both sides w.r.t. Mda gives 

2mG(q, q') = D(q') A^ 1 A^G^, q'). (10.1.11) 

Substitution of (10.1.11) in (10.1.9) then gives the symmetrical form 

(2m) 2 G(q, q') = D^A^A^G^q, q')D(q') A/" 1 A^ 1 (10.1.12) 



53 



Finally, integrating both sides of (10.1.8) w.r.t. Mda', we obtain a fully reduced 3D BSE 
for the 3D Green's function: 

(27r) 3 G(g, q') = D-\q J dq"K(q,q")G(q" ,q') (10.1.13) 

Eq. (10.1.12) which is valid near the bound state pole, expresses the desired connection 
between the 3D and 4D forms of the Green's functions; and eq(10.1.13) is the determining 
equation for the 3D form. A spectral analysis can now be made for either of the 3D or 
4D Green's functions in the standard manner, viz., 

G(q, q') = £ $ n (g; P)<S>* n (q'; P)/(P 2 + M 2 ) (10.1.14) 

n 

where $ is the 4D BS wave function. A similar expansion holds for the 3D G-function G 
in terms of ra (g). Substituting these expansions in (10.1.12), one immediately sees the 
connection between the 3D and 4D wave functions in the form: 

27ri$(g, P) = A 1 - 1 A 2 - 1 D(q)(P(q) (10.1.15) 

whence the BS vertex function becomes T = D x (f>/{2m) as found in [16]. We shall make 
free use of these results, taken as qq subsystems, for our study of the qqq G-functions in 
subsections 2 and 3. 



10.2 3D BSE Reduction for qqq G-fn 

As in the two-body case, and in an obvious notation for various 4-momenta (without 
the Greek suffixes), we consider the most general Green's function G{p\p 2 p3\p\p 2 p3) f° r 
3-quark scattering near the bound state pole (for simplicity) which allows us to drop the 
various inhomogeneous terms from the beginning. Again we take out an overall delta 
function 5(p\ + p 2 + P3 — P) from the G-function and work with two internal 4-momenta 
for each of the initial and final states defined as follows [10b]: 

V3£ 3 = p! - p 2 ; 3?73 = -2p 3 + pi + p 2 (10.2.1) 

P = Pl+P2+P3 = Pl'+P2+P3 (10.2.2) 

and two other sets £i,?7i and £2,^2 defined by cyclic permutations from (10.2.1). Fur- 
ther, as we shall consider pairwise kernels with 3D support, we define the effectively 3D 

momenta pi, as well as the three (cyclic) sets of internal momenta (i = 1,2,3) by 
[10b]: 

Pi=Pi- ViP ; ii = ii - SiP ; f]i - UP (10.2.3) 

Ui = (P Pt )/P 2 ; Si = (P.ii)/P 2 ; U = (P.<m)/P 2 (10.2.4) 

V3s3 = v\ — v 2 ; 3t3 = — IV3 + v\ + v 2 (+cyclicpermutations) (10.2.5) 
The space-like momenta pi and the time-like ones Vi satisfy [10b] 

Pi + P2 + P3 = ; vx + v 2 + 1/3 = 1 (10.2.6) 

Strictly speaking, in the spirit of covariant instantaneity, we should have taken the relative 
3D momenta £,17 to be in the instantaneous frames of the concerned pairs, i.e., w.r.t. the 
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rest frames of = Pi + Pj ; however the difference between the rest frames of P and Pij 
is small and calculable [10b], while the use of a common 3-body rest frame (P = 0) lends 
considerable simplicity and elegance to the formalism. 

We may now use the foregoing considerations to write down the BSE for the 6-point 
Green's function in terms of relative momenta, on closely parallel lines to the 2-body case. 
To that end note that the 2-body relative momenta are qij = (pi—pj)/2 = y/3£ k /2, where 
(ijk) are cyclic permutations of (123). Then for the reduced qqq Green's function, when 
the last interaction was in the (ij) pair, we may use the notation G(^kVk]Ck'i]k'), together 
with 'hat' notations on these 4-momenta when the corresponding time-like components 
are integrated out. Further, since the pair £k,Vk is permutation invariant as a whole, 
we may choose to drop the index notation from the complete G-function to emphasize 
this symmetry as and when needed. The G-function for the qqq system satisfies, in the 
neighbourhood of the bound state pole, the following (homogeneous) 4D BSE for pairwise 
qq kernels with 3D support: 

*(2tt) 4 G(^; grf) = £ Af^f 1 / dq'{ 2 Mda l2 " K{q 12 , # 2 )G(& V; &V) (10.2.7) 

123 

where we have employed a mixed notation (q 12 versus £ 3 ) to stress the two-body nature 
of the interaction with one spectator at a time, in a normalization directly comparable 
with eq.(10.1.3) for the corresponding two-body problem. Note also the connections 
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V3s 3 /2; gi2 = a/36/2; m = ~P3, etc (10.2.8) 



The next task is to reduce the 4D BSE (10.2.7) to a fully 3D form through a sequence 
of integrations w.r.t. the time-like momenta Sj, U applied to the different terms on the 
right hand side, provided both variables are simultaneously permuted. We now define the 
following fully 3D as well as mixed (hybrid) 3D-4D G-functions according as one or more 
of the time-like £,r) variables are integrated out: 



G(Zm?v')= I I I I dsdtds'd£G(&n; £'rf) (10.2.9) 

which is ^-symmetric. 

G^mi'V) = I I dt 3 dt 3 'G(^;ev'); (10.2.10) 



GattimN) = JI dssdss'G^ev'); (10.2.11) 



The last two equations are however not symmetric w.r.t. the permutation group S3, since 
both the variables £, rj are not simultaneously transformed; this fact has been indicated in 
eqs. (10.2. 10-11) by the suffix "3" on the corresponding (hybrid) G-functions, to emphasize 
that the 'asymmetry' is w.r.t. the index "3". We shall term such quantities "SVindexed" , 
to distinguish them from S'3-symmetric quantities as in eq. (10.2.9). The full 3D BSE 
for the G-function is obtained by integrating out both sides of (10.2.7) w.r.t. the st- 
pair variables dsids/dUdt/ (giving rise to an SVsymmetric quantity), and using (10.2.9) 
together with (10.2.8) as follows: 

WG{M'V') = J2D-\q l2 ) f d$l 2 K{q 12 Al 2 )G{iY-,m (10.2.12) 

123 
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This integral equation for G which is the 3-body counterpart of (10.1.13) for a qq system 
in the neighbourhood of the bound state pole, is the desired 3D BSE for the qqq system in 
a fully connected form, i.e., free from delta functions. Now using a spectral decomposition 
for G 

G(imi'v') = J2Miv;Pmi'v';P)/(P 2 + M 2 ) (10.2.13) 

n 

on both sides of (10.2.12) and equating the residues near a given pole P 2 = — M 2 , gives 
the desired equation for the 3D wave function <p for the bound state in the connected 
form: 

(2tt) 3 0(^;P) = J2D~\qn) [ dq'{ 2 K{q 12 A'i 2 )4>{i"v"\P) (10.2.14) 

123 J 

Now the ^-symmetry of <fi in the (£j, fji) pair is a very useful result for both the solution 
of (10.2.14) and for the reconstruction of the 4D BS wave function in terms of the 3D 
wave function (10.2.14), as is done in the subsection below. 

10.3 Reconstruction of 4D qqq Wave Function 

We now attempt to re-express the 4D G-function given by (10.2.7) in terms of the 3D 
G-function given by (10.2.12), as the qqq counterpart of the qq results (10.1.12-13). To 
that end we adapt the result (10.1.12) to the hybrid Green's function of the (12) sub- 
system given by G 3v , eq.(10.2.10), in which the 3-momenta fj 3 ,fj 3 play a parametric 
role reflecting the spectator status of quark #3, while the active roles are played by 
Q12, Q12 = V3(^3, C3O/2, fc> r which the analysis of subsec.(lO.l) applies directly. This gives 

(2m) 2 G 3v (^f l3 ; = D(q 12 ) A 1 ~ 1 A 2 ~ 1 G(i 3 f, 3 ; &ft WiaJV^A.,'- 1 (10.3.1) 

where on the right hand side, the 'hatted' G-function has full S^-symmetry, although (for 
purposes of book-keeping) we have not shown this fact explicitly by deleting the suffix 
'3' from its arguments. A second relation of this kind may be obtained from (10.2.7) by 
noting that the 3 terms on its right hand side may be expressed in terms of the hybrid 
G 3 £ functions vide their definitions (10.2.11), together with the 2-body interconnection 
between (£3, £3') and (£3, £3) expressed once again via (10.3.1), but without the 'hats' on 
n 3 and rjs- This gives 

(v^) 2 G(fcife;&V) = (v^) 2 G(ta£V) 

= ^Ar 1 A 2 - 1 ( 7 riv / 3) I d$l 2 Mda l2 "K{q l2 Al 2 )G{U'riz\Uriz) 

123 J 

= £ D ($i2) Ar 1 A 2 - 1 G 3e (£ 3 7 73 ; £3773') Ar'Ar 1 (10.3.2) 

123 

where the second form exploits the symmetry between £,77 and £',77'. 

At this stage, unlike the 2-body case, the reconstruction of the 4D Green's function 
is not yet complete for the 3-body case, as eq. (10.3.2) clearly shows. This is due to the 
truncation of Hilbert space implied in the ansatz of 3D support to the pairwise BSE 
kernel K which, while facilitating a 4D to 3D BSE reduction without extra charge, does 
not have the complete information to permit the reverse transition (3D to 4D) without 
additional assumptions. The physical reasons for the 3D ansatz for the BSE kernel have 
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been discussed in detail elsewhere [47], vis-a-vis contemporary approaches. Here we look 
upon this "inverse" problem as a purely mathematical one. 

We must now look for a suitable ansatz for G 3 ^ on the right hand side of (10.3.2) 
in terms of known quantities, so that the reconstructed 4D G-function satisfies the 3D 
equation (10.2.12) exactly, as a "check-point" for the entire exercise. We therefore seek a 
structure of the form 

Gac&i/s; £V) = G&rk; isVs) x F( P3 ,p 3 ') (10.3.3) 

where the unknown function F must involve only the momentum of the spectator quark 
#3. A part of the 773,773' dependence has been absorbed in the G function on the right, 
so as to satisfy the requirements of S^-symmetry for this 3D quantity [47]. 

As to the remaining factor F, it is necessary to choose its form in a careful manner so 
as to conform to the conservation of 4-momentum for the free propagation of the spectator 
between two neighbouring vertices, consistently with the symmetry between p 3 and p 3 . 
A possible choice consistent with these conditions is the form : 

F(p 3 ,P3) = CsAs -1 ^ - (10.3.4) 

Here A 3 _1 represents the "free" propagation of quark #3 between successive vertices, 
while C 3 represents some residual effects which may at most depend on the 3-momentum 
p 3 , but must satisfy the main constraint that the 3D BSE, (10.2.12), be explicitly satisfied. 

To check the self-consistency of the ansatz (10.3.4), integrate both sides of (10.3.2) 
w.r.t. ds 3 ds 3 dt 3 dt 3 to recover the 3D ^-invariant G-function on the left hand side. Next, 
in the first form on the right hand side, integrate w.r.t. ds 3 ds 3 ' on the G-function which 
alone involves these variables. This yields the quantity G 3 ^. At this stage, employ the 
ansatz (10.3.4) to integrate over dt 3 dt 3 . Consistency with the 3D BSE, eq.(10.2.12), now 
demands 

C 3 J J dv 3 dv 3 '/\ 3 - l 5(v 3 - v 3 ') = 1; (sincedt = dv) (10.3.5) 

The ID integration w.r.t. dv 3 may be evaluated as a contour integral over the propagator 
A -1 , which gives the pole at v 3 = u 3 /M, (see below for its definition). Evaluating the 
residue then gives 

C 3 = iit/{MCo 3 )- Col = m q 2 + p\ (10.3.6) 

which will reproduce the 3D BSE, eq. (10.2. 12), exactly] Substitution of (10.3.4) in the 
second form of (10.3.2) finally gives the desired 3-body generalization of (10.1.12) in the 
form 

3G(^;£V) = E^(^i2)AiFA 2F /J(g; 2 )A 1F / A 2F 'G(^3;&3)[A3F/(M7rcu 3 )] (10.3.7) 

123 

where for each index, A^ = — iA" 1 is the Feynman propagator. 

To find the effect of the ansatz (10.3.4) on the 4D BS wave function $(^77; P), we do 
a spectral reduction like (10.2.13) for the 4D Green's function G on the left hand side of 
(10.3.2). Equating the residues on both sides gives the desired 4D-3D connection between 
$ and 0: 

P ) = E D(qu) Ai-'A.-Vfe P) x ^ ( ^~^ M) (10.3.8) 
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defines the 4D wave fn in terms of piecewise vertex fns V*, as 



From (10.3.8-9), we infer the baryon-ggg vertex function V3 corresponding to the 'last' 
interaction in the 12-pair as 

V 3 = D(q 12( f>(i,f,) x v /2A 3 5(z/ 3 2 M 2 - <%) (10.3.10) 

and so on cyclically. (The argument of the (^-function inside the radical for V3 simplifies 
to p\ + m 2 ,). This expression is essentially the same as eq.(5.15) of ref. [10b], which had 
been obtained from largely intuitive considerations. 

To account for the appearance of the ID S-in under radical in (10.3.10), it is explained 
elsewhere [47] that it has nothing to do with connectedness [88] as such, but merely reflects 
a 'dimensional mismatch' due to the 3D nature of the pairwise kernel K [16] imbedded in 
a 4D Hilbert space. This in turn is the result of the 'contact' nature (in time dimension) 
of the pairwise interaction, somewhat analogous to a Fermi 5-fn potential to simulate 
the effect of the (short range) nuclear n — p interaction in the 'molecular' problem of 
(specular) neutron scattering by a hydrogen molecule [91]. As a further self-consistency 
check, it is instructive to compare (10.3.10) with one obtained by taking the limit of 
a point interaction, which amounts to setting K = Constant in the entire derivation 
above. This structure [47] which is worked out in Appendix C, is free from radicals, and 
explicitly 4D-invariant, in agreement with the so-called NJL-Faddeev (contact [4]) model 
[92] of 3-particle scattering. 



11 Fermion Quarks: QCD-Motivated qqq BSE 

We now turn to the more realistic case of fermion quarks for which we shall draw freely 
from a relatively recent analysis [29b] of a qqq baryon, which is basically a 3-body gener- 
alization of subsection 4.3 for the two-body case. For simplicity of description, without 
sacrificing the essential physics, we shall specialize to equal mass kinematics (mass=m g ). 

11.1 3D Reduction of 4D qqq BSE 

The starting 4D BSE has the form (c.f. (9.6)): 

(2n) 4 ^(p 1 p 2 p 3 ) = iS F ( Pl )S F (p 2 ) ]T / d 4 q' 12 K(q 12 , q' 12 M Pl p 2 Ps) (H-l) 

123 J 

where the kernels are given by eqs.(4.18-19) for each ij pair, except for the Casimir 
value of the color factor F 12 = Xi.X 2 /4 for the 3 state of a qq pair (to produce a color- 
singlet baryon), which is just half its 'singlet' value for a qq pair. And of course the DB\S 
mechanism is built-in as in the two-body case of Sec. 4. 3. The Gordon reduction of the 
product of two 7^-matrices [10] also goes through as in Sec. 4. 3, leading to [29b]: 

HpiP2 P3 ) = E , 9 \\ a 2 a / d% 2 VMvWV(hatq 12 , ^ 12 $(Mps) (H-2) 
123 (2tt) 4 AiA 2 J 
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where, following the steps of Sec. 4. 3, the 'bosonic' $-fn is related to the fermionic ^-fn, 
as in eq.(4.24), by [10b,29b]: 

*(p0 = U^S^i-p^ip,); A t = m\ + p 2 (11.3) 

while the 4- vectors V® are given by eq.(4.25). 

Next, for the 3D reduction of eq.(11.2), we need to define the transverse pi and longi- 
tudinal Vi components of the 4-momenta pi, as in Sec. 10.2, eqs.( ), and multiply the pairs 
of V^-fns, as in Sec. 4. 3, replacing in the process the longitudinal components Vi by their 
on- shell values Ui/M, where Cjf = m 2 +p 2 [29b], uniformly from such products [93a]. Now 
define the 3D wave function ip in the as in Sec. 10.2, viz., [10b, 29b]: 

^(PmPs) = J ds i dt i ®(p 1 p 2 p3); V3s 3 = u 1 -u 2 ; 3t 3 = -2v 3 + v x + v 2 (11.4) 

The product dsidti is cyclically invariant, so that the definition (11.4) can be taken over for 
all the three terms on the RHS of (11.2), with proper indexing. The rest of the procedure 
is straightforward, and follows closely the pattern laid out in the original formulations. 
Thus one integrates both sides of (11.2) w.r.t. ds 3 dt 3 , making use of (11.4) as well as 
the measure d A q' l2 = d 3 q' 12 Mds' 3 \/3/2 to give on its RHS / ds' 3 dt 3 <&' = i^{p'ip' 2 p 3 ). The 
additional 0^3-integration on the RHS is expressed by the result [10b]: 

^/^ = Du = -n u \[(M\l-v 3 )^uW 2 ]/2M\l-v 3 f (11.5) 

2/Q 12 = flu/uJ! + /*2lMl Al2;21 = ± (H-6) 

where v 3 has its on-shell value uj 3 /M in the foregoing equations, as befits a spectator 
quark in the first term of(11.2). The resultant 3D reduction of (11.2) now takes the form: 

HPmPs) = E , F , 3n I d% 2 V^.V^V(q 12 , q^ip'^p,) (11.7) 

123 <A) D 12 J 

11.2 Reduction to 6D Harmonic Basis 

The next task is to reduce eq.(11.7) to a more transparent form suitable for numerical 
treatment. To that end we base our procedure [29b] on the expected smallness of the 
S'3-invariant quantity 5= M — uj 1 — u 2 — uj 3 compared to Ui and/or M. This gives the 
crucial result: 

D l2 = -Aujiuj 2 5 + (5 2 ); -5 = u x + oj 2 + uj 3 - M (11.8) 

which ensures that in (11.7), all the three terms on its RHS have a common denominator 
S which, when transferred to the LHS, serves as a natural 'energy denominator' for the 
entire qqq equation. [Since the terms of 0(5 2 ) in (11.8) are fully calculable, any effect 
on their omission can be estimated perturbatively if necessary]. Next, from Sec4.3, the 
confining part of V(q± 2 , q' 12 ) is harmonic for -urf-quarks (Aq = 0), so that a perturbative 
treatment is possible, based on the (harmonic) confining part of V(q, q')\ 

V con = -(27r) 3 u qq 2 [Wl + Cq/^IW^ ~ <?'); ^W 2 = 4M 12 fi 12 fi 21 uj 2 a s {M 2 2 ) (11.9) 
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In this formula, the definitions (11.7) for the fractional momenta /ti 2 ,2i conform to their 
Wightman-Gaerding [56] definitions for unequal mass kinematics, a la eq.(4.1) of Section 
4, since the unequal masses arise from the mass-shifts m q — > Ui of the quarks (1, 2) in 
the presence of the spectator #3. Since such shifts are small, it is fairly accurate to 
approximate the fractional momenta as (1 — z/ 3 )/2 each, while M 12 ~ M — ou^ only. Now 
to emphasize the 3D character of the various momenta, define the pairwise items: 

Ly = -iq>j x Vij-; Vi2 = Vi-V 2 ; 2q 12 = Pi - p 2 ; Qu = 4q 2 2 V? 2 + 8qi 2 .Vi 2 + 6 

(H-10) 

Also to take full advantage of the HO form (11.9), recast the (small) energy denominator 
5 in the alternative form [29b]: 

-2Mb « (cj 1 + w 2 + cu 3 ) 2 -M 2 < 3(ujI + ujI+ujI) = A = 9mJ + 9(^ 2 + r] 2 )/2-M 2 (11.11) 

The resulting 'Master Equation' (11.7) is in pairwise notation [29b]: 

M = (W con +W OGE )^W con = M W 2 ^(l-z/ 3 ) 2 «i2 S M 12 x[V 2 2 + ^ + ^iVCT] (11.12) 

123 ^0 ^1^2 

(DC i 
NC T = ^ - - Li 2 (<ti + a 2 ) + -p 3 x Vi 2 (<Ti - a 2 ) - a x a 2 (11.13) 



and the OGi? term in a 'mixed' (r, p) representation is [29b]: 

Woge = ^E «*2[— + — (Qi2— qi2 + ^ 3 (r 12 )(l - 2^(72/3)) + etc] (11.14) 

d 123 r 12 ^1^2 ?" 12 

The OGE'-term is calculated perturbatively in a 6D HO basis given by the main confining 
term in (11.12) [29b], which in a common (£, rj) basis [29a], now provided by eqs. (10.2. 1-6) 
of Section 10, reads [29b]: 

.2/1 .-^2= -Arov72 , ow2 , C (M 2 — 3m 2 + A) Q B - 8JS + 18 



A^ = M W 2 (1 - uYa s (M - c)[2V 2 + 2V 2 + ^ + ™ 4 _ 2 ]V> 

(11.15) 

where the operators Qb, etc are defined in [29a], and the symbols (a,u,u) with 'bars' 
represent their 'average' values [29b]. From this H.O. equation, the scale parameter (3, 
analogous to the 2-body quantity (4.31) may be inferred as [29b] 

= 4Muj 2 & s (1-u) 2 (M -Q)/9; a s = 1/faJ 1 - 2C M(1 - vfj{M - Q)] (11.16) 

so that the basis function ipo in its ground state is exp[— |(£ 2 + i] 2 )/j3 2 ]; similar functions 
exist for L-excited states [40], providing a basis for perturbative treatment [29a] of the 
OGE terms (11.14). 

11.3 Complex HO Basis for qqq States 

It is however mathematically simpler to convert eq.(11.15) to a complex basis. To this 
end we define the complex dimensionless 3- vectors Zi, z*, and their (derivative) conjugate 
momenta d Zi , as 

V2p[ Zi ; z*} = & ± Vi ; V2(3~ l [d z - d z *\ = % =F d m (11.17) 
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For the construction of angular momenta in complex basis, see Appendix D. 

A more convenient basis for handling the various terms in (11.15) is provided by the 
creation/annihilation operator representation [29b, 46b] defined by two sets of complex 
operators 

V2a i = z t + d z: ; V2a* = z* + d Zi ; V2a\ = z* - d Zi ; yfia? = z t - d* Zi (11.18) 

which satisfy the commutation relations 

[a i ,a]] = [a* i ,a*i] = 6 ij (11.19) 

with all other pairs commuting. In the next subsection 11.4, we define the number oper- 
ators N c and N* which now play the role of N% and N v , but unlike the latter, the former 
can be simultaneously diagonalized; so their sum N and difference N a are both constants 
of motion. Together with certain two-step operators, they form several sets of SO(2, 1) 
algebras (described below), which diagonalize the momentum dependent operators Qb, 
etc, in terms of their respective Casimirs [46b], so that the solution of eq.(11.15) takes a 
simple algegraic form [29b]: 

F(M, N) = F con (M, N) + Foge(M.N) = N + 3 (11.20) 

where the first term is given by Eq.(4.19) of [29b], while the second term lends itself to 
a simple perturbative treatment (see [29b] for details). Appendix D gives a summary of 
the normalized SU(6) x 0(3) structures of the 3D ip-fias in the complex basis, which are 
needed for calculating the F OGE term of eq. (11.20) above. 

11.4 50(2, 1) Algebras of Bilinear Operators 

We start by defining the number operators N c , N*, and the mixed quantities N m) 
[46b] 

N c = a\ ai - N* c =a?a^ N m = a i af; rf m = N* m = a*a\ (11.21) 
and their linear combinations 

N — N c + N* = + jV„; N a = N c - iV c * (1 1.22) 

Note that both N and iV a are simultaneously diagonal in this (complex) representation, 
while in the real (£, rj) basis, only their sum is diagonal. Next define the two-step operators 
(and their h.c.'s) [46b]: 

A = 2a ia *; = 0*0*; C* = a*a* (11.23) 
^ = 2ajaf; C* = ajaj; = 

Now the trio A, A^ and iV form an S^-symmetric set, whose normalized forms 

Q+ = A^/2; Q^ = -A/2; Q 3 = (N + 3)/2 (11.24) 

form an 5*0(2, 1) algebra (bounded from below) with the Casimir [93, 46b] 

u(u + 1) = Q 2 = -(AA ] + A^A)/8 + (N + 3) 2 /4 (11.25) 
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where u(u + 1) = +3/4 for even N and +2 for odd N, while the eigenvalues of Qs are 
— u + k, (with A; = 0, 1, 2, ...). These imply that u = —3/2 and -u = — 2 for even and odd 
N respectively. Similarly the mixed symmetric set (C, C\ N c ) form an 5*0(2, 1) algebra 
in the normalized form [46b] 

Q c+ = C7+/2; Q c - = -C/2; Q c3 = ^(iV c + 3/2) (11.26) 

with the corresponding Casimir [46b] 

u c (u c + 1) = Q 2 C = -(CC* + CtC)/8 + (N c + 3/2) 2 /4 (11.27) 

This spectrum is again bounded from below [93], with the eigenvalues Q c3 = —u c + k, 
where u c = —3/16 for even N c and u c = +5/16 for odd iV c . An identical structure 
holds for the 'starred' operators (C*, C*' , N*), with the same eigenvalues. Finally the trio 
(N m ,Nj n = N^,N a ) which is 5 3 -antimmetric, satisfy a 'normal' SO (3) algebra [46b]: 

[N m , N a ] = 2N m ; [Nl N a ] = -2^; [N m , <] = -N a (11.28) 

with spectra bounded from both above and below. The corresponding Casimir is 

s(s + 1) = (N m Nl + NlN m )/2 + N 2 J A (11.29) 

The spectrum is here determined from the condition that both N c and N* are non-negative 
integers. The result is [46b] 

— N < N a < N; s = N/2 (11.30) 
11.5 "Exotic" qqq States 

The comparison of eq. (11.20) with the baryon spectra is described at some length in [29b], 
and it is not our intention here to go into these details afresh. Instead, we shall end this 
Section with some qualitative analysis on the capacity of this model to identify some 
exotic baryonic states which have for long remained elusive. The main reason for such 
optimism stems from the precise predictions on the 'spectroscopic' locations of the states 
on the one hand, and the possibility of making more reliable SU(Q)<S>0(3) assignments for 
such states on the basis of their decay characteristics which the model also allows within 
its broad framework. To see the logic, a good calibration is first provided by the fairly 
accurate location of several 'known' states in a parameter-free manner; see Table I of [29b] 
for comparison. With this first check, a more sensitive test is now a comparison of the 
alternative SU(6) <E>0(3) assignments for the mass locations of the same states; see Table 
II of (29b). Specifically the competition is between the 56, odd" and 70, odd~ assignments 
for A-like states for the same total quantum number N. The question is clearly of physical 
interest since in the entire history of baryon spectroscopy 56, odd~ states have suffered 
from popular perceptions of elusiveness, despite occasional attempts to the contrary [91]. 
The analysis in [29b] suggests that the 56 assignment has a slight edge over 70, at least 
for a couple of odd-patity states by virtue of 'location', but a more sensitive test requires 
a more detailed comparative study of the decay and/or production characteristics that 
these alternative assignments provide, vis-a-vis the data (which are still elusive). In this 
respect it was shown in [29b] that the general mechanism of 'direct' versus 'recoil' modes 
of single-quark transitions [45b], do not inhibit in any way the production of natural parity 
56~ states w.r.t. the corresponding 70 _ states, (perhaps contrary to popular beliefs). 
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11.6 CIA vs CNPA for Fermionic qqq Dynamics 

In the foregoing we have mostly described the CIA predictions [29b] on the baryon spec- 
tra. How about the corresponding ggg-scenario with the other MYTP-governed CNPA 
dynamics whose qq counterpart has been employed in Sections 4-6 ? The reason for avoid- 
ing this exercise for the qqq problem is one of pedagogical necessity. For, from the results 
of Sections 4-6, it has been fairly clear that the earlier NPA treatment [40] based on the 
old-fashioned NP-language [35] formally provides the same CNPA structure of 3D BSE as 
well as 4D vertex functions for qq systems, so that a similar qqq structure should be ex- 
pected. In this respect, the baryon spectral results [29a] based on the old-fashioned NPA 
treatment are already available in detail [40], and the comparison with the CIA treatment 
[29b] shows considerable overlap therewith. As for the reconstruction of the qqq vertex 
function under CNPA, a closely analogous treatment akin to Section 10 formally leads to 
almost identical results, with the CIA-CNPA correspondence already indicated in Section 
4. 

How about the reconstruction of the 4D qqq vertex function for fermion quarks ? 
Again the treatment, which is analogous for both CIA and CNPA, consists in reducing 
the fermionic structure to an effectively scalar problem via eq.(11.3) which relates the 
fermionic BS wave function ^ to the 'scalar' function $, fits in smoothly with the treat- 
ment outlined in Section 10 for spin zero quarks, with almost no change, thus rendering 
unnecessary another fresh formulation for fermions. As for quark loop applications to 
the qqq problem, the general problem of 'Lorentz mismatch' of 3D wave functions in a 
quark-loop integral, that had led us to abandon the CIA treatment in favour of CNPA 
for the qq problem (see, e.g., Sections 5-6), is also encountered in the qqq case, so that it 
is profitable to adopt CNPA [41] for baryonic transition amplitudes as well. 

The only exceptions are two-loop integrals, as in the self-energy problem (Sect. 7), or 
one-loop integrals, as in the vacuum condensate problem (see Sect. 8), where this pathology 
is just avoided. Full-fledged baryon-loop calculations are still being developed, so such 
topics are not intended for a detailed coverage in this article, except for indicating the 
results of a recent calculation of SU (2) n— p mass splitting [95] analogous to the treatment 
of pseudoscalar mass splittings (Sect. 7) [32a] by this method. Thus, using the same value 
(4MeV) of the 'current' d — u mass difference, the total n — p mass difference works out 
as 1.28MeV [95], to be compared with the experimental value of 1.29MeV [13], except 
for possible QED gauge corrections [63]. On this last item, an indication of the expected 
correction is available from its effect on the Kaon e.m. mass difference, which yields a 
~ 60% upward revision on its (uncorrected) value of about lMeV [32a]; see Appendix 
C for an estimation of this correction. If this value for the kaon case is taken as rough 
indication of the seme effect expected for the nucleon case, then (on a proportionate basis) 
the QED gauge corrected value for n — p mass difference comes down to ~ lMeV. For 
details of this methodology, see [95]. 

12 Summary And Conclusions 

In this article an attempt has been made to present a somewhat 'less than conventional' 
BSE-SDE formalism based on the Markov- Yukawa Transversality Principle (MYTP) [15] 
on the one hand, and a strongly QCD motivated 4-fermion Lagrangian which generates 
the BSE-SDE framework by breaking its chiral symmetry dynamically (DB\S) [23-27], 
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on the other. The MYTP mechanism provides an exact interconnection between the 
3D and 4D forms of the BSE, so that both can be used interchangeably, a facility which 
does not seem to exist in other alternative 3D BSE formalisms [39], or the null-plane 
formulations-both non-covariant [35] and covariant [36-38]. This twin property of the 
MYTP-governed BSE formalism [16], termed 3D-4D BSE for short, gives rise to a natural 
'two-tier' description [40], the 3D sector (with its relativistic Schroedinger-like BSE) being 
appropriate for making contact with the hadron spectra [13], while the reconstructed 4D 
BSE yields a vertex function which allows the direct use of the language of Feynman 
diagrams for evaluating transition amplitudes as 4D loop integrals. (This contrasts with 
other 3D formulations [35-39] which require specialized versions of Feynman diagrams [37] 
for calculating loop integrals). 

At a more quantitative dynamical level, both qq and qqq hadrons are amenable to a 
unified treatment, since their respective BSE's emanate from a common (input) chiral 4- 
fermion Lagrangian with a gluon-like propagator whose 'color-factor' has the right relative 
strengths for both systems. And while the 3D-4D structure of the qq BSE [16,28], as well 
as the 3D reduction of the qqq BSE [29], have been around for some time, the missing 
link of a reconstructed 4D BS wave function for the qqq system (only conjectured in [10b]) 
has now (hopefully) been supplied through a formal derivation in Sect. 10 via Green's 
Function techniques [47]. Indeed the main emphasis in this Article has been on the 
'second stage' of this two-tier formalism, relating to the calculation of 4D quark-loop 
integrals, of which some selected examples have been given in Sects. 5-7, to bring out 
the feasibility of its applications to the meson sector. The corresponding applications 
to baryonic amplitudes via loop integrals are still being developed, and only a 'pilot' 
example, relating to SU(2) mass splittings [95], is as yet available. However the scope 
(and feasibility) of such applications is quite substantial [96]. 

The capacity of this BSE-SDE formalism to relate its parameters to the 'vacuum con- 
densates' of QCD-SR theory [2] has been sought to be brought out in Sect. 8, wherein it 
has been shown that several types of condensates (both direct and induced) lend them- 
selves with great ease to this simple treatment [30], while the corresponding QCD-SR 
treatments [82-85] often need additional ansatze for their evaluation. This facility it owes 
to its (input) gluon propagator on the one hand, and the (derived) mass function m{p) 
from the SDE solution [11] on the other. The two fundamental parameters [11] of the in- 
frared gluon propagator are not calculable within this (Bethe's Second Principle oriented) 
framework, but they are firmly rooted in spectroscopy, as their contact [28-29] with data 
[13] reveals. Indeed the performance of this spectroscopy-oriented BSE-SDE framework 
in predicting the vacuum condensates, can be directly attributed to its off-shell structure. 

An inportant (new) aspect of this Study has been a demonstration of the powers of 
MYTP extending from its original mandate [15] of transversality in terms of 'Covaiant 
Inst ant aneity' (CIA) [16,17], to a wider 'transversality' on a Covariant Null Plane (CNPA) 
[41], thus vastly enhancing the applicability of this important Principle. In this article we 
have tried to present both CIA and CNPA on very similar lines, but the mathematical 
viability of the latter [41] seems to exceed that of the former [16], inasmuch as a CIA 
treatment of triangle (and higher) loop integrals is fraught with problems of 'Lorentz 
mismatch' of different CIA wave functions, leading to ill-defined integrals due to the pres- 
ence of time-like momentum components in the exponential/gaussian factors inside the 
integrals concerned [57]. This problem, which has been known since the FKR paper [25], 
is properly circumvented in CNPA, except for the (less serious) problem of dependence 
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on the 'null-plane orientation' which can be tackled through other means, e.g., a sim- 
ple device of 'Lorentz completion' which yields an explicitly Lorentz-invariant structure. 
This has been illustrated in Sect. 5 for the pion form factor which shows the expected 
high energy behaviour as well as very reasonable results [58-62] in both the high and low 
energy regimes. For more general three-hadron amplitudes [31] too, similar calculations 
in Sect. 6 show that the anomalies of ill-defined 4D loop integrals are absent in a CNPA 
treatment. The only exceptions are two-quark loops [32] (Sect. 7), where both methods, 
CIA and CNPA, work. 

Clearly, the MYTP is a very powerful Principle which helps organize a whole spectrum 
of phenomena under a single umbrella. It has been possible to study only a very few 
(though illustrative) examples to bring out its powers, but its potential is vast, and 
warrants many more of such applications. More importantly, the 3D-4D structure of BS 
dynamics provided by MYTP takes in its stride the spectroscopy sector as an integral 
part of the dynamics, as envisaged long ago by Feynman et al [25]. 

A good part of the logic behind this Article was evolved during my tenure of an INSA 
Professorship (1989-94), while the actual contents of this Article include both published 
and unpublished unpublished material developed subsequently, in my capacity as a free- 
lance worker (unattached to any Institution), as part of an ongoing research process. 



Appendix A: Derivation of F(k 2 ) and Nh for P-meson 

In this Appendix we outline the main steps to the derivation of the P-meson form factor 
(5.9), as well as the Normalizer (5.8), given in Sec. 5 of Text. Collecting the various pieces 
after p 2 n-pole integration, gives for (5.1) 

F(k 2 ) = 2{2TifN n {P)N n {P')m l J d 2 q ± dz 2 P.ng(z 2 )e~' i y^~^^ 2 + [1=>2]; (A.l) 

f(z 2 ) = M\f[9 k zl - z 2 k 2 m 2 + 6 k m 2 2 k 2 /4}; (A.2) 
D n + D' n = AP.n[ql + M 2 (z 2 2 - z 2 k 2 m 2 /2 + m 2 2 k 2 /4)/eta k - A/4M 2 ]; (A.3) 

D n + D' n M 2 + 5m 2 

9(Z2)= 4 WTWl + {Z2)] ( } 

h(z 2 ) = 2P.n(m 2 - z 2 )[M 2 - 5m 2 + m 2 M 2 (5m 2 — M 2 — k 2 /2)/(M 2 + k 2 /A)] (A.5) 

The integration over q± and z 2 are both routine, the latter with a translation z 2 — > 
z 2 + \m 2 l : 2 /9 k , to reduce the gaussian factor to the standard form. Note that, unlike the 
conventional (Weinberg) form [39a] of light-front dynamics, the present 4D form which 
permits off-shellness of the internal momenta, does not restrict in principle the limits of 
z 2 integration. Thus after the translation, the odd-z 2 terms can be dropped, and f(z 2 ) 
reduces to 

f(z 2 ) = M 2 z 2 2 9 k / V 2 + (Mm 2 k) 2 /A9 k (A.6) 
while the (/-function is a sum of two pieces g± + g 2 : 

9l = V k[qi + M 2 z 2 2 / Vk + - A M 2 m 2 2 k 2 {l + 3k 2 /4)/6 2 k - A/4M 2 ](1 + 5m 2 /M 2 ); (A.7) 

g 2 = 2 Vk (M 2 - 5m 2 )m 2 /6 k + 2(5m 2 — M 2 — k 2 /2)m 2 2 r ]k 2 /9 k (A.8) 
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Before writing the final result for F(k 2 ), it is instructive at this stage to infer the 
normalizer N H of the hadron, obtained by setting = 0, and demanding that -F(O) = 1. 
This gives after some routine steps: 

N n (P)~ 2 = 2M(2n) 3 (P.n/M) 2 J d z qe~^' p '' '^(O); (A.9) 

G(0) = [(1 + 5m 2 /M 2 )(q 2 - A/4M 2 ) + 2m x m 2 [M 2 - 5m 2 )] (A.10) 

where q = (q±,Mz 2 ) is effectively a 3-vector, in conformity with the requirements of 
the angular condition [35d,38], which gives a formal meaning to its third component q 3 = 
Mq.n/P.n = Mz 2 . The normalization factor N n (P) is also seen to vary inversely as 
P.n, while the multiplying integral is clearly independent of the NP-orientation n^. To 
exhibit this P n independence more explicitly, define a 'reduced normalizer' Nh which 
equals N n (P) x P.n/M and gives for Njj 1 the Lorentz-invariant result, eq.(5.8) of Text. 

Now insert the result N n (P)= MN H /P.n on the RHS of (A.l), and note, via eq.(5.3), 
that 

M 2 /overP.nP'.n = M 2 /overP.nn k ; n k = 1 - k 2 /A. (A.ll) 

One now checks that the factors P.n cancel out completely, and the evaluation of the 
gaussian integrals leads after a modest algebra to eq.(5.9) of Text, where G(k), after 
collecting from eqs.(A.6-8), is given by 

G(k) = (1 + 5m 2 /M 2 )h(k) + 2(M 2 - 5m 2 )m 2 /6 k + 2m 2 2 n k 6 k 1 (5m 2 — M 2 — k 2 /2); (A.12) 
h(k) = (1 + n k 2 /29 k )[3 2 - A/4M 2 + (Mm 2 k/2d k ) 2 (l + ^P); 5m = m x - m 2 . (A.13) 

Appendix B: Gauge Corrections to Kaon E.M. Mass 

We outline here a practical procedure to evaluate the gauge corrections to the e.m. self- 
energy of a qq system, vide Fig. lb of [18a]. For brevity we shall refer to the figures of KL 
[63] in their notation without drawing them anew. Thus Fig, lb of [18a] corresponds to 
fig la of KL [63], except for the presence of the hadron lines at the two ends. We shall 
call this simply 'la', with the understanding that the hadron lines are 'attached' to la. 
For the actual mathematical symbols (including phase conventions) we shall draw freely 
from [18a], without explanation. In [18a], only la of KL [63] was calculated, but now one 
must add 2(a,b,c,d,e) of KL [63], all with hadron lines understood at the two ends of each. 
There is no need to calculate lb or lc of [63] which are mere e.m. self-energies of single 
quarks (g.i. by themselves), and are routinely absorbed in quark mass renormalization 
(of little significance in this study which has these masses as inputs). 

A new ingredient is a 4-point vertex in each of 2(a,b,c,d), and two 4-point vertices in 
2e, except that the word 'point' is now understood as an extended structure characterized 
by the hadron-quark vertex function D(q)<f)(q) where one must insert a photon line in 
each such Hqq blob. Since it is not a standard point vertex, the method [63] of inserting 
exponential phase integrals with each current is not technically feasible; instead we may 
resort to the simple-minded substition — CiA[xi) for each 4-momentum pi (in a mixed 
p, x representation) occurring in the structure of the vertex function, which has the same 
physical content, at least up to first order in the e.m. field, without further comment. 
This amounts to replacing each q^ occurring in T(q) = D(q)(f)(q), by g M — e q A^, where 
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e q = m 2 e\ — m\e 2 . The net result in the first order in is a first order correction to 
T(q) of amount e q j(q).A defined by 

j(q).A = -4Mq.A<f>(q)(l - D(q) / (4M[3 2 )) (B.l) 

(The effect of the hat structure of q on the e.m. substitution is ignored in this approximate 
treatment). This effective 4-point vertex function is operative at one end in each of 
2a, 2b, 2c, 2d of KL [63] and at both ends of 2e. For the e.m. vertex at the quark lines of 
2(a,b,c,d), we use simply ie^.A, as in [18a]. The matrix elements can now be written 
down on exactly the same lines, and the same phase convention as in [18a] , to keep proper 
track of the gauge corrections with sign. We need write these down only for 2a and 2e, 
noting the equalities 2a=2b, as also 2c=2d, and the further substitutions (1) — > (2) and 
vice versa to generate 2c(=2d) from 2a(=2b). The contribution from 2a [63] to the e.m. 
quadratic self-energy of a kaon is expressible as 



Ml = N 2 H (27iy 5 e ieq J j{q)»D{q , )(t>{q')Tr[ lb D F ^{k) 

S F (pi - m 1 k)ie 1 -f 1/ S F (p' 1 )-f 5 S F (-p' 2 )]d 4 qd 4 k (B.2) 

where p[ = p 1 +rh 2 k and p 2 = p 2 = p 2 — rh 2 k are the 4-momenta of the quarks at the other 
(right-hand) end, and the photon propagator in the Landau gauge is —i(5 flu — k fJi k u /k 2 )/k 2 . 
It is now convenient to change the variable from k^ to q' noting that q' = q + rh 2 k, which 
gives d 4 k = d 4 q'/m 2 , etc. This shows that fig 2a(=2b), where the photon line ends on 
the heavier quark mi, gives a bigger contribution than does fig.2c(=2d) which would give 
m{ 4 arising from the (i 4 /c-measure. Evaluating the traces, and integrating over the poles 
of the two time-like momenta go and q' gives for the sum of the contributions from 2a-2d 
to the quadratic mass difference between K° and K~ as a product of two 3D quadratures 
after some simplifications with factorable approximations a la [18a]: 

2 _ 6N 2 H M5(e ieq ) r r 00' D(q) 

[(<f (2 - 4/tt) - qq'/3)(M 2 - 5m 2 + Dtf)^ 1 /2 + D^J' 1 /2) 

+-m 2 qq'(D(q')uo'f l /2 + M 2 — 5m 2 )] + [1 <-> 2] (B.3) 
3 

Here S(eie q ) is the K° minus K~ difference between the indicated charge factors associated 
with line 'i', while uj i2 2 = rn 2 2 + q' 2 and uj\ k = mf + (q — rhik) 2 . 

Next the contribution to 5M 2 arising from fig 2e of KL [63] which involves the product 
of two vertex blobs like (B.l) is given by 

5M 2 e = iN 2 H (27i)- 5 e q 2 J d 4 qd 4 kD F ^(k)j(q)^(q) u Tr[ l5 S F ( Pl - mi k) l5 S F (-p 2 + m 2 k)} 

' (B.4) 
This integral is somewhat different in structure from (B.2) in as much as k^ is fully 
decoupled from either wave function 0, 0', both of which have the same argument q. This 
makes it possible to integrate first over d A k as well as the time-like component q of q^ 
neither of which is involved in the vertex function. The relevant integral after tracing and 
rearranging has the form 

F(q) = 3H) 2 Jd 4 kjdq k- 2 (5, u -k,k u /k 2 ) 

[q 2 -q 2 + m x m 2 - m x m 2 {P - k) 2 \j{^ 2 ) (B.5) 
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where Aj = m 2 + (pi — rhik) 2 . The integral which is entirely convergent works out after 
some standard manipulations involving Feynman techniques as well as differentiation 
under integral signs as 

F(q) = 67r 3 [miTO 2 + f + A][VX - yf A - m x m 2 M' 2 \ / \m x m 2 M) 2 (B.6) 
where A = rhim 2 M 2 + D(q)/2M. And the final expression for (B.4) in terms of (B.6) is 

5M 2 2e = N 2 H (2ny 5 5(e 2 q ) J d s qj(q) 2 F(q) (B.7) 

Further evaluation of (B.3) and (B.7) can be made a la [18a] in a straightforward way. 
The key ingredients are 

5e ieq = 0.236e 2 ; 5e 2 e g = 0.139e 2 ; be\ = -0.0294e 2 . (B.8) 

The break-up of the final results for the diagrams 2(a-e) after dividing the results of(B.3) 
and (B.7) by 2M, since 5M 2 = 2M5M, is (in MeV): 

5M 2a+2b = -0.6996; SM 2c+2d = +0.1358; (B.9) 
5M 2e = -0.0481; SM tot = -0M2MeV. 

All these corrections, which reinforce one another due to a complex interplay of signs, 
add up to a figure which increases the value — 1.032MeV due to Fig. 1(b) of [18a], to 
— 1.644MeV , roughly a 60 percent (negative) increase, which is a rough indication of the 
type of QED gauge correction to be expected from such diagrams. 

Appendix C: A 4D NJL-Faddeev Model 

We summarize here the results of a simplified 4D NJL-Faddeev bound state problem 
[47,89], with 3 scalar-isoscalar quarks interacting pairwise in a contact fashion, a la NJL 
[4]. It is merely a special case of 3D-4D-BSE when its kernel K becomes a constant A. For 
ease of comparison, we employ the same notation and phase convention for the various 
quantities as in Sees. (4,9), but in view of the bound state nature of the problem it is 
enough to work with the 4D BSE for the wave function only, without invoking Green's 
functions. We start with the qq problem as a prerequisite for the solution of the qqq 
problem. 

C.l qq Bound State in NJL Model 

The BSE for the 4D wave function $ for a qq system may be written down for the NJL 
model: 

i(2n) 4 $(q 12 P 12 ) = Ar 1 A 2 - 1 A J d% 2 '^{q l2 'P 12 ) (C.l) 

where A is the strength of the contact NJL interaction for any pair of (scalar) quarks. 
The solution of this equation simply reads as [87b] 

$(912^12) = AAr 1 A 2 - 1 (C.2) 
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When plugged back into (C.I), one gets an 'eigenvalue' equation for the invariant mass 
M d = —P 2 2 of an isolated bound qq pair in the implicit form of a determining equation 
for A: 

A" 1 = -i(2n)- 4 J d 4 qA 1 ~ 1 A 2 - 1 = h(M d ) (C.3) 

where Ai )2 = m 2 + q 2 — M\ /4 ± q.P\ 2) and we have indicated the result of integration by a 
function h(Md) of the mass Md of the composite bound state (diquark). Unfortunately the 
integral (C.3) is logarithmically divergent, but it can be regularized with a 4D ultraviolet 
cut-off A, together with a Wick rotation, i.e., go - ► iqo, which is allowed by the singularities 
of the two propagators. The exact result is: 



J Am 2 - M3 

lQ^X- 1 = 1 + \n(4A 2 /M 2 ) - 2^ q - -axcsin(M d /2m,) (C.4) 

M d 

under the condition Md < 2m q . A slightly less accurate but much simpler form which 
is also easier to adapt to the qqq problem to follow, may be obtained by the Feynman 
method of introducing an auxiliary integration variable u (0 < u < 1) to combine the two 
propagators, followed by a Wick rotation and a translation to integrate over d A q (ignoring 
surface terms which formally arise due to the logarithmic divergence) : 

6A 2 

lG^A- 1 = In - - ^ 2 - 1 = 16n 2 h(M d ), (C.5) 
bra* — M d 

thus defining a diquark 'self-energy' function h(M) where the 'on-shell' value is M = Md- 
Eq.(C.5) also provides a determining equation for the NJL strength parameter A in terms 
of the 'diquark' mass Md and the cut-off parameter A, in a clearly 4D invariant form. 



C.2 NJL-qqq Bound State Problem 

We now set up the corresponding NJL-ggg problem under the same q — q contact inter- 
action strength A. Using the same notation for the various 4-momenta and propagators 
as listed in Sec. 9, the 4D wave function $(£,?7;P) expressed in terms of any of the S3 
invariant pairs (Ci,Vi) °f internal 4-momenta satisfies the BSE: 

*(2vr) 4 $(£, m P) = ^ AA 1 _1 A 2 _1 / d 4 q[ 2 <S>(&, V3 ; P) (C.6) 

123 J 

where the arguments of $ on the LHS are not-indexed since it is SVsymmetric as a whole, 
while those on the RHS are indexed in order to indicate which subsystem is in pairwise 
interaction (see explanation in Sec. 9). The solution of this equation may be read off from 
the observation that the integration w.r.t. q' l2 = ^£3/2 leaves the respective integrals as 
functions of rji only, where % = 1,2, 3. Thus [87b] 

^,r ] ;P) = J2^ 1 ^- 1 F(vs) (C.7) 

123 

where F is a function of a single variable r/j. Next, plugging back the solution (C.7) into 
the main equation (C.6), gives the following integral equation in a single variable i] 3 , as a 
routine procedure applicable to separable potentials [87b]: 

(h(M d ) - h(M 12 ))F(rj 3 ) = -^(27r)- 4 A 3 - 1 J d'q'^F^)^ 1 + (1 <- 2)] (C.8) 
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Note that the cut-off parameter A drops out from the LHS, as checked by substitution 
for h(M) from (C.5). This means that the 4D diquark propagator (h(M d ) — /i(M 12 ))^ 1 is 
formally independent of the cut-off A, in this simple NJL model. 

Next, the meaning of the function F(rj) can be inferred from an inspection of eq.(C.8), 
on similar lines to 3D [87b] or 4D [89] studies: F(i] 3 ) is the 4D 'quark(3)-diquark(12)' wave 
function which is generated by an exchange force represented by the propagators A'^ 1 and 
AfT 1 in the first and second terms on the RHS respectively. And the barjon-qqq vertex 
function V3 corresponding to a break-up of the baryon into quark(3) and diquark (12) may 
be identified by multiplying this quantity with the product of the inverse propagators of 
quark(3) and diquark(12): 



V 3 = V( m ) = A 3 f( m )F( m ) (C.9) 
where the diquark inverse propagator is reexpressed as 

/(») = MMj) - h(M 12 ) = (4.)- In 6m "^|:^ /9 . (CIO) 

making use of eq.(C.5) and the kinematical relation A, = m 2 q +rjf — M§/9, where M B is the 
mass of the bound qqq state, and i — 1, 2, 3. The quantity V3 of eq(C.9) may be compared 
directly (except for normalization) with the corresponding '3D-4D-BSE' quantity (9.3.10). 



C.3 Solution of the Bound qqq State Eq.(C.8) 

We now turn to the Lorentz structure of the NJL-ggg equation (C.8), as well as an 
approximate analytic solution for the energy eigenvalues of the bound qqq states. To that 
end we substitute (C.9) in (C.8) to give an integral equation for Vfa), with r)' 2 = rj for 
short: 

V( V3 ) = -2*(2vr)- 4 J dVfa)/- 1 ^) x (ml + v 2 - M%/%)-\m\ + ( m + V f - Ml/9)- 1 

(Cll) 

where the factor 2 on the RHS signifies equal effects of the two terms on the RHS of 
(C.8). For a bound state solution of this equation, with M B < M d + m q , the singularity 
structures permit a Wick rotation r] — > ii] which converts rj into a Euclidean variable t)e- 
This shows without further ado that eq.(C.ll) is 4D-invariant just like its qq counterpart 
eq.(C.3). This is not quite the same thing as the old result [13] on 0(4)-like spectra with 
harmonic confinement in the limit of infinite quark mass [14], since this NJL-Faddeev 
model of contact interaction, patterned after similar approaches [89], lacks a confining 
interaction, so that although in principle eq.(C.ll) predicts a spectrum of bound states at 
the qqq level (starting with NJL (contact) pairwise interactions), such spectra cannot be 
a realistic representation of the actual hadron spectra [13]. We now show how this comes 
about via Wick rotation in (Cll). 

For an approximate analytic solution of eq.(C.ll), note that the logarithmic function 
f(rj) in the integral appearing on the right is slowly varying, so that not much error is 
incurred by taking it out of the integral and replacing it with an average value < f(r)) >. 
It is now possible to 'match' both sides with an effectively constant V(i]), provided any 
further logarithmic dependence on rj is also similarly treated for consistency. The integral 
is now exactly of the type (C.3), i.e., logarithmically divergent, and can be handled 
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successively by Wick rotation, Feynman auxiliary variable u, and a translation. The 
result is again of the form (C.5), and after cancelling out the factors V(r] 3 ) and V(rf) from 
both sides, the eigenvalue equation reads: 

</ W > = 2(4,)-> <i?|/6>+ A ^_ M|/27 - 1 ] (C.12) 

To simplify this equation, we express all quantities in terms of the h(M) functions given in 
(C.5) and (C.10), and ignore the difference between 77 = r]' 2 and 773 inside the logarithms, 
to give 

h(M d ) - h(M 12 ) = 2h(M 12 ); A" 1 = h(M d ) = 3h(M 12 ) (C.13) 

The last equation brings out clearly the fact that the baryon binding comes about from 
three pairs of qq interaction, albeit off-shell, since the function < M\ 2 >=< n\ > — 
still depends on the (average) value of n 2 . The qualitative features are thus on expected 
lines, but this oversimplified model is not intended for a realistic fit to the nucleon/Delta 
masses (which at minimum require the introduction of spin-isospin d.o.f.), beyond the 
general feature of a quark-diquark structure that characterizes an NJL-Faddeev approach 
[92], as expected from any separable potential [87b], of which the NJL model is a special 
case. 



C.4 Comparison of NJL-Faddeev with 3D-4D-BSE 

We end this Appendix with a comparison between the vertex functions (10.3.10) and 
(C.9). The NJL-Faddeev form (C.9) of V3 is Lorentz invariant, being derived from a BSE 
with a constant kernel, viz the K = const limit of 3D-4D BSE [16]. Its quark-diquark form 
merely reflects the 'separable' nature of the NJL model [4]. There is no motivation here 
for a 3D BSE reduction, or 4D reconstruction, since 4D invariance is in-built throughout. 

In contrast, the vertex function (10.3.10), obtained from 3D-4D-BSE [47], is merely 
Lorentz covariant due to the 3D kernel support, but the derivation is otherwise more 
general than NJL-Faddeev, since it is valid for any spatial form of the kernel as long as 
it is 3D in content. This leads to an exact 3D reduction of the (4D) BSE whose formal 
solution is a 3D wave function 0(£, 17), a function of two independent 3-momenta [10b], in 
contrast to its NJL counterpart ^(773) in (C.9) which is a function of a single 4-momentum 
773 only. The denominator function D(qi 2 ) of (10.3.10) similarly is a 3D counterpart of the 
corresponding 4D inverse propagator 7(773 ) in (C.9). Finally the big radical in (10.3.10) 
corresponds to the inverse propagator A 3 in (C.9), except for its more involved structure, 
which we now seek to explain. 

While the 'zero extention' in the temporal direction is common to both approaches, 
NJL-Faddeev has also a zero spatial extension, while 3D-4D-BSE has a 'normal' spatial 
extension. Thus any ambiguity in the reconstruction of the 4D wave function from the 
3D form of the 3D-4D BSE, vanishes in the K = const limit, so that the same is directly 
attributable to the (mere Lorentz covariant) 3D form of the BSE kernel. Indeed the ID 5- 
function in (10.3.10) fills up an information gap in the reconstruction from a truncated 3D 
to the full 4D Hilbert space in the simplest possible manner, while satisfying a vital self- 
consistency check by reproducing the full structure (10.2.12) of the 3D BSE. This already 
lends sufficiency to the ansatz (10.3.4) which leads to (10.3.10). As to its 'necessity', 
this ansatz has certain desirable properties like on-shell propagation of the spectator in 
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between two successive interactions, as well as an explicit symmetry in the p 3 and p% 
momenta. There is a fair chance of its uniqueness within some general constraints, but 
this is still short of a formal 'proof of necessity'. 

The other question concerns the compatibility of the ID S- function in (10.3.10) with 
the standard requirement of connectedness [87]. Both the 5-function and the A 3 ^ prop- 
agator appear in rational forms in the 4D Green's function, eq. (10.3.7), reflecting a free 
on-shell propagation of the spectator between two vertex points. The square root feature 
in the baryon-ggg vertex function (10.3.10) is a technical artefact corresponding to an 
equal distribution of this singulariity between the initial and final state vertex points, and 
has no deeper significance. Furthermore, as the steps in Sec. 10.2 indicate, the three-body 
connectedness has already been achieved at the 3D level of reduction, so the 'physics' of 
this singularity, generated via eq. (10.3.4), must be traced to some mechanism other than 
a lack of connectedness [87] in the 3-body scattering amplitude. A plausible analogy is 
to a sort of (Fermi-like) 'pseudopotential' of the type employed to simulate the effect of 
chemical binding in the coherent scattering of neutrons from a hydrogen molecule in con- 
nection with the determination of the singlet n — p scattering length [92] . Such 5-function 
potentials have no deeper significance other than depicting the vast mismatch in the fre- 
quency scales of nuclear and molecular interactions. In the present case, the instantaniety 
in time of the pairwise interaction kernel in an otherwise 4D Hilbert space causes a similar 
mismatch, needing a ID ]delta- function to fill the gap. And just as the 'pseudo-potential' 
in the above example [92] does not have any observable effect, the singularity under rad- 
icals in (10.3.10) will not show up in any physical amplitude for hadronic transitions via 
quark loops, since the Green's functions (10.3.7) involve both the 5- function and the 
propagator A 3 p in rational forms before the relevant quark loop integrations over them 
are performed. 

Appendix D: 5(7(6)0 0(3) Wave Fns In Complex Basis 

In this Appendix, we outline a general method ofexpressing the qqq wave functions in a 
complex basis [29b. 46b], as an alternative to the 'real' representation given in eqs.(9.1-5). 
Such a basis gives a compact realization of the doublet representation of the permutation 
group S3, with the two complex vectors z, z* substituting for the real pair £, 77. The action 
of the permutations on this basis in the order (12); (31); (23) is [94a] 

P ijZ = [1; e 2iw/3 ; e~ 2i7T/3 }z*; P ijZ * = [1; e" 2m/3 ; e 2i7T /3]z (D.l) 

Identical doublet representations hold for the orbital ip, spin \ an d isospin d.o.f.'s, in 
the notation of Sections 9-11. To that end, define the corresponding complex quantities 
(except for an overall i-factor) 

v% c ; Xc ; C ] = W - x" - *x'; 0" - #1 (D.2) 

togrther with a second set of complex conjugate relations. Using these definitions, the 
action of the permutation group on the full wave function is 

Pa [fa ®Xc®<P; Vc ®x* c ® €} = Wc ®xl®€^c®Xc®<P] (D.3) 

Another important result concerns the action of Pij on any pair of component wave 
functions (A c , //*), where (A c ,/i c ) are any two out of the full set (-0,x, 0) of three [94a]: 

Pij[X c <S> /i*; A* <g> fx c ] = [A* <g> fi c ; A c <g> //*] (D .4) 
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As to the quantities (i/j s ;i/j a ), they are S^-singlets by themselves, with eigenvalues ±1 for 
the Pij operators. The properly symmetrized SU(6) x 0(3) states, eqs.(9.1-5), are now: 

|56 > q = ^ W; |56 > d = {xc€ + xM/V2; (D.5) 
1 70 > q = x s (^c€ + rM/V2; \70> d = (^ cXc c + fe^)/^ ( D -6) 



|20 > q = iiax 8 ^ |20 > d = MXc€ - XM/V2 (D.7) 



D.l Construction of ^-Fns in Complex Basis 

We now turn to the construction of the orbital -^-functions in terms of (zi, Z*), so as to 
preserve the total angular momentum adapted to the complex language. To that end, the 
angular momenta (both diagonal and 'mixed') in the complex basis are given by 

L z = -iz x V z ; L z * = +iz* x V*; L c = -iz* x V z ; L* = +iz* x V* (D.8) 

which obey the connections 

L = L z + L z . = Li£ + L^; L a = L a = L z + L z * (D.9) 

These quantities transform according to eq.(D.l) under the elements of S3. 

To construct angular momentum states of correct S3 symmetry, it is useful to take 
those of highest seniority [94b], now expressed in appropriate powers of z + and z + *, and to 
note that zbfz*, z + z + * and z+ 3 are all ^-invariant . The angular momenta carried by these 
basic units are easily checked to be in conformity with the above (complex) definitions 
(D.8-9) of the angular momenta. Using these basic building blocks, the natural parity 
states of highest seniority [94b] for a given angular momentum are compactly written in 
an HO basis as [46b]: 

|56 + ; 70"; 70+; 56" > = (2z+z+*) e [l; z + ; z 2 + ; zjje" 2 ""'* (D.10) 

The superscripts ± on the various states on the LHS correctly describe their parity struc- 
tures, by noting that z + n has parity (—1)™, while zz* is a 3-scalar. The L p - values of the 
states (D.10) in this order are [46b, 29b]: 

L p = (2£) + ; (21 + I)"; (2£ + 2) + ; (2£ + 3)"; (D.ll) 

while i goes through the values 0, 1,2,3..., thus bringing out the naturalness of the re- 
spective parity structures. 

In a similar way it is possible to systematically span all the "unnatural" parity states 
in the same representation [46b], noting that the main carrier of unnatural parity is the 
axial vector ( = iz x z*, which is a fully antisymmetric SVsinglet. The L p -structures of 
such states of highest seniority, corresponding to the series (D.10) are [46b, 29b]: 

|20 + ; 70"; 70 + ; 20" > = ( + (2z + z + *)% z + ; z\- z 3 + }e~ 2zz * (D.12) 

together with the respective L p -values 

L p = (2£ + 1)+; (2£ + 2)"; (2£ + 3) + ; (2£ + 4)~; (D.13) 

thus bringing out the 'unnaturalness' of their respective parities. For the construction of 
more involved states on these lines, see [46b]. 

A similar construction is possible for the 'spin' wave functions in the complex basis; 
see [46b] for details. 
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D.2 Normalization of Natural and Unnatural Parity Baryons 

We now outline a new method of integration for the normalization of the spatial wave 
functions (D.10) and (D.12) in the 6D (z,bfz*) space, which is rather well-suited to the 
(complex) variables on hand [29b,46b]. The volume measure in this 6D space is expressible 
in the spherical basis as 

d 6 r = d 3 zd 3 z* = (dz + dz*_)(dz^d* + )(dz 3 dz* 3 ) (D.14) 

where the six elements on the RHS of (D.14) have been rearranged into 3 sets of real 2D 
volumes, since the three pairs on the RHS, each form complex conjugate pairs. Now put 

V2(z + ;z*_)=n 1 e ± ^; V2(z_; z* + ) = H 2 e^- V2(z 3 ; z*) = K^ 3 ; (D.15) 

Then the volume element (D.14) becomes 

d e T = n 1 dn 1 d®{k 2 dn 2 d® 2 il z dll?,d'& z \ (D.16) 

< 7^1,2,3 < inf; < 1>2j3 < 2tt; 'R% + 'R% + 'R% = 2zz* = TZ 2 (D.17) 

Since the phase angles (not quite Euler angles) will not appear in the squared modulii of 
the wave functions, these are integrated out to give 

d T = it*dn\dn\dnl (d.is) 

Now the natural parity sequence (D.10) is compactly expressed as 

i, = N en (2z + z* + Yz + n e- n2 / 2 (D.19) 
where n — 0, 1, 2, 3 for the states (D.10) in sequence, and the normalizer is 

AT" 2 = J d\\n\ll\^\ll\l2) n e-^ = n 3 T(£ + 1)T(£ + n + l)/2 ra , (D.20) 

which agrees with the result for the £, ^-representation [29a]. For the unnatural parity 
sequence (D.12), the extra ^-factor gives 

(C + C) = ^ 4 /4 - nt/ A - K\K\ (D.21) 

Denoting the corresponding normalizers by N^ n , similar integration now leads to the result 

K 2 = " 3r(£ + ^| 2 ( ! ) +n+1) x [(£ + n + l)(n + 2) + (i + 1)(£ + 4)] (D.22) 

Radial excitations can be similarly handled. E.g., one radial excitation gives an extra 
multiplicative factor in the normalization integral (D.20) for natural parity states, giving 
rise to an extra factor (2£ + n + 4) in (D.20). Further, the reciprocity between the 
momentum and coordinate spaces implied in an HO description as above, allows the 
same formulation to be adapted in the dual space, a result which is useful for evaluating 
the OGE corrections to the mass formula (11.20). 
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